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Synopsis 


Combustion of fuel is a highly complex phenomena involving interactive heat, 
momentum and mass transfer processes both inside and outside each individ- 
ual droplet. It becomes even more complex when the fuel is multicomponent 
by nature. An understanding of the combustion process is very much needed 
for it’s wide industrial applications. The present study deals with the anal- 
ysis of a laminar multicomponent single droplet fuel burning in a quiescent 
atmosphere. The spray is assumed to be very dilute so that a droplet is not 
influenced by the presence of the neighbouring droplets. The effect of radi- 
ation, sorret and Dufour effects are neglected. The properties of both gas 

* 

and liquid phases are assumed to be constant in general. However, later on, 
for one of the cases variable specific heat for the gas phase was also consid- 
ered. There is no internal circulation or flow in the liquid droplet due to the 
absence of any external flow. The Raoult’s law and the Clausius-Clapeyron 
relations are used for coupling the liquid and gas phase concentrations at 
the liquid-vapour interface. Micro-explosion is neglected. The stationary 
coordinate system is transformed in such a way that the droplet surface is 
always maintained at r=l. This converts moving boundary problem into 
a fixed boundary problem. However, due to this transformation, a convec- 
tive term appears in the liquid phase heat and mass transfer equations. At 
the droplet center, radially symmetric boundary conditions are applied for 
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both heat and mass transfer equations. At the droplet surface (r=l), the 
overall mass conservation relation, gas-side flux balance relation for each of 
the species, liquid-side flux balance relation for each fuel, heat flux balance 
relation and the Clausius-Clapeyron relation for phase equilibrium at inter- 
face are applied for both boundary conditions as well as for retrieving the 
coupling parameters. The far boundary conditions for the solution domain 
are assumed to be at far stream conditions. The gas phase properties are 
assumed to be same as those of free stream conditions. The study is carried 
out in two parts. 

In the first part, semi-analytical solution has been developed for the 
spherico-symmetric burning of binary fuel droplets. The gas phase is as- 
sumed to be quasi-steady. Thin flame combustion, unit Lewis number ( for 
all gas species ) are taken to be valid. For rendering the liquid phase equa- 
tions to be solved semi-analytically, the convective term, which appears in 
the liquid phase transport equations due to the coordinate transformation 
and which takes into account the surface regression rate, is neglected. The 
highly non-linear reaction term is eliminated by using Schvab-Zeldovich vari- 
ables and then the quasi-steady gas phase equations are solved analytically 
for 1 < r < Tf and rj < r < r^a- Liquid phase transient energy and species 
equations are solved in infinite series form using 100 terms. The liquid phase 
and the gas phase equations are solved in a sequential manner through an 
iterative scheme. The liquid phase transient equations are solved analytically 



in space and numerically in time direction. Thus, the name ’semi-analytical’ 
appears. The semi-analytical solution is compared with the finite difference 
solution of the same governing differential equations and excellent agreement 
(upto 4 places of decimals) is observed between them. It was found that the 
semi-analytical solution was about 10 times faster than the numerical solu- 
tion. For the finite difference procedure, using a uniform mesh size, it had 
become necessary to use atleast a mesh size of 1/1000 within the droplet for 
achieving a good match witlj the analytical results. Ignition and extinction 
can not be analysed by this model. 

Results are preseirted for combinations of n-heptane, n-octane and do- 
decane fuels. While presenting the results, certain salient features are high- 
lighted. The time history of mass-fraction and temperature profiles of heptane- 
dodecane droplet for the liquid phase are presented. Transient variation of 
the important surface parameters such as mass-gasification rate, surface tem- 
perature, square of the diameter, fractional mass-gasification rate etc. are 
also looked into. These results show the exsistence of two <P — law periods 
as shown by experiments. A transition zone exists between the two cP — law 
periods. An initial transition also exists. This initial transient and the inter- 
mediate transition zone are characterised by vigorous droplet heating which 
can not be neglected. The effect of initial mass-fraction combination and 
volatility differential on the transient nature of multicomponent fuel burning 
is also presented. As the initial mass-fraction of the more volatile fuel is 
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gradually increased, the first (P — law period becomes more prominent and 
the second — law period gradually decreases. As the volatility differential 
between the two fuels is increased progressively, the transition zone becomes 
sharper, and the difference of the two <P — law periods becomes progressively 
more distinct. 

In the second part, full transient heat and mass transfer equations for both 
gas and liquid phases are solved numerically. Finite rate flame chemistry is 
used to analyse the nature of the flame. In addition to this, dissociation of 
products at high temperature near the flame zone is also included. Finite 
difference method is used to solve the liquid phase equations by implicit 
scheme in the time direction. Galerkin Weighted Residual Method is used 
to discretize the gas phase transport equations. Both liquid and gas phase 
discretized equations are solved by TDM algorithm. The gas phase solution 
for the present combustion system contains regions in which variables exhibit 
high spatial activity resulting in sharp peaks. For such cases, these high 
activity regions need to be resolved adaptively for ensuring efficient solution 
to the problem. Therefore, an adaptive grid generation technique is used. 
It is effected statistically by equidistributing a positive weight function in 
such a way that the local discretization errors are minimized. Both the gas 
and liquid phases are solved iteratively in time direction to trace out the 
entire burning history of the droplet. The present model can analyse both 
the transient phenomena namely ignition and extiction in detail with high 
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efficiency. 

The coupling parameters for the liquid and the gas phases are the evapo- 
ration velocity (itj), heat flux rate coming from the gas pha^e to liquid phase 
{Hi) and the fractional mass gasification rate (e,). Initially, these parameters 
are guessed. Based on these guess values, both liquid and gas phase solutions 
are obtained, from which new values of the coupling parameters are obtained 
and checked for convergence. 

Initially, at a given time, a coarse grid is chosen as a guess, and for this 
guess, a unique solution for the entire domain is reached by iteration. Then 
based on the solution obtained, total discretization error is checked against 
global tolerance limit. If it exceeds the latter, fresh grid is adapted according 
to the local error distributions. This loop, for a given time, is repeated until 
the optimum grid distribution is reached. Thus the solution proceeds in time 
direction. The adaptive grid algorithm is tested against a known analytical 
solution. It showed excellent agreement. 

In this part of the thesis, a full numerical model is developed for the case 
of a transient multicomponent single droplet burning in an otherwise quies- 
cent atmosphere. An adaptive grid generation technique is used based on the 
equidistribution of a positive weight function such that the local discretiza- 
tion error is minimised. This adaptive grid generation technique helps to 
record accurately the flame structure as well as the position and movement 
of the flame. This study gives a critical analysis of how a diffusion flame 
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around the liquid droplet develops during the ignition period, how it moves 
away from the droplet, how it moves towards the droplet and finally, how it 
is extinguished. The study reveals that the flame zone is initially established 
very near the droplet surface and then as the chemical reaction gathers mo- 
mentum, the flame temperature rises very fast and almost instantaneously 
reaches it’s peak temperature. Simultaneously, the flame moves and at the 
same time expands towards the droplet surface in search of fuel vapour for 
sustaining the enhanced chemical reaction at the flame. As the flame moves 
and expands towards the droplet surface, the temperature gradient at the 
surface becomes very high as a result of which the droplet surface is heated 
up very fast and the surface temperature approaches the boiling point of 
the more volatile fuel. So, the fuel starts evaporating vigorously. For this 
high evaporation rate, the flame starts moving away from the droplet surface 
due to fuel vapour accumulation. But as the burning proceeds, the droplet 
volume becomes progressively smaller for which the overall mass-gasification 
rate gradually falls. Therefore, fuel vapour accumulation between the droplet 
surface and the flame becomes progressively lesser. So, the flame then be- 
comes steady maintaining a fixed distance from the droplet centre. Finally, 
towards the end of burning, the flame moves back towards the droplet surface 
almost instantaneously. 
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The flame exists till the droplet becomes very small. The total burn- 
ing time for a heptane-hexadecane droplet for 0.9-0. 1 initial mass fraction 
combination is about 0.0658048 s. Ignition becomes complete within 0.003s 
and extinction takes about 260 pis. Results show that the pre- ignition pe- 
riod, initial transient ( ignition period ), first (P — law period, 2’^'^ transient 
( transition ), second <P — law period and extinction cover approximately 
2.28 %, 2.28 %, 35.6 %, 19.5 %, 40 % and 0.36 % respectively of the total 
burning time. This study brings out the two (P — law periods quite explicitly. 
The multicomponent feature of burning is predominant when the volatility 
differential of the two components is high, and also, when the initial liquid 
side mass fraction of the more volatile fuel is around 0.9. Results show that 
the flame is not thick. Comparison shows that the surface regression rate is 
neither too high nor too low. Therefore, we can neither neglect the effect of 
surface regression, nor can we assume the batch distillation model to be valid. 
Gas phase dissociation, taking gas phase specific heat to be variable and fuel 
vapour accumulation are all important considerations and they should not 
be ignored. Through parametric study of important parameters such as XI, 
D, (5p, Cp, it is found that the total burning time as well as other burning 
characteristics are very much sensitive to the property values being used. 
Comparison of the quasi-steady model proposed in this thesis, it is observed 
that the quasi-steady assumption of the gas phase is quite reasonable bar- 
ring the initial transient and the final extinction. The present theoretical 
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results shows very good agreement with the experimental results available in 
literature for hepane fuel burning under zero gravity condition with normal 
ambient conditions. The comparison with the experimental results available 
for the burning of heptane-hexadecane fuel with 0.7-0. 3 initial liquid side 
mass fraction combination, the agreement is quite satisfactory considering 
the fact that the present model deals with spherically symmetric combustion 
whereas, experiments were performed on freely falling droplets under natural 
convection conditions. These experiments show three periods, the two peri- 
ods correspond to the preferential gasification of the two fuels and the third 
period corresponds to the short sudden intermediate transition in between 
during which vigorous droplet heating takes place. The present model also 
exhibits three dominant zones, two (P — law periods separated by a transition 
zone which exhibit a gradual change. The present model does not include 
the elFect of thermal expansion for the liquid phase which may cause this 
deviation between the theoretical and experimental results. 
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1.1 Preface 

Combustion has become very important and relevant to the present day civilization, spe- 
cially when, the rate of energy consumption is growing alarmingly, day by day, due to 
progressive technological development. The modern day society can not survive even for a 
moment without energy. Life has become very fast. In order to meet this extremely high 
demand, a progerssively large variety of fuels are being used. So far, fossil fuels namely 
coal and petroleum have been used extensively, but because of relatively easy handling and 
transportation, the use of the latter i.e. petroleum has been on the increase. Therefore, 
research on gasification, finite rate flame chemistry, flame-movement tracking, dynamics 
of liquid and gas phases, dissociation due to high temperature of fuel droplets etc. have 
become very important as far as both practical and fundamental research in energy and 
combustion science is concerned. These studies have lead to modelling of practical com- 
bustion situations, and their verification with experiments which have revealed various 
burning characteristics. Desired objectives in the combustion research have been to estab- 
lish design criteria for efficient and stable combustors, determine the heat transfer rates 
to combustion chamber surfaces and examine the formation of pollutants such as soot, 
unburnt hydrocarbons, oxides of nitrogen and carbon. While many studies are purely ex- 
perimental, an underlying theme has always been to develop predictive models for spray 
combustion processes in order to reduce the cost of development by cut and try methods. 

Fundamentally, droplet combustion involves chemically reacting multicomponent two- 
phase flows with phase change. It involves complex physical and chemical phenomena, such 
as, heat, mass and momentum transfer processes in gas and liquid phases and their coupling 
at the interface. Kinetic effects are also important, specially, for critical phenomena such 
as ignition and extinction. In typical combustors, it is found that the gas temperatures are 
rather high and the gas is fax from a fuel- vapor-saturated state so that regression rates of 
the liquid-gas interface are high. , 

The fuels are sprayed into the combustion chamber through the injectors. Injectors are 



classified into two major categories (Faeth,1977): (1) pressure atomization, where only the 
liquid passes through the injector; and (2) twin fluid injection, where atomization of the 
liquid is accelerated by flow of a high velocity gas through the injector passage. 

Injectors take care that drop sizes are properly distributed. It’s properties are characterised 
by the angle of spray and by it’s ability to produce spray structures with full cones and 
hollow cones. 

In pressure atomization, the atomization is achieved through pressure drop in an orifice 
and in twin fluid injection system, one fluid is blasted through another. There are many 
steps in between leading to the final atomization. They are break-up of primary jet, break- 
up of secondary droplet and collision between droplets. These sequences lead to an evolving 
distribution of drop sizes and velocities which are normally measured some distance away 
from the injector. Spray characterization is generally carried out under cold conditions, 
but there are factors that influence this characterization. Thus, exact description of the 
velocity, shape and size of the fully atomized droplets become very difficult to estimate 
and we must be satisfied with average spray characteristics ( Mellor et el., 1970). 

The fuels may be either mixed with or kept separate from the oxidiser before it is 
injected into the combustion chamber. When they are thoroughly mixed before being in- 
jected, the resulting flame is known as a pre-mixed flame. It moves at a characteristic flame 
speed. This type of flame is observed in the stabilization of spray flames on flame holders, 
rocket engine combustion, and the development of deflagration and detonation waves in 
spray fields, Burgoyne and Cohen( 1954 ) inferred through their classical experiment on 
pre-mixed flame combustion that when drop sizes are less than 10 microns, the resulting 
flame speeds are not very diflferent from the pure gas case. But when the droplet size 
was increased, the flame speed decreased due to the time required to evaporate the fuel. 
Williams( 1965 ) agreed with this theory. Mizutani( 1973 ) pointed out that the burning 
velocities, at a fixed overall fuel-air ratio, are increased due to the presence of drops in the 
gaseous fuels. It also extends the region of stable burning to leaner overall mixture ratios. 
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A diffusion flame occurs when the fuel(s) and the oxidiser are kept separate before 
injection. After injection of the fuel, evaporation, diffusion of one species into another and 
chemical reaction between the fuel vapor and the oxidant give rise to a luminous diffusion 
flame around the droplets. Diffusion flame can exist as a part of pre*mixed spray in the 
sense that the spray is injected from a single point into an air stream and the combustion 
process proceeds simultaneously with the mixing process. This happens because diffusion 
flames can surround each droplet. Burke and Schuman( 1965 ) described the simplest type 
of diffusion flame. For such a case, the concentration of all species becomes zero at the 
flame. This enables the flame position to be determined as a function of distance from the 
source( Williams, 1965 ). Bilger( 1976 ) dealt with more practical situation where reaction 
between fuel and oxidiser did not take place at a fixed radial position rather, there was a 
broad region of overlap of the fuel and the oxidiser concentrations. The diffusion flame is 
automatically taken care of depending upon the availability of relative fuel and oxidiser. 
For such cases, a lot of pollutants, such as oxides of nitrogen, axe produced due to very high 
reaction rate and high flame temperature. For pre-mixed flame, this control can be achieved 
by controlling the mixing of fuel and air. Thus, generation of high temperature at flame, 
and correspondingly, formation of hazardous oxides can be avoided for a premixed flame, 
though this mixing needs a high level of sophistication. But the pre-mixing causes danger 
of undesirable burning and explosion which makes this mode of combustion practically 
still unviable and in most of the industrial combustors, diffusion flame combustion mode 
is used. 

1.2 Literature Survey 

Combustion phenomenon is a very complex process involving a very high degree of non- 
linearity. So, initially, the study of the combustion process started with simplified models 
like quasi-steady burning of a single component, individual fuel droplet burning in a qui- 
escent atmosphere in absence of any external flow. Later, in due course of time, complex- 
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ities in the modelling were gradually incorporated in order to simulate realistic situations. 
These are; external convection, droplet heating, transient nature, internal convection, fuel 
vapour accumulation, finite rate flame chemistry, variable property, micro-explosion, ef- 
fect of neighbouring droplets, presence of more than one fuel in the liquid phase etc. The 
present literature survey mainly deals with combustion in the diflfusion flame mode relevant 
to the present study. 

1.2.1 Single Component Fuel 

(1) Spherically Symmetric Burning 

Agafanova et el.( 1958 ), Godsave( 1949, 1953, 1950, 1951 ), Goldsmith and Penner( 
1954 ), Spalding( 1950-, 1953 ), Varshavskii( 1945, 1957, 1962 ), Wise et el. ( 1955 ) and 
Williams( 1972 ) are the pioneers who first formulated the basic droplet combustion model 
for an isolated, single-component, quasi-steady droplet burning in a stagnant, oxidising 
environment. This basic model has become famous by the name ’d-square law model’. 

It predicts that the surface regression rate of a droplet is proportional to the burning 
time elapsed. In this model, the chemical reaction rate is assumed to be infinitely fast 
and thermal diffusion & radiation losses are neglected. This has been shown to be the 
case except in the case of heavy fuel oils ( Probert, 1946 ; Friedman and Churchill, 1965; 
Wolfhard, and Parker, 1947, 1949 ). But later many authors questioned the validity of the 
above <P—law model. Krier and Wronkiewicz(1972) showed a general relation (d”) between ; 
surface regression rate and the time elapsed. This relation is based on their experimental 
observation that the flame stand-off ratio ( ratio of flame radius to droplet radius ) grows | 
linearly with time. The value of n lies between 1 and 2 when the chemical reaction rate | 
becomes the controlling factor and it equals 2 when diffusion becomes the rate controlling 
process. 

(i)Quasi-Steady Combustion and Infinite Rate Kinetics: 

Williams( 1960 ) assumed quasi-steadiness, which implies that soon after the ignition | 
period, the burning reaches steady state condition i.e. the plot of (P against time remains 
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linear after the initial transient. The quasi-steadiness also assumes that the liquid side tem- 
perature becomes spatially uniform and invariant with time. This temperature was taken 
to be the boiling point of the fuel. He also assumed the droplet radius to be fixed. This 
means that the velocity of the fuel vapour flowing away from the surface is much greater 
than the surface regression velocity i.e. the effect of surface regression was neglected. The 
general continuity equation for predicting the burning rate coefficient ( defined as the area 
of the droplet surface regressed per unit time i.e. (d^))) temperature and composition 
field surrounding a burning droplet were developed by Hirschfelder et el.( 1954 ). Most 
analyses are based on constant transport properties. Wise and Agoston( 1972 ) gave a 
typical and widely accepted solution based on temperature independent transport proper- 
ties. Solutions depended upon the choice of gas side specific heat, Cp and conductivity, A. 
Both conductivity and specific heat were generally calculated at the arithmatic mean or 
the logarithmic mean of flame temperature and droplet surface temperature. 

In practical situations, both temperature and the composition of the gas mixture vary 
over wide ranges. More sophisticated methods have been used by Sioni and Roblee( 1969 
) and Annamalai et el.( 1971 ). Sioni and Roblee used a computer calculation of A while, 
Annamalai et el. used a more accurate analytic type calculation for evaluating conductivity, 
A. Temperature dependent expressions for A and Cp are also used. Goldsmith and Penner( 
1954 ) and Faeth and Lazar(1971) assumed linear variations of A and Cp of the gas mixture 
with temperature. Kassoy and Williams ( 1968a, b ) and Fendell( 1966 ) have also developed 
analytic solutions based on singular perturbation technique, allowing for a general transport 
property variation with temperature and variable Lewis number and have shown their 
effects on flame temperature, droplet radius and burning rates. They were obtained in 
terms of asymptotic expansions and are given for the limit of the Damkohler number, 
which is defined as the ratio of chemical reaction rate to heat conduction rate, approaching 
infinity. Analytic forms are available for simple temperature dependence( 1968b ) but for 
more complex forms, numerical methods have to be used. An alternative diffusion theory of 
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combustion applicable to the special case of combustion of a liquid hydrogen droplet in still 
air was given by Varshavskii and Germeier( 1967 ). Brzustowski( 1965, 1966 ) developed 
a dimensionless form of the quasi-steady droplet combustion theory. Raghunandan and 
Mukunda ( 1977 ) in their qutisi-steady analysis considered the variation of specific heat 
and thermal conductivity as a function of temperature and concentration of various species. 
(ii)Finite Rate Kinetics 

Very precise theoretical indications of the rate of chemical reaction rates in droplet 
burning are very difficult to quantify. Brzustowski( 1965 ) pointed out the general restric- 
tions imposed by chemical reaction and also the evaporation rates. Agafanova et al.( 1958 
), Lorell et aJ.( 1956 ), Tarifa et al.( 1962 ), Williams( 1961 ), and others 
used a simplified chemical reaction mechanism for their combustion analysis involving fi- 
nite rate kinetics. For such cases, an overall reaction rate expression is adopted. It is seen 
that using such an ( overall reaction rate ) expression, for typical hydrocarbon fuel burning 
under normal conditions, little error is introduced in the mass burning rate relationship by 
assuming infinite reaction rate except under certain extreme conditions. Inclusion of finite 
rate kinetics gives a more realistic information regarding the structure of the flame zone 
or the reaction zone in the sense that it lies nearer to the droplet surface compared to that 
indicated by the infinite reaction rate, though the evaporation rate hardly changes. This 
inclusion also brings down the flame temperature to a more realistic value. Agafanova( 
1958 ) first made this realistic consideration which is supported by experimental results [ 
Kumagai and co-workers ( 1957, 1971, 1975) , Wang et al. (1984) etc. ]. These effects are 
highlighted by Lorell et al. ( 1956 ) in their calculation of flame structure. Their numerical 
results show that the mass burning rate and the location of the flame front are hardly 
affected by variation in the activation energy (hence, the reaction rate) of the chemical re- 
action. Bracco( 1973 ) has also used this type of analysis to determine the rate of formation 
of NO in hydrocarbon droplet flames. Chung and Law(1984) worked with constant, but 
mutually different gas properties for the two regions across the flame for their quasi-steady 
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model. They considered finite rate chemical kinetics and has shown that the gas-phase 
dissociation results in a lower flame temperature. These predictions are in conformity with 
the earlier analyses of Lorell et al.(1956) and Agafanova et al.(1957). 

(iii) Transient Combustion 

Experimental results indicate that the heat and mass transfer involved in and around 
the droplet during it’s combustion are in a transient state during the major part of the 
droplet lifetime. Kumagai and Isoda( 1957a ) were the persons who first tried to solve 
the transient energy equation. Their analysis considers only gas-phase unsteadiness while 
neglecting radial convection terms. Spalding( 1959 ), Chervinsky( 1969a ), Kotake and 
Okazaki( 1969 ), and Krier and Wronkiewicz( 1972 ) carried out further analyses on solving 
transient equations. The transient theory of Spalding(1959) makes use of a point source 
at the origin where the fuel is injected into the field to simulate the burning drop. The 
radial convective velocity is assumed to be zero. By considering the gas-phase unsteadiness, 
Spalding ( 1959 ) has shown that the flame initially moves away from and later towards 
the droplet. The transient analysis as well as the results of Chervinsky(1969) constitute 
a modification of Spalding’s work. Williams( 1960 ) discussed the effects of unsteadiness 
during the pre-ignition period and established the region of validity of the quasi-steady 
assumption. Parks et al.( 1966 ) gave an account of temperature distribution within the 
droplet, while Shyn et al.( 1972 ) investigated the liquid phase cracking in a heavy fuel oil 
during combustion. In all the above cases, time-dependent terms were retained while the 
analyses were carried out. Kotake and Okazaki( 1969 ) solved an unsteady state model by 
numerically integrating the equations. However, the spherical symmetry was retained. The 
equations of spcies, momentum and energy of the gas surrounding the droplet and that of 
energy within the droplet were written as the time-dependent differential equations and 
integrated. Hubbard et al.(1975) pointed out that gas-phase unsteadiness is insignificant 
in the evaporation process. This gas-phase unsteadiness exists for a very short time. 
Waldman(1975) and Crespo and Linan(1975) used singular perturbation technique to study 
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the far-field unsteadiness of the g«is phase for unsteady state combustion. Crespo and Linan 
used the ratio of densities of the ambient gas and the liquid as the small perturbation 
parameter, while Waldmau used the ratio of droplet radius to the characteristic diffusion 
field length as the small perturbation parameter. Both reveal that the flame initially moves 
away from the droplet surface and then comes back towards the droplet. This feature was 
also reported by Kumagai and Isoda(1957) and Kumagai et el.(1971). 

(2) Moving Droplets 

(i) Low Reynolds Number Flow; 

It is seen that in almost all practical situations, droplet burning takes place in an 
external velocity field. Calculation of the mass burning rate of a droplet in a convective 
field from fundamental principles would require quite an elaborate numerical scheme. 

The first realistic model of combustion in a convective field was put forward by Brzus- 
towski and Natarajan(1966). Essentially, Natarajan and Brzustowski(1966), and Brzus- 
towski and Natarajaii(1970) put forward a two-film model. This model assumes a spherical 
inner stagnant film through which diffusion and conduction occur, while the convective and 
diffusive mass transfers of the oxidant occur through the outer film. Due to the presence 
of convection, the oxidiser mass-flux towards the inner film is enhanced and is expressed 
empirically through dimensionless parameters. It is reported that the augmentation, for 
free convection, is proportional to ( Gr )? and for forced convection, it was reported to vary 
as ( Rea ). The evaporation rate was evaluated in terms of droplet surface heat transfer 
Nusselt number. A more detailed and a more fundamental analysis of droplet combustion 
in a slow viscous flow was carried out by Fendell et al.(1966). They assumed that the 
entire flow field was in stokes regime, the internal circulation was zero and the flame was 
thin. They solved the governing differential equations by singular perturbation technique. 
They used small values of droplet Reynolds number based on far-stream conditions as the 
perturbation parameter. The enhancement of burning rate and the distortion of flame due 
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to external flow field were evaluated. Sangiovanni(1978) used an integral formulation to 
see the effect of droplet motion. He calculated the thickness of the convective film around 
the droplet ( within which both heat and ma.ss transfer were assumed to be important ) by 
using Ranz and Marshall’s (1952) correlation. He used the Green’s function technique to 
evaluate droplet burning rate and flame stand-off ratio. Law and Williams(1972) tried to 
link the earlier theoretical and experimental analyses. Okazima and Kumagai(1975,1982) 
conducted a series of experiments to see the burning of both stationary and moving droplets 
under normal and micro-gravity conditions. It was reported that burning rate enhance- 
ment due to droplet motion was insignificant. But, in the micro-gravity experiments, the 
free-fall height was inadequate. The effect of internal circulation inside the liquid droplet 
in the creeping flow regime was looked into by Ayyaswamy and co-workers(1983,1986, 
1988,1990). Here the effect of radial flow was considered to be dominating. From another 
study of Sadhal and Ayyaswami(1983), it was observed that the internal circulation within 
the droplet causes the coldest region to be shifted from the droplet centre to a ring of 
radius approximately equal to two-third of the droplet radius. This reduces the character- 
istic diffusion length, but high rate of mass transfer(i.e. high surface regression rate) at the 
surface reduces internal circulation and may, eventually, cause a reversal of flow. Gogos 
and Ayyaswamy(1988) cited in their variable property analysis that a reduction of viscosity 
at higher temperature may also lower the internal circulation. Gogos et al.(1986) showed 
that the total drag comes from momentum transfer to the droplet due to evaporation. In 
this work, they accounted for the initial droplet heat up period. 

(ii) Intermediate Reynold Number Flow: 

For the droplet motion in the intermediate Reynold’s number zone, various numerical 
schemes are used to analyse the typical burning of liquid fuels. This is because in this 
Reynolds number zone, neither a boundary layer analysis ( generally applied for high Re. 
flow ) nor an asymptotic analysis (generally employed for creeping flow regime) can be ap- 
plied. So, one has to depend only on various eleborate numerical schemes. The works in this 
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area are, in general, quite recent. Dwyer and Sanders(1984,1986) considered the influence 
of drag on the deceleration of the droplet while applying the finite volume stream- vorticity 
function formulation. They took into account the effect of transients in both the phases 
and used non-orthogonaJ grids to track down the motion of the shrinking droplet. The 
drag results deviated from experimental data. Renksizbulnt and Ynen(1983) used steady 
state correlations and variable thermo-physical properties in their finite difference analysis 
to calculate drag. Their results indicated close agreement with experiment. Dwyer(1989) 
analysed an isolated droplet by finite volume method. He used variable property and 
continuity equation to correct pressure through a direct solver. For solving the transport 
equations, a predicator-corrector algorithm of the Alternating Direction Implicit ( ADI ) 
sweeping type was used. They did not measure drag. Recently, Renksizbulnt and Hay- 
wood(1988), Haywood et al.(1989) and Huang and Ayyaswamy(1990) assumed gas-phase 
quasi-steadiness in their variable property analysis for an evaporating drop. A Lagrangian 
frame of reference was used for balancing the forces on the droplet. The gas-phase analysis 
was coupled with the droplet motion. Huang and Ayyaswamy considered the effect of 
gravitational field. They adapted a hybrid solution scheme like that of Gogos et al.(1986). 
The SIMPLEC solution procedure developed by Vandoormaal and Raithby(1984) was used 
by Haywood and co-workers. All these studies of droplet evaporation belonging to the in- 
termediate Reynolds number zone point out that the evaporation drag contributes to only 
a small fraction of the total drag. This is quite contrasting to the low Reynolds flow 
situation. 

(iii) High Reynolds Number Flow: 

It is seen that boundary layer approximations are generally used for studying high 
Reynolds number(of the order of 100) flow situations. Prakash and Sirignano(1978, 1980), 
Tong and Sirignano( 1982a, b) used the boundary layer approximations for their study in 
high Re. flow regime. Prakash and Sirignano established the existence of boundary layer on 
either side of the droplet surface with the help of their order of magnitude analysis. They 



pointed out that the droplet core resembled Hill spherical vortex and that the gas-phase 
boundary layer merged with the potential flow. Tong and Sirignano(1983) developed a 
simplified model for the vaporization of a single component droplet. The model accounts for 
the liquid-phase internal circulation and axi-symmetric gas-phaise convection. The solutions 
have been obtained by coupling the gas-phase similarity solution to the simplified liquid- 
phase solution. Gakkhar and Prakash(1990) varied the ambient conditions and studied 
their effects on burning. 

In high Reynolds number flow situations, both for single component and multicompo- 
nent fuel droplets, it is noted that the characteristic diffusion length inside the droplet 
is reduced due to internal circulation. And at the same time, the droplet lifetime is also 
reduced by a similar scale due to enhanced evaporation rate. So, inspite of enhanced 
transport rate due to droplet motion, the transient heating remains almost throughout the 
droplet lifetime. 

(3) Interaction Between Spray Droplets 

In practical situations, it is found that there axe many droplets present in a spray system 
and the average distance between droplets can become as low as few droplet diameters. 
So, a typical droplet, in real situation, will not behave as an isolated droplet, rather it 
will be very much influenced by neighbouring droplets. Basically, there are two kinds of 
models used for describing the interaction of spray droplets. The first one is the group 
combustion model introduced by Suzuki and Chiu (1971), and the alternative approach is 
the discrete droplet burning model or droplet array model. In the latter, a few droplets 
in a well defined geometry or a large number of droplets in a periodic configuration are 
examined. The droplet group combustion theory is distinct from the array theory in that 
a statistical description of droplet spacing, rather than a precise geometrical description, 
is used. As a matter of fact, group theory can deal with many more droplets than array 
theory. Array theory gives a detailed field analysis i.e. it gives a more detailed picture of the 
burning characteristics in the vicinity of a typical droplet whereas, an overall averaging 
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or collective burning effects of a cloud of droplets is given more stress in the statistical 
approach of group combustion model. Results of the array theory can be employed as an 
input to the group combustion approach. The two approaches have relative advantages 
and disadvantages over each other depending upon different practical situations. 

Sirignano (1986) discussed different types of formulations of the governing equations 
for modelling spray combustion situations. In the two-continua approach where both the 
gas and the liquid phases are eissumed to be co-existing, group combustion model is suit- 
able. He formulated the governing equations in Eulerian frame, Lagrangian frame and 
discussed their relative advantages and disadvantages. For discrete droplet burning model 
or the droplet array model, Eulerian description can be used for solving the gas phase, 
and both Lagrangian description or Eulerian description can used for solving the liquid 
phase. Eulerian- Eulerian form is not advantageous when resolution is desired on scale as 
small as or smaller than the average distance between the droplets. In such situations, 
Eulerian-Lagrangian form is advantageous. Eulerian form is simple and easy to apply, 
but the Lagrangian frame can eliminate numerical diffusion, resolve multivalued solution 
and has the ability to achieve resolution on a scale smaller than the droplet spacing. He 
also outlined a probabilistic (or, distribution) function approach to analyse spray system. 
When the scale of resolution is even smaller or finer than to the average droplet spacing, 
the practical situation becomes probabilistic rather than deterministic, since it is difficult 
to locate any droplet. In such situations, probabilistic function approach becomes more 
realistic for modelling. 

(i) Discrete Droplet Burning: 

In this model, each droplet is considered as a fuel source as well as a heat sink. Twardus 
and Brzustowski( 1977 ), Brzustowski et al. (1979 ) and Umemura et al. (1981 ) employed 
a bi-spherical coordinate system to facilitate their analysis of the interaction of two burning 
droplets. Twardus and Brzustowski ( 1977 ) did not consider Stefan convection and forced 
convection. They indicated that two droplets may burn separately only when the ratio of 



the distance between droplet centers to the droplet radius is above a certain critical value, 
otherwise they burn with one single envelope flame. Chigier ( 1976 ) ha.s cited experimental 
evidence that the droplet in combustors rarely burn in an isolated fashion. Rather, a flame 
envelopes many droplets. Umemura et al.(1981) used a potential flow approximation and 
Brzustowski et al. (1979) used a Nusselt number correlation to account for the convective 
effect. All of them noted the decrease of burning rate due to the influence of neighbouring 
droplets. Umemura et al.(1981) saw that when the two droplets are of different sizes, the 
smaller droplet is affected more. 

Labowsky(1976, 1978, 1980a, 1980b) and Labowsky and Rosner(1978) used a modified 
superposition principle for solving a situation having simple arrangements of droplet array. 
Labowskk(1976) used this modified image method to determine vapour density field in 
arrays of spherically symmetric and isothermal evaporating droplets. In this study, Stefan 
flow was neglected which enabled the equations to be linear and consequently, method of 
images or superposition could be applied. Later, Labowsky(1978) included Stefan flow in 
his analysis, but the equations were reduced to Laplace equations through a suitable trans- 
formation. After this, method of superposition could be applied to solve the transformed 
linear equations. Labowsky ( 1980b ) extended his work (1978) to the case of interact- 
ing droplet combustion. He(1980a) also relaxed the isothermal condition of the droplet 
in his analysis. In the superposition method, an infinite succession of sources and sinks 
are placed at the image points of the earlier generations and thus, the flow pattern was 
simulated. Labowsky(1976) proposed an extrapolation formula which enabled to approxi- 
mate the boundary conditions within bsources. Labowsky and Rosner (1978) pointed out 
that the difference in results between a continuous distribution of monopole sources and 
the discrete distribution of monopole sources becomes insignificant when the number of 
droplets in the cloud is large. Reduction of burning rate was noted and this became more 

severe when the number of interacting droplets increased in the system. 

Ray and Davis(1980, 1984) have presented a generalised treatment of heat and mass 
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transport between an assemblage of particles and a surrounding continuum. This analysis 
is based on the solution of unsteady state diffusive transport equations with unknown point 
sources located at the centers of the particles. Integral equations mathching the boundary 
conditions at the surface of the particles are used to describe these point sources. They 
compared their results with Labowsky(1978) quite successfully. This formulation was ex- 
tended by Merberry et al.(1984) to the case of interacting droplet combustion. The theory 
is applied to burning droplets of different sizes as well as monodisperse droplets in multi- 
particle arrays. In their quasi-steady analysis, they developed a simpler solution technique 
by using a modified Laplace equation with point sources. They provided correction factor 
by which, burning rate for multiple particle can be retrieved from the burning rate of a 
single particle. Sangiovanni and Labowsky(1982) found that with decrease in droplet spac- 
ing, the burning time of monodisperse droplets in a linear array increased gradually. They 
initiated the quasi-steady halo model to simulate the transient nature of droplet interaction 
and suggested only the droplet having overlapping halos could interact directly. 

Miyaska and Law ( 1981 ) experimented on linear arrays of two and three droplets in 
low pressure environments. Their, both natural convection and low pressure effects were 
shown to be negligible. They indicated that the droplet vaporization rate and droplet 
heating were significantly retarded due to the interaction effect. Their result showed that 
the superposition calculation for a pair of droplets situated far apart were in qualitative 
agreement with this study. But when droplets were close to each other, the results deviated 
significantly. This mismatch was attributed to the longer heat-up period resulting from 
droplet interaction. Ryan et al.(1990) utilised the group combustion theory which involved 
a large number of droplets, assumed top hat profiles in the interstitial space between the 
drops, and accounted for their variation in the radial direction. Their study accounted for 
the change in transport rate due to proximity of other droplets and found close agreement 
between the results obtained by group combustion approach and droplet array study. Sam- 
son et al.(1978), in their statistical model, established the concept of a screening distance 



within the radius of which, the interaction of that droplet was significant. 

(ii) Cell Model Approach: 

Unit cell model approach due to Happel{1958) and Kuwabara(1959) is yet another 
powerful technique applied for the hydrodynamic analysis of a multiparticle system. In 
this approach, a given spray system is assumed to consist of a number of cells which are 
isolated from each other through suitable boundary conditions at the cell surface. Each 
of these cells has one droplet at it’s centre. Thus, by studying one representative cell, 
one can greatly simplify analysis of practical spray situations. Happel (1958) applied zero 
shear stress at the cell surface to isolate the cell from the surrounding. Kuwabara (1959) 
used zero vorticity condition at the cell boundary. Jaiswal (1990) used Happel model in his 
finite element analysis. In order to study the transient effect, Zung(1967) applied zero mass 
transfer condition in order to isolate the cell, in his pure meiss-transfer analysis. Tishkoff 
(1979) considered the heat transfer aspect in addition to mass transfer, and applied zero 
heat flux condition at the cell boundary. Bellan and Cuff el (1983) used spherical cells 
which were touching each other. The global equations were modified to accommodate the 
heat and mass transfer that take place in between the cells. Bellan and Harstad (1987a) 
used a correlation to account for droplet motion and later they( 1987b) modified the global 
conservation equations to account for convection. El-Kaissy and Homsy (1973) used regular 
perturbation technique while using Happel and Kuwabara models for their analysis. Tal 
and co-workers (1982, 1983, 1984a, 1984b) used cylindrical cell model to take advantage 
of the periodic structure of the array. The finite-difference calculations gave excellent 
agreement with experimental results for drag-coefficients, indicating that drag increased as 
streamwise and transverse spacings decreased in identical fashion. Tsai and Sterling (1990) 
used the grid-embedding technique to solve for velocity and temperature field in their 
cylindrical cell model study. Chen and Tong (1988) predicted stronger influence of droplet 
interaction for high Reynolds number flow in their three drops-in-ceU arrangement study. 
Whereas, Shuen (1988) in his planner droplet arrays oriented normal to the flow, indicated 
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that interaction effect became less with increasing Reynolds number. Ramachandran et 
al.(1987) gave a numerical model to study the transient characteristics of a train of three 
vertically falling droplets. These results comply with their earlier(1985) experimental data. 
Ramachandran et al.(1989) used three spheres-in-cell model in their finite element analysis 
to see heat-transfer and flow efi’ects. Wang (1989) used the above analysis to see interaction 
effects on vaporizing droplets. Recently, Sanyal and Sundararajan (1991) applied this unit 
cell approach in their study of quasi-steady, laminar, single component fuel combustion. 
They considered the droplet to burn in a conceptual spherical cell. They studied this 
spherical cell by regular perturbation technique in the creeping flow regime and by a full 
numerical solution using FEM in the intermediate Reynolds number regime. They assumed 
the cells to be touching each other. 

(iii) Group Combustion Model: 

Suzuki and Chiu(1971) proposed the first theoretical model on group combustion, fol- 
lowed by experimental studies conducted by many researchers such as McCreath and 
Chigier (1973), Hayshi and Kumagai (1975), Mizutani et al.(1977) and Khalil and Whitelaw 
(1977). Suzuki and Chiu (1971) established that when the group combustion number G 
exceeds a critical value, group combustion mode starts. Chiu and Liu (1977) and Chiu et 
al.(1983) used the characteristic group combustion number in order to isolate the two basic 
burning models, namely, the group combustion model and the discrete droplet burning(or, 
droplet array ) model. Nusselt number correlations were used to account for the droplet 
motion. Depending upon the value of the group combustion number (G), the burning of 
droplets was categorised as follows: 

Discrete droplet burning for G < 0.01 
Internal group combustion for 0.01 < G < 1.0 
External group combustion for 1 < G < 100 
External sheath combustion for G > 100 
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Annamalai and Ramalingam (1987) presented the analysis of a spherically symmetric 
cloud of particles burning in a quiescent atmosphere in order to study the group combustion 
behaviour. They predicted results for burning rate and the flame structure as function of 
G. They found that for char/carbon particles, group flame occurs i.e. the change from 
discrete droplet burning to the group combustion mode takes place at G > 5 while for 
a cloud of liquid drops, this occurs at G > 0.1. So, it is noted that this critical value 
of G for the liquid drops is one order lower compared to that for the carbon particles. 
This happens due to the lower burning rate of carbon particles inside the cloud. They 
established the existence of individual particle combustion, group combustion and sheath 
combustion. Later, they (1988) performed this study in cylindrical clouds. Lotters(1989) 
has indicated that the group combustion structure broke down due to various transport 
processes. 

1.2.2 Multicomponent Fuel 

It is observed in the practical situations that the fuels used are mostly multicomponent in 
nature. These fuel constituents have widely different properties. This variation of charac- 
teristics within the fuel alters the combustion characteristics as compared to that of the 
pure liquids. So, in order to simulate the practical situations for fuel burning in a more 
realistic manner, the study of multicomponent fuel burning is very essential. Therefore, in 
recent years, significant attention is being given to the burning of multicomponent fuels. 
Multicomponent droplet burning is more complex than single component fuel combustion 
due to mass-transfer effects in the liquid phase. The blend of the multicomponent fuel can 
be of two types: miscible and immiscible. The phase change characteristics of a multi- 
component fuel depend upon the extent of miscibility of the constituents. Two simplified 
models corresponding to batch distillation and mass diffusion linaited mode of burning 
have been developed to look into the multicomponent aspect of combustion. For diffusion 
limit model, in absence of internal circulation, diffusion is considered to be the only mode 
of transport inside the liquid droplet. For the batch distillation mode of burning, at any 



instant, the most volatile component has the fastest gasification rate relative to it’s con- 
centration such that the depletion of the constituents from the droplet composition occurs 
approximately in the sequence of their relative volatilities. The batch distillation mode of 
burning was initially looked into by Faeth ( 1970 ), Newbold and Amundsion ( 1973 ) and 
Law ( 1976 ). They adopted the batch distillation limit concept and assumed the droplet 
temperature and composition to be perpetually uniform such that they varied with time 
only. They did this in order to incorporate the effects of internal circulation. 

Based on the diffusion limit model, by numerically integrating the conservation equa- 
tions, it was shown by Law ( 1978a ), and Landis and Mills ( 1974 ) that the evaporation 
process consisted of an initial transient regime, an intermediate quasi-steady regime and 
a final volatility-dominated regime. This behaviour was attributed to the high diffusional 
resistance of the liquid phase. 

Law and Law (1982) examined the diffusion-limited quasi-steady regime with the as- 
sumption that no changes occured in the temperature or the mass-fraction profiles on the 
liquid side. Using an asymptotic analysis, they showed that the fraction of mass-flux leav- 
ing the droplet surface for any species was equal to the initial liquid mass-fraction of that 
species prior to evaporation. Shaw (1990) extended the study of Law and Law (1982), and 
found that the mass-flux fraction of a species at droplet surface was equal to the volume- 
averaged mass-fraction within the droplet at that time. Wang et al.(1984) and Niooka and 
Sato(1986) reported the occurence of sudden and rapid heating of the droplet after the 
constant surface regression rate was over. At this point in time, flame moves back towards 
the droplet surface and the burning rate falls. 

Considering a trace amount of low volatility component in the liquid phase, Shaw and 
WilBams (1990) developed an asymptotic analysis for heat and mass transport during the 
burning of a multicomponent fuel drop. For the entire process, evaporation of the impurity 
material was neglected. It was indicated by Makino and Law (1988) that a.s the diameter 
of the shrinking droplet approached the characteristic mass-diffusion length, the volatility 
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differential of the constituents became the rate controlling factor and thus the vaporization 
followed the batch-distillation mode. Law et al.(1977) analysed vaporization of a multi- 
component droplet in a strong convective gas stream. They demonstrated that boundary 
layer existed in both liquid and gas phases near the droplet surface over most of the droplet 
lifetime, and that convection was an important transport mechanism within the droplet. 
Lara-Urbaneja and Sirignano (1981) found that droplet heating effects were important for 
most part of the droplet life time , especially if the boiling points of the components dif- 
fered widely. They applied similar boundary layer analysis. Beisically, they followed the 
same arguments used in Prakash and Sirignano’s single-component droplet study (1980) 
and extended the analysis to the multicomponent droplet case. The gasification and the 
burning of the individual components are strongly influenced by the coupled transport pro- 
cesses occuring in the liquid and gas phases. The studies of prakash and Sirignano ( 1978, 
1980 ) and Lara-Urbaneja and Sirignano ( 1981 )have clarified that when the gasification 
rate is not too slow, rapid internal circulation can not cause perpetual spatial uniformity 
throughout the droplet, and therefore, the batch distillation limit can never be approached. 
The cause is that though the internal circulation is infinite and therefore having perpet- 
ual streamwise uniformity, transport in the direction normal to the stream lines is still 
effected through the slow diffusion process. They have shown that since the reduction in 
characteristic diffusion and burning time are of the same order even for a highly convective 
case, the transient mass diffusion should prevail throughout the droplet life time. Tong 
and Sirignano(1986) studied transient vaporization of a multicomponent droplet in a hot 
convective environment. They accounted for liquid-phase internal circulation, transient 
droplet diffusion, and axi-symmetric gas-phase convection in their model. This simplified 
model is an extension of their previous single-component droplet model(1983). 

Very recently, a fully transient numerical study for both the liquid and gas phases was 
presented by Mawid and Aggarwal ( 1991 ) for spherically symmetric multicomponent fuel 
droplet burning. The authors observed that for certain time intervals, the mass evapo- 
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ration from the droplet surface was suppressed. They attributed such a trend to vapour 
accumulation in the gas phase and the high diffusional resistance of the liquid phase. 

Mi cro- Explosi on : 

It may happen during the course of fuel burning that the droplet surface becomes more 
concentrated with the less volatile ( high boiling point ) component while, the droplet core 
has a higher concentration of the more volatile ( lower boiling point ) components. Thus 
the more volatile component in the droplet interior may accumulate substantial amount 
of superheat beyond it’s boiling point. This may start homogeneous nucleation, whose ex- 
tremely fast gasification rate may lead to sufficient pressure build-up and thus the droplet 
may eventually fragment, causing what is known as micro-explosion. This phenomenon can 
best be explained by diffusion-limit model ( Law, 1978 ). They accounted for liquid-phase 
internal e environment. They accounted for liquid-phase internal 5:qmicro-explosion can be 
controlled properly, then large-scale mixing can be acheived through spraying and penetra- 
tion, with some what larger droplets, which is then followed by instantaneous gasification 
and local mixing through micro-explosion. Law ( 1978 ) and Law( 1981 ) pointed out the 
possibility of increasing the occurence of micro-explosion by increasing the system pressure. 
This is attractive for fuel burning in high pressure environments. Law( 1981a ) pointed out 
that high pressure may also decrease the possibility of micro-explosion because, at higher 
pressure, the dynamics of bubble expansion is weakened. Higher pressure also increases 
internal motion. Thus, the concentration gradient is reduced, reducing the possibility of 
nucleation. Law et al.( 1979 ) predicted that micro-explosion is favoured in an optimum 
range of mixture composition. The reason for this is that the more volatile component 
facilitates homogeneous nucleation while, the less volatile component raises the droplet 
temperature. This was later experimentally verified by Lasheras et al.( 1980,1981 ). It was 
experimentally established by Lasheras et al.( 1980 ) that the intensity of micro-explosion 
for miscible multicomponent droplets was less. 

For the case of immiscible fuel undergoing equilibrium evaporation, the partial pressures 
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of the individual components are independent of the pressure of the other components. 
Therefore, the surface temperature of the droplet is limited to the boiling point of the 
more volatile fuel. This was first indicated by Law(1981). By making use of the above 
property, the surface temperature of the fuel can be suppressed by adding a small quantity 
of the low boiling point liquid like water. This can reduce pollutant formation as well as 
incendiary ignition. By adding a small quantity (a few percent) of surfactants to normally 
immiscible water and oil, a stable emulsion can be established. Cook and Law(1978) 
showed that the addition of water generally does not exceed 20 percent for smooth engine 
operations. Water-oil emulsions have been tested in different combustors of diesel engines, 
gas turbines, and furnaces and boilers ( Dryer, 1977 ). Wheatherford et al.( 1979 ) and 
Law( 1981 ) has come up with the fact that water/diesel emulsions are fire-safe because 
upon spillage and incendiary ignition, the initial fire ball will extinguish itself. Lasheras et 
al.( 1979 ), in their experiment, have shown that the emulsion- droplets explode with much 
more intensity than those of miscible fuel blends. This takes place due to the presence of 
larger amount of superheated mass contained by individual micro-droplets of water when 
the nucleation starts. Lee and Law(1991) have shown that the particles agglomerate into 
a rigid porous sphere during initial d-square law period in case of gasification for slurry 
droplets. This was an experimental verification of their own proposition. 

Experimental observations by Wang et al.(1984) indicate that two <P-law regimes occur 
during the burning of binary fuel mixtures. After the initial heating, the droplet follows 
cf^-law burning for the more volatile fuel until it is nearly depleted. Then rapid droplet 
heating takes place upto the boiling point of the less volatile fuel, leading to the d^-law 
burning regime of this fuel. 

1.2.3 Adaptive Grid Generation 

In Modelling environments or situations, where regions of high spatial activity are a func- 
tion of time, initial boundary value problems mixed with moving adaptive grids are needed 
to be solved. Several approaches are available. Miller and Miller(1981) used a moving finite 
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element :iniethod. In their method, they expanded the solution in piecewise linear polyno- 
mials vk- ere both time dependent coefficients of the series and the grid points are unknown 
functianss of time. These quantities were obtained by minimising the partial differential 
equati >oiu ’s residue in a least square sense. David and Flaherty (1982), in their adaptive 
finite •element method, did not couple together the calculation of the solution and the grid 
points; , but kept them seperate. Different approaches or methods were used by Dwyer et 
a], (190^, Tscharnuter and Winkler(1979) and Bolstad(1982) to solve the above problem. 
Wbite:(l f982), in his study, calculated the solution and the grid together by using finite dif- 
ferexcxe method. He generalised the idea of equidistributing the arc-length of the solution 
of a Wco. point boundary value problem to one-dimensional mixed initial-boundary value 
problens.s, Pereyra and Sewell(1975) equi distributed the local truncation error. Apart 
from iBcO'lstad (1982), all have used methods that use a fixed number of grid-points so that 
they concentrated in the regions where the spatial activity of the solution is highest. 
Essen.' tisilly, here, positive weight functions are equidistributed over a fixed number of spa- 
tial iciteirvals and the positive weight functions are chosen in such a way that it reduces 
the IqcsiI discretization error. Either finite element or finite difference approach is used 
to dicre'tize space in these cases. The weight function may depend on the solution either 
impli*«:it,Il 7 or explicitly. In case of implicit dependence, the solution components and the 
grid are; calculated together, where as for explicit dependence, the grid is determined hy 
using a previously calculated solution. Later, Smooke and Koszykowski(1986) developed a 
fully a.l.aptive method for solving one-dimensional mixed initial boundary value problem. 
Their ronethod had the capability to change the total number of grid points according to 
then, eeod, during adaption. In addition, when a fixed number of grid points were enough to 
accurafc-ely describe a solution, the mesh was moved by extrapolating the nodal positions 
from esiilier time levels. 
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1.3 Scope of Present Study 

The present study deals with the analysis of a laminar mtalticomponeiit sirngle droplet 
fuel burning in a stationary atmosphere. The spary is assumed to be v&rp Jliilute so that 
the effect of neighbouring droplets on the droplet under consideration is n •egCtgibly small. 
Liquid phase mass diffusion is considered to be the only means for Ii«(iiaii<l. phase mass 
transport. Liquid phase and gas phase properties are considered to be eo»:iis’tJtmts. Effect 
of radiation is neglected. Equilibrium condition is assumed to be always p trewveiling at the 
liquid-gas interface. Therefore, the Raoults law and Clausius Clapeyron lellftl tion hold true 
at the interface. At the droplet surface, overall mass balance, species flu>cb<)»,l!aa.ace for both 
liquid and gas sides and heat flux balance are used as boundary conditi-onnsHjtc couple the 
'two phases. At the center, radial symmetry is applied as boundary conditriioDo. Far stream 
conditions are assumed to be valid at the outer boundary of the solution <li.sMGia.in { r — > oo 
). Chapter 2 deals with the above general formulation part where general 'aeq[uations are 
derived in context of the present problem. The study is carried out in t'urooj-sarts. 

In the first part, a semi-analytical solution is sought. The additional asaanraptions made 
are the gas phase quasi-steadiness and unity gas phase Lewis nunnber. The effilect of surface 
regression which appears as the convective term in the liquid phase equa-tiooKS is neglected. 
Thin film approximation as well as Schvab-Zeldovich variables are used 41o eliminate the 
non-linear reaction term. A detailed parametric study is attempted to briiu.5croTit the salient 
features of the multicomponent fuel burning. The details of the forniu!a»..tiOD» and results 
and discussion for this first part are presented in chapter 3. 

In the second part, the full transient heat and mass transport eq^uaitioors are solved 
numerically for both the liquid and the gas phases. Finite rate chetaiisbiry, gas phase 
dissociation, use of variable gas phase specific heat are considered to ca-idte finer details. 
Unity Lewis number for the gas phase has been relaxed- The liquid equations are 

solved by using FDM and the gas phase equations are solved by using KEiM. Fixed grid 
distribution with unequal grid spacing has been used for the liqnid phasfs; and adaptive 

t. 
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grid generation technique has been developed to track down the steep flame region in the 
gas phase. The detailed formulation for this second part is included in chapter 4 and the 
details of the results and discussion are presented in chapter 5. 



Chapter 2 

General Formulation For 
Multicomponent Combustion 
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2.1 Introduction 

Combustion is a highly complex phenomena involving interactive momentum, heat, and 
mass transport processes both inside and outside of each individual fuel droplet. The 
complexity of analysis is further enhanced when the fuel droplet is of a multicomponent 
nature. Very high temperature, smallness of each fuel droplet, unpredictable flame chem- 
istry, complex gas-liquid phase coupling, flame movement, highly transient nature of the 
process are some of the key features of the combustion phenomena. Moreover the effect 
of the neighbouring droplets is almost unpredictable. This causes problem in specifying 
proper boundary conditions. Relative velocity between the droplet and the gas adds to 
the complexity. In view of such complexities, the complete modelhng of spray combustion 
becomes extremely involved and is the subject matter of the present chapter. 

2.2 Liquid Phase and Gas Phase Transport 

As soon as the atomised fuel droplets are injected into the combustion chamber, the liquid 
fuel immediately starts evaporating due to a concentration gradient at the droplet surface. 
The fuel vapour meets the oxidiser in stoichiometric proportion somewhere near the droplet 
surface, where ignition takes place. Ignition can be achieved by an ignition stimulus or if the 
environment is sufl&ciently hot. At the flame, fuels and oxidiser undergo chemical reaction, 
as a result of which, intense heat is produced. So, immediately, heat starts flowing from 
the gas phase to the liquid phase, thereby, effecting rise in droplet temperature and more 
evaporation of fuel at the surface. Due to the volatility-differential between the species 
of the liquid droplet, one fuel is preferencially evaporated than the other. This develops 
concentration gradient within the droplet and, the gradient is maintained throughout the 
droplet lifetime. Therefore, both heat and mass transport within the droplet are very 
important for multicomponent fuel droplet burning. Momentum transport inside the liquid 
droplet is also important because, as soon as the atomised droplets are introduced into the 
combustion chamber, interfacial shear stress is generated due to the relative motion between 





the droplets and the gas phase which creates circulatory motion inside the droplet. But 
the strength of this circulation depends on the Reynolds Number. If the relative velocity 
between the liquid and gas phase is low, the liquid phase circulation is negligibly small. 
So, for accurately predicting combustion behaviour, mass, momentum and heat transport 
phenomena generally become important for the liquid phase. 

As a result of chemical reactions at the flame, a number of products like oxides of carbon 
and ntrogen, water vapour and also many other intermediate species are generated. Infact, 
there are innumerable steps in between, before the final products are produced. But, a 
detailed modelling, which includes all the intermediate steps, is impossible because of the 
complexity of the system as well as due to the lack of complete knowledge of most of the 
intermediate reactions. In most of the combustion studies, a single step global irreversible 
reaction is generally used for a given fuel. Dissociation of the products at high temperature 
is another feature which influences the droplet burning. These dissociation reactions are 
endothermic in nature. So, it does not allow the flame temperature to rise beyond a certain 
limit. It is seen that unless the prevailing temperature is quite high, i.e. RuT > activation 
energy {Ea) of the fuel, the molecules cannot take part in chemical reaction. Therefore, 
the combustion reactions remain frozen in most of the places of the gas phase. Only in 
a very narrow zone, the molecules are generally, able to penetrate this high activation 
energy barrier and, therefore, flame exists only in this narrow zone. In this region, reaction 
kinetics becomes very important and the reaction proceeds to equilibrium at a finite rate. 
Finite rate flame chemistry is very important in analysing the structure of the flame. 

Assumptions of co-existing media are assumed to be valid for deducing the conservation 
equations for a multicomponent system ( both reacting and non-reacting ). In this concept, 
the entire volume is assumed to be accessible to all the components in the mixture. The 
various chemical species in a diffusing mixture move at diffrent velocities. Let denote 
the velocity of the species with respect to stationary coordinate axes. So, for a mixture 
of N species, the mass-averaged velocity is defined as 
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!> =!;&>; ( 2 . 1 ) 

t=l 

where, mass-fraction of the species. 

Diffusion velocity, V^, of the species with respect to the mass-average velocity of the 
mixture (u,) is defined as 

Y.i=m-v (2.2) 

The diffusion velocity, Vi, is evaluated from Pick’s law of binary diffusion as 


36 = 


D^Yi 

Yi 


(2.3) 


where, D is the coefficient of binary diffusion. In the above expression, the effect of temper- 
ature and pressure gradients as well as that of external forced fields axe neglected because, 
their effect on diffusion velocity are very small compared to the concentration gradient. 
The diffusion coefficients for liquids are much smaller compared to the low pressure gases 
owing to the densely packed liquid molecules. Pick’s law of diffusion can be assumed to 
be valid for multiphase systems also, if the binary diffusion coefficients of all the pairs of 
species are equal. This is valid if the molecular weights of all the species are close to each 
other. This assumption is generally valid for most of the combustion situations. Pick’s law 
of diffusion can also be used for multiphase systems in situations where dilute reactants are 
present in an almost inert environment. The above equations and relations are generally 
valid for both gas and liquid phases. 


2.3 Derivation of General Conservation Equations 

2.3.1 Species Conservation Equations 

The species conservation equation for each species, Y,, can be derived by making a species 
balance [Kuo, 1986] over an arbitrary differential fluid element of volume Ax Ay A^ in a 
mixture as 
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or, in vector notation, 


dpi dihix dthiy drhi^ _ . 

dt ^ dx dy dz 


dpi 

dt 


+ V-^i 




where, p is the mixture density, />,• is the density of the species, 

^ = rate of change of mass per unit volume of species i, 
an^ + = net rate of change of mass-flux vector in the 3 directions, 

rhia:, TTiiy and m ,-2 are the rectangular components of mass flux vector, rh^ and 
d>,- is the rate of production of the species i in kg/m^s. 

Using, Pi = pYi, Vi = Zi -I- u. 


(2.4) 


(2.5) 


m,- = piVi 

= py;(Z, + u) (2.6) 

Using pi = pYi and Eqs. 2.5, 2.6, one can write the general conservation equation for the 
i**' species as, 

+ V- [pYi QU + £)] = (2.7) 

(i) for liquid phase, the species conservation equations can be written as, 

+ 2 - [pY (£; + «)]= 4 ( 2 . 8 ) 

(ii) for gas phase, the species conservation equations can be written as, 

+ 2 [pYi (Vs + !l)l = Wi (2.9) 

The symbol, * stands for liquid phase quantities, and those without this symbol stand for 
the gas phase quantities. 
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2.3.2 Overall Mass Conservation Equations 

Adding the conservation Eq. 2.5 for all the N species, the general overall mass conservation 
is obtained as, 

^ do ^ N 

i=i i=l ,=i 

or, 

a N N 

^E^^ + v-E^. = o 

•=! t=l 

or, 

^ + v-m = o ( 2 . 10 ) 

(i) for liquid phase, the continuity equation can be written as, 

dp 

^ + V-ffi = 0 (2.11) 

(ii) for gas phase, the continuity equation can be written as, 

= 0 ( 2 . 12 ) 

2.3.3 Momentum Conservation Equations 


Momentum conservation equations can be derived by making momentum balance over a 
control volume as follows ( Kuo, 1986 ) 


A = B-\-C-VD 


(2.13) 


where, A = rate of increase of momentum inside the control volume ( c.v. ), 

B = Net influx of momentum into the control volume through control surface ( c.s. ), 
C = Net body force on control volume and 
D = Net surface force on control surface. 

So, the momentum balance can be expressed as in integral form as 


dt 


= - j j + j J J i + J J g-iA 

= -j j j xdv + J J g_.dA 

= ~ J J J (2 /^) + III xdv + 



or, 


III ^ (V-/=>VV) - V-£ - ^ 


dv = 0 


where, 


Therefore, 


or. 


N 


^ = pJ2^iL 


i=l 


_a 

dt 


N 

(pv) = - (v-pvv) + V-£ + P E YiLi 


t=l 


d 


0 ^ {pv) + (v-/ 0 Vv) = V-£ + P E >^7.- (2.14) 

t=l 

where, i? is the velocity vector, p is the density, g, is the total stress tensor, is the 
contribution due to the body force on the species and, N is the number of species. 
Here, 

V-£ = - VP + V- [v (w + (vn)^) - I (v-li) l] 

Therefore, Eq. 2.14 can be written as 


dt 


(psi) + (v-p w) = />$ + 3!^ + E (s. «) + «' 2 E 


(2.15) 


Since, 


Therefore, 


+ (v-PP.) = 0 ( from continuity eq. 2.12) 


p^ + pv-^v = -VP + V- \n - 1 (v-n) i] + /’E 


or. 


P^ + Pi^-W = -^P + 2- [»? (sii + (V2l)^) - 3 (v-li) l] - PS. (2-16) 

where, g is the gravitational force which acts on the system as body force, r) is the dynamic 
viscosity and, X is the unit tensor. For constant p and rj, the momentum conservation 
becomes 

^ ^ ^ (2.17) 
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(i) for liquid phase, the general momentum conservation equations can be written as, 

+■ pv.^v = -^p + + (5;u) ^) - I (s-i) Z] - Pg (2-18) 

(ii) for gas phase, the general momentum conservation equations can be written as, 

P^ + pv-Vv = -VP + V- [»? (vn + (V^) - I (v-v) Z] - Pg (2-19) 

2.3.4 Energy Conservation Equations 

Starting with the stationary differential volume through which a pure fluid is flowing, the 
conservation of energy for the fluid contained in the differential volume, can be written as 
[Kuo, 1986] 

A = B + C + D + E (2.20) 

where, A = rate of accumulation of internal and kinetic energy, 

B = net rate of influx of internal and kinetic energy by convection, 

C = net rate of heat addition due to heat flux by conduction, 

D = rate at which heat is added by the heat source and, 

E = net rate of work done on the system by the surrounding, which includes the work 
done by body force and the total surface stress tensor (£) acting on the boundary of the 
control volume. So, the energy balance can be written down mathematically as 

o . N 

(pe) = —2- (P^Sl) — V-g + Q Ag: Vjl + P^ (2.21) 

dt ,=i 

where, e is the specific internal energy, p is the mixture density of the fluid, v_ is the mass- 

average velocity, q is the heat flux due to conduction, Q is the rate of heat generation, £ 

is the total stress tensor, Yi is mass-fraction of the i^^species, £. is the external body force 

on the species and, V^- is the mass-diffusion velocity of the species. Here, specific 

heat, Cp of the overall mixture at constant pressure is defined as 

Cp - Efei Cp,iYi and standard heat of formation, for the overall gas mixture is defined 
as h° = Eili where, Cp,,- and ft? are the specific heat at constant pressure and the 
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standard heat of formation at the reference temperature, T° respectively for the individual 
species. In general, the Mach number is well below unity for most of the combustion 
situations. Therefore, the flow work contribution (^) arising in the energy balance can be 
neglected. Therefore, 



( 2 . 22 ) 


= 0 as there is no heat generation. Contributions of viscous dissipation and pressure 
work are also zero unless the Reynold’s number for the flow is very high i.e. £ : VH = 0- 
The external body force (/.) is only due to gravity. 

Therefore, 

N N 


p E = -Pill = -PI 


and since, 


YiVi = Yi{^-v), 


So, 


Now, 


1=1 i=l t=l 


1=1 


N 

«=1 


= 0 


£ = - + pY,^iYiYj 

i=l 

= - \VT + pfyiYu{hU ^ 


(2.23) 
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The first term in Eq. 2.23 represents the heat flow due to conduction and the second term 
corresponds to the total enthalpy of all the species per unit area per unit time flowing 
relative to the mixture. This energy is carried across a surface when the average velocity 
of the species differs from the mass average velocity of the mixture. Using the Eqs. 2.22 
and 2.23 and all the above conditions, the Eq. 2.21 for energy balance can be written as 

T 


dt 

or, 


+ Jv 


Cr,dT 


+^. 


+ jC ('“” + /! Cr.fiT] - Ayr 


= 0 


dt 


pj! + S' U C C,dT + r C,,dT - \^T 

J L t=:l 


d \ ^ 


t=i 


V- 


i=l i=l 


N 


t=l 




(2.24) 


Using species Eq. 2.9, the above equation can be written as 

rT N fT 


Cp<£rJ + 2- [«£ c,dT ■¥pY.y<y^ C,.idT - AyrJ =-■£ ( 2 . 25 ) 

Using Cp = I^ili Cp^iYi, the following simplification can be made: 


3.j^c,dT = v/Ji;c,,wr 

f — 1 

= jCs. (r dT + 

= [f:sy^£cv.>dT+c,vT 

or, 

J2-^yiJ^Cp^i<iT = -^J^^CpdT + CpVT (2.26) 

using Fick’s law of diffusion (eq.2.3), 

pt,i^Yij^C,^dT = -pDs.l^C^-hpDC,^T (2.27) 
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Substituting the above expression (Eq. 2.27) in Eq. 2.25 yields the general equation for 
the conservation of energy as 


dt 




+ {pDCj>-\)^T 



t=i 

(2.28) 


which for liquid and gas phases can be written as 
(i) liquid phase 


dt 




+v- 


n (/^ c^dtj - pbvijl 


N 

'Lh% 

’(2.29) 


(ii) gas phase 


dt 


» / Cj>dT 

Jt° 


+V- 


pv 


IjK) 


CpdT - pD^ / CpdT + {pDCp - A) vT 


lT° 


N 

i=l 

(2.30) 


2.4 Modeling for Multicomponent Fuel Droplet Com- 
bustion 


For the present problem, it is assumed that a single multicomponent fuel droplet is burning 
in a hot quiscent atmosphere. The only flow that is generated in the gas phase is due to the 
evaporation velocity in the radial direction that arises due to the evaporation of fuels from 
the droplet surface. Since, there is no transverse flow in the gas phase, the shear stress at 
droplet surface is zero. Therefore, momentum transport within the droplet is completely 
absent. So, for the present study of multicomponent droplet burning, only mass transfer 
and heat transfer aspects are considered. 

2.4.1 Liquid Phase 

(i)Species conservation equations 

Substituting the expression for diffusion velocity (F,) from Pick’s law of diffusion, Eq. 
2.3 into Eq. 2.7 and since there is no generation or consumption of species inside the liquid 



droplet i.e. Gji=0 and also the mass average velocity, 0 = 0 within the droplet, the species, 
conservation equations for the liquid phase can be expressed as, 


S {pY) 

dt 


V- [pD^\) = 0 


(2.31) 


where i stands for the two species of fuels fl(fuell) and f2(fuel2). 

The above equation for fuel 1 and fuel 2 for spherically symmetric case, will be read as 


= 0 

dt r2 dr V dr ) 


(2.32a) 


dt r2 dr V dr j 

(ii) Energy conservation equation 
The properties including Cp for the liquid phase are constant. Therefore, the liquid phase 
conservation Eq. 2.29 can be rewritten by, first evaluating the definite integrals and then 
by simplifying the algebra as 



Dividing both sides by pCp, and writing for the spherically symmetric case. 


(2.33) 


dT ad 


dt r2 dr 1 dr 




(2.34) 


2.4.2 Gas Phase 

(i) Species conservation equations 

The species conservation Eq. 2.9 for the gas phase can be simplified in a manner similar 
to that adopted for the liquid phase, and by using the constant property assumption, it 
can be written as 


dipY^) 


dt 


'■^+v-{pvyi-pDvy^ = ui 


(2.35) 



for spherically symmetric case, 


dt dr \ ' dr j p 


(2.36) 


where, i stands for fuell(/l), fuel2(/2) and the oxidiser(a:) for the gas phase. Replacing i 
by /I, /2 and x in Eq. 2.36 yields 


dY,, 1 9 

dt dr 


+ r‘^VrYf2 - r^D 


— + — — (r^v Y -r^D 
dt^r^drV^ 



_ ‘hi 

dr ) 

p 

9Y,2\ 

- hi 

dr J 

p 

dYx\ 

- h. 

dr j 

p 


(2.37a) 

(2.376) 

(2.37c) 


It is assumed here that all the species produced during combustion are represented by one 
single product, p such that 


t=l 


Y, = l-Yf,-Yf2-Y, 


(2.38) 


Once Yfi,Yf 2 and Yx are solved by the three species Eqs. 2.37a,b,c ,the average mass- 
fraction for the product can be obtained from the above relation. 

(ii) Energy conservation equation 
The gas phase energy Eq. 2.30 is 

I ^ t +2- [« ( /J CvdTj - pDs. ( /J + (pDC, - A) vt] = - E '>?“< 


Using spherical symmetry, 
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2.4.3 Boundary Conditions 

Two factors control or influence any physical system. They include the processes within 
the system, as well as the interaction of the system with it’s surrounding. So, unless the 
boundary conditions are properly specified or implemented, the analysis of any physical 
system remains incomplete. The physical system chosen for the present problem includes 
both liquid and gas phases. The balance of species, mass and energy fluxes at the interface 
of the liquid and the gas phases lead to the coupling conditions between the two phases. 
The coupling parameters are, obviously, functions of time. The boundary condition at 
the droplet centre is obtained from the condition of symmetry. The condition of symmetry 
signifies that there is no transfer of species or energy across the droplet centre. The 
boundary conditions at the outer boundary of the gas phase solution domain come from 
free stream conditions. 

(a) Droplet centre (r = 0) : Using the condition of symmetry at the droplet centre, 
boundary conditions can be written as 


f ' = « 

(2.40a) 

to 

II 

o 

(2.406) 

o 

II 

(2.40c) 


(b) Interface of liquid and gas phases i.e. Droplet surface (r = 1) 

{i)rnass balance: 

In absence of any flow in 9 or, (j) direction, all the tangential velocity and shear stress 
components at the liquid- gas interface are zero. The radial flow, which is generated in 
the gas phase due to the surface regression of the liquid droplet, is the only component of 
velocity that prevails in the gas phase. Therefore, if we neglect the value of p in comparison 
to p because p » p, the mass -balance at the interface, at a given time gives rise to the 
condition 

J-f-d 

ptid = - pUd = ~ P~^ 


(2.41) 



where vj, is the droplet radius at the given instant, (^) is the droplet regression rate, 
uj, the evaporation velocity of the evaporating fuel vapour at the droplet surface at that 
given instant. Through out the thesis, the positive value of uj, is always used. Therefore 
whenever uj, is used in equations, this information is kept in mind. 

{ii)species balance 

In both the phases, species are trairsported due to both convection as well as diffusion. 
Therefore, species flux, n,- for the species in the gas phase is given as 

ni = piv + l. ( 2 . 42 ) 


where, piV is the convective component and ji is the diffusive component of the species 
flux, V is the bulk velocity and pi is the density for the species. For the present problem, 
since there is only radial component, the above expression becomes 


= PiVr + it,r 

where, j,-,,. = piK-.r = (using Pick’s law) 

or. 


n,> = pYiVr - pD 


Similarly, for the liquid phase. 


hi^r = pYiVr - pD 


dr 


dYi 

dr 



( 2 . 43 ) 


( 2 . 44 ) 


Elamtntal volumi in sphtrical coordinatts 

Consider a differential volume ( as shown above ) . r'^dO sinO dcj) dr in the spherical 
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polar coordinates ( see Fig. 2.2 ). The net efflux of the species i in r-direction is 

= ["‘.r Ir+Ar “ Ir-Af] ded<t> 

= {^YiVr — — ^pYiVr — pD-^^ r^sinO d9d<j) 

Now, efflux + storage = production. In the limiting case, as Ar — > OandAr — 0, the 
control volume approaches a control surface. Assuming that the storage and production 
of species are zero at the interface, the expression for the net efflux of species i can be 
expressed as 


pYiVr - pD 


dr 


r=Td 




r=:rdj 


r^sinO dOdcj) = 0 


or. 


/ BY' 

^’^sinO ded(j>p \YiVr -D^ 


T=Ti 


( -dYi 


r^sind dddift 


(2.45) 


r=rd 


The L.H.S. of the above expression is the species flux, m, for the component at the gas 
side interface and R.H.S. is the species flux, rh,- for the component at the liquid side 
interface. Since, there is no production or storage at the interface. 


rrii = m,- 


(2.46) 


Thus, by physical analysis also, one can arrive at the same expression ( see Eq. 2.45 ) for 
species conservation. Rewriting the species flux, rhi for the gas phase, 

DdY 


r^sind d6d<f>pVr 

Since, Vr = wj and if = K, at r = rj, 

r^sind dO d4>pud 


Yi- 


Vr dr 


= mf 


j rz=:rd 


Yis- 


DdYi 


Ud dr 


mi 


J rt=rd 


Again, m = r^sinS dO d4> pui = total mass evaporation rate. So, 

„ D dYi . 


m 


Ud dr m 


7 ^ = Cfi[say) 


(2.47) 


(2.48) 
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Final expression for gas side species conservation equation can be written as 

0 


dYi Ui 


(2.49) 


where, i stands for fuell (/I), fuel2 (/2) and oxidiser (x) of the gas phase and e,- is called 
the fractional mass-gasification rate for the species. So, gas side species boundary 
conditions for different species can be written as 

dYnl 


dr 

dYf, 


r=rd 


+ ^ (^/i “ Yfu) = 0 


dr 


r=ra 


+ ^ (e/2 - Yf2s) - 0 


aK 


dr 


r=rj 




(2.50a) 

(2.506) 

(2.50c) 


Similarly, the liquid side species boundary conditions for different species can be 
written as 


dY 




dr 

dYf2 


+ 1 (e,. - y,.,) = 0 


r=rd 


(2.51a) 


dr 


r=rd 


+ ^ (c /2 - = 0 (2.516) 

Here tfi and e/j represents the fractional mass gasification rates for fuel 1 and fuel 2 
respectively. 

(iv)heat balance 

In evaluating the overall heat balance at the liquid-gas interface, we find that in absence 
of any storage or production of energy at the interface, the total amount of heat available 
from the gas phase is used up for droplet heating as well as for fuel evaporation at the 
droplet surface. Therefore, the heat balance at the interface can be written as 


^ dr 


r=rj 


^ dr 


r=rd 


N 

+ P'^dYl 

t=l 


(2.52) 


where, p is the mixture density for the gas phase, Ud is the evaporation velocity, e,- and X,- 
are respectively the fractional mass-gasification rate and the latent heat of vaporization for f 
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the species, A and A are the thermal conductivities for the gas phase and liquid phase 
respectively. 

{v)phase equilibrium at interface 

For the present problem, the component in the liquid phase are assumed to form a miscible 
mixture. For a miscible mixture, the vapour pressure (P,,) of species i depends on it’s liquid 
phase mole fraction (Xis) and on it’s vapour pressure, when it is pure [ Law, 1982 ]. It 
is assumed that the molecular structure for the liquid phase constituents are sufficiently 
similar. In such a case, the mixture can be assumed to be ideal so that the Raoult’s Law 
remains valid. Therefore, 

Pis = 


or. 


■■x,sXl 


(2.53) 


where, Xia is the mole fraction for the i^^ species on the gas phase side and X,® is the 
concentration for the species on the gas phase side when that species is pure. Now, X? 
is given by Clausius-Clapeyron relation as 




(2.54) 

where, T, is the surface temperature, Tbi is the normal boiling point evaluated at the 
normal atmospheric pressure, P„, P is the prevailing system pressure and Wi and Li are 
the molecular weight and latent heat of vaporization respectively for the z*'" species, and 
Ru is the universal gas constant. Three fundamental assumptions are involved in Eq. 2.54. 
Firstly, the gas near the liquid surface is ideal. Secondly, the vapour pressure of the fuel is 
not affected by the presence of any other species. Finally, finite mass-transfer rate at the 
droplet surface has negligible influence on the equilibrium condition. Now, by using the 
Raoult’s Law and the Clausius-Clapeyron relation, the gas side surface mole fractions for 
the two fuels can be expressed as 

= (2-55) 



1 

--'ll 

Ru 

U/r 

T.)\ 
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Xf2, = A'/2«exp 


Wj2Li2 ( 1 


Ts. 


(2.55b) 


Ru \Tbf2 

(c) Outer boundary of the gas phase domain(r oo) 

This outer boundary of the gas phase domain is far away from the reaction zone. Therefore, 
there is no influnce of combustion at this boundary. So, it can be safely assumed that the 
free stream conditions prevail at the outer boundary. Therefore, the boundary conditions 
at r oo can be written as 


T = T^ 

Yfi = 0, 

Yj2 = 0, 

Y. = r.,00 

where the subscript, oo refers to the free stream conditions. 


(2.56a) 

(2.56t) 

(2.56c) 

(2.56d) 


2.4.4 Initial Conditions 


When micron-size droplets axe sprayed into the combustion chamber, combustion can staxt 
in two ways. When the initial temperature of the combustion chamber is high (e.g. 1000 
K) and uniform, combustion will staxt near the droplet surface where the stoichiometric 
ratio of the fuel and the oxidiser is prevalent. Combustion may also start at a fixed point by 
supplying high temperature spaxk (>1500 K) at that point in an otherwise normal ambient 
(300 K) atmosphere. Initially, the gas phase is completely devoid of any of the fuels. Also, | 
the initial mass-fraction profiles, as well as the temperature profile for the liquid phase are I 
uniform. Therefore, the initial conditions for both the liquid and the gas phases can be | 
expressed as | 

(a) Liquid phase 


t (r,0) = fin, 

(2.57a) 

Yn {r,0) = Yin.n, 

(2.566) 

Y„(r,n) = Yin,,2 

(2.57c) 



(b) Gas phase 


T{r,0) = T,„, 

(2.58a) 

Yn{r,0) = Yi„,n, 

(2.586) 


(2.58c) 

K (r, 0) = Kin,, 

(2.58d) 


Present study uses both of the above mentioned methods for starting ignition. 

2.4.5 Non-dimensionalization 


Non-dimensionalization of the governing differential equations clearly reveals the relative 
importance of the various aspects of the transport processes. This gives rise to some very 
importapit non-dimensional quantities, such as, Schmidt number, Prandtl number, Lewis 
number etc. Schmidt number indicates the relative importance of species diffusion over 
viscous diffusion, Prandtl number indicates the relative importance of heat diffusion over 
viscous diffusion and Lewis number indicates the effect of heat transfer over mass transfer. 
Non-dimensional quantities related to the present problem are indicated by the asterisk 
(♦) mark as 
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(2.59) 


For the sake of convenience, from now onwards, the asterisk (*) mark will be omitted if 


otherwise not mentioned. 

Using the above parameters, the governing differential eqs. 2.32, 2.34, 2.37 and 2.39, 
the boundary conditions given by Eqs. 2.40, 2.41, 2.50, 2.51, 2.52, 2.55 and 2.56, and 



initial conditions given by Eqs. 2.57 and 2.58 can be nondimensionalized and the following 
equations are obtained. 

A. Non-Dimensional Governing Differential Equations 


(a) Liquid phase 


dt p dr \ dr ) 

dt p dr r^ dr I dr J 


dT p ^ dT Lea d 
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where Yr = 

B. Non-Dimensional Boundary Conditions 

(a) Droplet centre (r = 0) 


dYn 

dr 

dYn 

dr 


?L 

dr 


0 , 

0 . 

= 0 


(2.60a) 

(2.606) 

(2.60c) 


(2.61a) 

(2.616) 

(2.61c) 


(2.61d) 


(2.62a) 
(2.626) 
(2.62 c) 


.dr^ 

pui - -pUi ^ 


(2.63) 


(b) Droplet surface (r = 1) 
(i)mass balance 



(ii)gas side species balance 


BY,, 
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BY, 2 


dr 
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dYJ 


dr 


+ (e/i - Yfu) = 0 

+ (^/2 - Yj 2 s ) = 0 

+ ^ (-K.) = 0 


r=l 


(iii)liquid side species balance 

dY,i 


dr 

BY, 2 


+ (^/I “ ^/I^) = 0 


r=l 


{iv)heat balance 
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{v)phase equilibrium at interface 

Xfu =Xfisexp 


Xf2s = ^/2aexp 
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where, 
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Wf^L 
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and, /^/2 = 


jRu 


(c) At outer boundary (r -+ oo) 


T = Too, 

1/1 = 0 

1/2 = 0 

F, = F.,oo 

C. Non-Dimensional Initial Conditions 


(2.64a) 

(2.646) 

(2.64c) 

(2.65a) 

(2.656) 

( 2 . 66 ) 

(2.67a) 

(2.676) 

(2.67c) 

(2.68a) 

( 2 . 686 ) 

(2.68 c) 

(2.68d) 
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(a) Liquid phase 


(b) Gas phase 


r(r,o) = r,„. 

(2.69a) 


(2.696) 

Vy, (r,0)=r,n,f2 

(2.69c) 

T(r,0} = Tin 

(2.70a) 

Vfi (r,0) = y;-„.yi 

(2.706) 

Vf 2 (r,0)=V,nJ2 

(2.70c) 

K(r,0) = Vi„,,. 

(2.70d) 


2.5 Thermo- chemistry 

Chemistry involved in any combustion process is practically very very complicated. Though, 
the reactions here are represented by a single step global reaction, the factual chemical re- 
actions axe almost infinite in number before the final products are generated. Therefore, 
whatever simplified model one may select, the knowledge of heat of reaction and rate 
of chemical reaction are essential. 

2.5.1 Heat of Reaction 

The science of thermo-chemistry is concerned with the energy change associated with; 
chemical reaction or in other words, it deals with the conversion of chemical energy into) 
heat energy and vice-versa. A chemical reaction is merely a recombination of elements; 
to form new products. This reaction may be endothermic or exothermic, depending on: 
whether heat is being absorbed or released in the reaction. Obviously, a combustion process! 
is necessarily exothermic. There are a number of ways in which heat of reaction can be 
defined or measured. In general, the heat of reaction is defined as the energy released 
by the system when the products are restored to the same temperature as that of the 
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reactants. This amount depends on reactant temperature and the variation of the system 
pressure. It is customary to measure the heat realesed at standard atmospheric pressure. 
For a general chemical reaction involving N species can be written as 

( 2 . 71 ) 

t=:l tz=l 

The heat of reaction equals the change in enthalpy of the system under constant pressure. 
It is calculated from heat of formation of the reactants and the products as 

Ah = '£v;'h°-f:y;h1 ( 2 . 72 ) 

t = l 1=1 

at standard state at Tq where h° = heat of formation of species, i, at temperature, 
v\ and v'l = stoichiometric coefficients of reactants and products respectively. 

Mi = arbitrary specification of of chemical species. The heat of formation, hi of a species 
i at any temperature T , is defined as the enthalpy of the substance relative to that of 
its constituents elements at the same temperature. It is to be noted here that the heat of 
reaction is used to find out the ADIABATIC FLAME TEMPERATURE. If a combustion 
process takes place adiabatically such that there is no work or change in kinetic or potential 
energy then for such a process, the temperature of the products is referred to as the 
adiabatic flame temperature. This is the maximum temperature that can be achieved 
ideally for given reactants. This is because any heat transfer from the reacting substances 
and any incomplete combustion would tend to lower the temperature of the products. Thisi 
adabatic temperature can be controlled with the amount of excess air used for combustion.: 

2.5.2 Chemical Kinetics 

It is that part of chemical science that deals with the quantitative study of the rated 
of chemical reactions and the factors upon which they depend. It also deals with the 
interpretation of the empirical kinetic laws in terms of reaction mechanisms. All chemical 
reactions take place at a definite rate, depending on concentration, temperature, pressurf 
catalyst etc. For a general combustion reaction involving N species such as given by Eq 
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2.71, the individual rate of reaction (wi) may be expressed as the rate of decrease of 
concentration (or, consumption) of a reactant i in moles/m^/sec as 

0,. = (2.73) 

where, u; is the overall reaction rate for the above reaction. From conservation of mass, 

N 

= 0 (2-74) 

t=l 

The law of mass action states that the rate of disappearance of a chemical species is 
proportional to the products of the concentrations of the reacting chemical species, each 
concentration being raised to a power equal to the corresponding stiochiometric coefficients. 
Thus the reaction rate for a combustion reaction is given as (Kuo, 1986) 

M 

^==K]l{CMir' (2.75) 

t=i 

where, M is the number of reacting species and K is the proportionality constant known 
35 the specific reaction rate constant with a strong temperature dependence. The above 
relation is valid only for simple reactions. But for complex reactions represented by a single 
step, the above relation can also be used to serve the practical purpose. In such a case v[ 
is known as the order of reaction. This order of reaction is understood to be the number of 
colliding particles which participate at the molecular level. For a simple gas phase reaction 
taking place at a normal pressure, the order of reaction is not expected to be more than 2. i 
The specific rate constant, K, can be estimated from Arrhenius equation as (Kuo, 1986) ; 

if = XT“exp(-^) (2.76) 

where, E is known as the activation energy, ot and A are constants, a lies between 0 and I,; 
E/R indicates the temperature below which K is small. This law reveals the fact that only] 
•those molecules which have sufficient energy to overcome the opposing activation energy 
(E) and approach each other from preferred directions can partcipate in the combustion 
reaction. Here we see that this generation term (w) is highly non-linear which makes the 

OTRAL library : I 

proper analysis of combustion process quite difficult. I I. T., KANPUR 

*«**^aA 123594 

■ .... .*f ;■<- . 
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2.5.3 Dissociation Reactions 


At high temperature, dissociation of the products becomes very important, as it brings 
down the flame temperature in practical situations. The dissociation reactions are repre- 
sented as 

C02^C0 + (2.77) 

+ (2.78) 

1 

Kpi = (2.79) 

■^C02 

1 

P p2 

Kp2 = (2.80) 

Where, Kpi and A'p 2 are the two equilibrium constants for the above two reactions which 
are functions of temperature(Kuo,1986). Considering the reaction for fuell, ^coji, 

are the stoichiometric coefficients. They are dependent on the equilib- 
rium constants and also on the type of fuel being burnt. 

r'r/i = + vcoji + ^o^fi + (2.81) 


Relations for Kpi and Kp 2 can be simplified to yield 


K. 
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r'co/i 
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VHiOJl, 

VTfl . 
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11/2 


7l//i = Vcofl + l^COa/l 


(2.82) 

(2.83[ 

(2.84i 


2vxfi — 2vo2f\ + I'coji + ^cOifi + ^HiOji (2.85; 

16vfi = 2uH20n + 2uH2fi ( 2.80 

Solving Eqs. 2.81 through 2.86 by Newton-Raphson iterative scheme, the stoichiometrj 
coeflScients corresponding to different temperatures, can be determined. 
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2.6 Closure 

This chapter deals with the modeling of spray combustion. The governing differential 
equations, the boundary conditions and the initial conditions in the nondimensional form 
for liquid and gas phases are given by Eqs. 2.60 through 2.70 which need to be solved for 
determining the mass fractions of the fuels fl and f2 in the liquid and gas phases and that 
of oxidant in the gas phase. The temperature time history can also be determined. 



Chapter 3 

Semi Analytical Solution 



3.1 Introduction 


[n practical situations, atomised fuel droplets are generally very small in size. In such a 
case, when they are sprayed into a combustion chamber, they move with velocities almost 
equal to that of the gas. Keeping such situations in mind, present study stresses upon 
the study of transient burning for a multicomponent single droplet fuel in an otherwise 
quiescent atmosphere. The initial radius of the droplet is of the order of 100 microns. 
The fuel is atomised and injected into the combustion chamber in a fairly hot (1000 K) 
atmosphere. There is no generation of species inside the droplet. Therefore, the net mass 
transfer inside the droplet is zero. The spray is assumed to be dilute in nature, so that 
the burning of one single droplet is not influenced by the presence of the neighbouring 
droplets. Stress is also put on the multicomponent aspect of droplet burning. Though 
the neglect of the external flow as well as the influence of neighbouring droplets is not 
physically realistic, it still gives enough physical insight into the fundamental aspects of 
multicomponent fuel burning. These two above mentioned major assumptions do not affect 
the fundamental aspects of fuel burning in general. A schematic diagram of such a droplet 
burning in a quiescent atmosphere is shown in fig. 2.1. Since, there is no external flow 
around the droplet and hence, the evaporation of fuel vapour at the droplet surface is 
only radial, the combustion process becomes symmetric in 9 and 4 > directions. So, the 
process becomes a function of radial distance (r) and time (t) only. A coordinate system 
(r), fixed to the droplet centre, is considered here. In absence of any external flow, the 
droplet surface as well as the surrounding flame remain spherical. There is an inherent 
tendency of the droplet to remain spherical when the size of the droplet is very small. 
Due to continuous evaporation, the droplet surface regresses continually giving rise to a 
moving boundary. When the stationary coordinate is transformed into moving coordinate 
such that the droplet surface always remains at r=l, this moving boundary gives rise to a 
convective term in the liquid phase transport equations. For the present study, transient 
feature of liquid phase species and energy equations are retained, while the gas phase is 
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.Tied to be quasi-steady because, any change in the gas phase occurs in a much smaller 
time-scale compared to that for the liquid phase. Properties for both the phases are taken 
to be constants. The combustion takes place at normal atmospheric pressure and the 
liquid fuel has two species. The phenomenon of micro-explosion is not considered here. 
The assumptions of thin flame and unit Lewis number for the gas phase are retained. 
Dissociation of the products at high temperature near the flame is neglected. Schvab- 
Zeldovich variables ( to be described later ) were incorporated to eliminate the non-linear 
generation term (uj) in the gas phase. The quasi-steady gas phase is solved analytically. 
The free stream conditions are applied at the outer boundary of the gas phase domain. 
The liquid phase equations are solved in a semi-analytical manner. The general solution 
for the liquid phase equations is constructed in a series form for each discrete time step At, 
assuming that the gradients at the droplet surface remains constant over this time period 
At. Then the overall solution is obtained in a quasi-linearized form using superposition 
principle. 

3.2 Mathematical Formulation 


For the liquid phase 
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Now, using y/i = and equations (3.1a-b) can be rewritten as 
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Since p » p, the convective part arising out of droplet surface regression (^) is small for 
cosiderable part of the droplet burning. And also, for the sake of simplicity, this convective 
part of liquid phase equations are neglected. The gas phase Lewis number is assumed to 
be unity and, using = 3 ^, and Therefore, eqs. 3.2a, 3.2b and 3.2c can be 

written as: 
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(3.3a) 


(3.36) 


(3.3c) 


For the quasi-steady gas phase ; r > 1 

The gas phase Eqs. 2.61a through 2.61d for fuell, fuel2, oxidiser, and temperature can also 
be simplified in a similar manner as that used for simplifying the liquid phase equations 
and by using Ud = (from mass conservation), they can be written as 
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PooDoO p 


Ah 


where^ 
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‘^/i? ^/ 2 ? — rate of consumptions of species fuell, fuel2 and oxidiser respectively in 

moles /m^/ sec. 

Now, using 


^/1 5 ^ f2 


Yf 


y. 


and from Eq. 2.73, 


u 


d’/I _ 0Jf2 (^x 


(3.5a) 


(3.56) 


^/1 ^f2 

where, Ufi, z// 2 , and are the stoichiometric coefficients for fuell, fuel2, and oxidiser, the 
species equations for fuell, fuel2 and the oxidiser can be written as 


dY, 


dY;, 1 d 


and 


+ T^c r-^ + 

at p or or 

dYf, , 3Y;, , 1 d 

+ ’•-aT + ^a;: 


dt '‘ dr r2 ar 


Sc UdYJj — Dr 




Sc UdYj2 — Dr 


dr 


dr 

dY* 

ScudY:-Dr^-^ 


== LO 


UJ 


= to 


(3.6a) 

(3.66) 

(3.6c) 


For the gas phase energy equation. 


N 


y~) = h^iUjfi + hPj2'^f2 + h^x^x + ^pd;p 


»=i 


and since from Eq. 2.73, 


u = 


^/i a>/2 _ ujx _ ^ (^p 


+■ 


Vfl Vf2 ^p 


Therefore, 


Where, 


N 

— h’^uji = — ta [^'p6p — ~ ^'/ 26 y 2 ” “ d;A6 

1=1 


Ah = —Uphp + Ufik'fi+ Vf2h°f2 + ^p^x 

and i/p, i//i, r// 2 ,i/x are the stoichiometric coefficients for products, p and the reactants, fl, 

/ 

f2 and x. Therefore, the above Eqs. 3.6a-d becomes 

mM 

ch* (3.6d) 


dt p dr r^dr 


Sc Ud (CpT) — a Le r^ 
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where, ch* = ^ ^ = overall reaction rate for the above reaction in moles/m^/sec. 
For reasons given earlier, the convective part containing Udj is neglected in Eqs. 3.6a-d, 
and incorporating the gas phase quasi-steadiness, and using a — D because gas phase 
Lewis number is unity, the species and energy equation for the gas phase become 


dr 


i-A 

dr 

LA. 

dr 


f jy* 

Sc — r^D — ^ 


Sc UdYj2 — 


LL 

r^ dr 


Sc Ud {CpT) - r^D 


dr y 

r 


1 - 

dr ) 

1 

dY:\ 


dr J 


d{C,T)\ 


u 


u* 


dr 


u 


(3.7a) 


(3.7b) 


(3.7c) 


(3.7d) 


The non-dimensional boundary and initial conditions given by Eqs. 2.62 through 2.70 can 
be rewritten in terms of variables defined by Eq. 3.5 as: 
a. At the droplet centre (r = 0), 


II 

o 

(3.8a) 

dYT^ ^ 

-^ = 0 
or 

(3.86) 

^ = 0 

(3.8 c) 



60 


b. At droplet surface (7- = 1) , 


(i)Mass conservation 


puj, = -pui = -p— 
at 

where, Ud is the evaporatioir velocity and Ud is the radius regression rate. 


(3.9) 


(ii) Gas-side flux balance for each of the species ( two fuels and an oxidiser ), 

+ ^ (n,i£/. - YJ„) = 0 (3.10a) 

+ - y/as) = 0 (3.106) 

+ ^(-y;.) = o (3.10c) 

r=l 

where, O/i = — 0/2 = ^/i ^/2 the fractional mass gasification rate for 

fuell and fuel2 respectively. 


BY 


n 


dv 


r=l 


ay;, 


dr 


r=l 


(iii)Liquid-side flux balance for the fuels, 


dYj\ 

dr 

aa Sc 
^ D * 

{piicji - y/ij 

d}% 

dr 

T* — ! 

aa Sc 
+ D ' 

r=l 

- y;,.) 


where, aa = pud/p 


0 

0 


(3.11a) 

(3.116) 
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(iv) From heat balance at the interface, 

. Qf 

= pud{tjiLfi tfiL ^ 2 ) (3-12) 

r=:i or 

where, Lji and Lf 2 are the latent of evaporation for fuel 1 and 2 respectively. 

(v) phase equilibrium at interface : 

Assuming that the liquid fuel behaves as an ideal mi.xture and the fuel vapour is saturated 
corresponding to the surface temperature, the Raoult's law and the Clasius-Clapeyron 
relation can be considered valid. Accordingly, the liquid -phase and the gas-phase mole- 
fractions are related through the expressions; 
at r = 1.0 , 



Xfu = Xn,exp{/3fr{l/%n-l/Z)) (3.13a) 

Xf2s = Xf2sexp(pf2{llT,j2 - IfZ)) (3.136) 



T — T 

j- — -^00 

(3.14a) 

II 

0 

0 

(3.146) 

11 

0 

0 

(3.14c) 

y = Y* 

(3.14d) 



62 


(d) The initial conditions are: 


r(r,o) = rj„ 

(3.15a) 


(3.156) 


(3.15c) 


3.3 Analytical Solutions 

3.3.1 Gas-Phase 

It is seen that due to the presence of the generation term (w*), the gas phase Eqs. 3.7 
are highly non-linear. Analytical solutions are not possible as long as Cj* is present in 
the equations. This generation term is eliminated by using thin flame approximation 
and consequently, applying the Schvab-Zeldovich variables (Kuo, 1986). In thin flame 
approximation, the fuels and the oxidiser are assumed to be fully consumed at the flame 
surface so that the oxidiser mass-fraction is always zero inside the flame and the fuel mass-: 
fractions are always zero outside the flame. The Schvab-Zeldovich variables are defined 
as ; 


(jri ~ 

' (3.16a) 

g^ = y;,-y: 

(3.166; 

h = c,t-y: 

(3.16d 

Now, subtracting eq. 3.7c from eqs. 3.7a, 3.7b and 3.7d and then applying the Schvab 

Zeldovich variables, we have 


(3.176 

--(scu,G, 

(3.1 7i 
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Sc UdH - r^D- 


(3.17c) 


By virtue of thin flame approximation, for 1 < r < rj, Y* = 0. Therefore, in this 
range, 

G’l = y/i (3.18a) 


= y;. 


H = CpT 


(3.18a) 

(3.185) 

(3.18c) 


Again, for rj<r< co, y)*j = 0 andV^j == 0- Therefore, in this range, 

= G2 = -y; 

and 

h = c,t-y: 


(3.19a) 


H = C^T- y; ( 3 . 195 ) 

The boundary conditions for the ga.s phase in terms of Schvab-Zeldovich variables can be 
written as 
at r=l 


dG2 UdSc 


At r — + oo. 


- [fl/iCyi - 

(3.20a) 


(3.205) 


(3.21) 

— V 

"2 07,00 

(3.22 a) 

Hx) 5^,00 

(3.225) 


Equations 3.17a through 3.17c along with the boundary conditions represented by Eqs. 
3.20a through 3.22b can be solved analytically in terms of the Schvab-Zeldovich variables. 
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Gi, G 2 and H. Then, by using the relations 3.16a,b,c, the final analytical expressions for 
species and temperature of the gas phase can be given as 


for 1 < r < r/, 


y;, = 


Yj 2 — 0,f 2 


'Y 


tfi - 

C/2 - 




n 


/I 




n 


J2 


+ ^J2 


N / UdSc\ 


nn ^ m V* 

X — Joo Jx,oo 




c/iT/i + C/2T/2 + 




dr 

r=l 


UdSc 


exp 


UdSc 

~D~ 


exp 


/ UdSc\ 

\'W7) 


for Too > r > r/, 


y: = -n 


n 


(Yx,oo , ^ f ‘^dSc 

Cfi - { — h Cfi j exp 


Cl 


n 


D r 


T = 

+ ^flLfi + €/ 2-^'/2 + 




r=l 


UdSc 


(3.23a) 

(3.236) 


(3.23c) 


(3.24a) 




3.3.2 Liquid-Phase 


Equations 3.3a, b,c along with boundary and initial conditions given by Eqs. 3.8, 3.11, 3.12 
and 3.15, can be transformed defining new variables of the form : 

Ui = rY; and ?h = rf ' (3.25) 

Where, i corresponds to the two components of the liquid phase. The resulting governing 


equations are : 

dui _ 

(3.26a 


dt dr^ 


dw ^d'^w 

~di 

(3.26i 


The transformed boundary conditions from Eqs. 3.11 and 3.12 are given by 
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at r = 1 : 


where, 


and 


and at r = 0 : 


Ci{t) 


dui 
dr 
dw 

^-Crd) 


C< it) = [l + «... - 


Cr(t) = l(x§ 


- +'^S 

r=l 1 


(3.27a) 

(3.276) 


(3.28a) 


(3.286) 


u = 0;u; = 0. (3.29) 

The initial conditions from Eq. 3.15 are: 

Ui (r, 0) = rT’in . (r) (3.30a) 

tt; (r,0) = (r) (3.306) i 

The general solutions for the above system can be constructed in a series form for each 
discrete time step At, using the solutions for the constant gradient case . The fundamental 
solutions, corresponding to constant interfacial gradients, can be obtained in the following i 
way. In solving eq. 3.26a, let us assume that the solution, consists of two parts namely | 
steady state and transient parts as given by 

Uj = «(r) + v(r, <) (3.31) ; 

Then, = 0 and vt = Dvrr 

Let us assume the steady state boundary conditions as I 

at r=l, ; 

= Ci (3.32a) j 

r=l i 


du 

dr 
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and, at r=0, 


u (r) = 0 


(3.326) 


Therefore, applying the above steady state boundary conditions, we get the steady state 
solution as 

u (r) = Cir (3.33) 

The corresponding transient boundary conditions for Vt = Dvrr yields 


v,{r = l,t) = Ci-Ci = Q (3.34a) 

u (r = 0, i) = 0 (3.346) 

and the initial condition becomes 

u(r,i = 0) = -r (C, -Kn) (3.35) 


Analytical solution for the above transient Eq. 3.26 can be obtained by the method of 
separation of variables by using the above boundary conditions and can be written in a 
series form as ( Myers, 1971 ) 

OO 

v{r,t) = ^ A„sin(p7r r) exp (p^ir^Dtj (3.36) 

71 = 0 ) 


Now, satisfying the initial condition, we get 

— 2(7 /— li" /■! - 

+ 2 / (r) r sin (p7rr) dr (3.37) 

Jo ■ 

where, p — i = /I (for fuell) and /2 (for fuel2). Therefore, the complete solutions 
for Ui (r, t) can be expressed as 


tli(r,<) = C.r-I- AnSinpTT r exp(pV^i9<) (3.38) 

n=0) 

since, u,- = rY*, 

1 OO 

Y*{r,t) = Ci + -^A„iSm{pirr)exp(^-p^'ir^bty (3.39a) 

n=0 
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where, 

+ 2 ^ i^in,. ^ sin (p7rr) c/r (3.396) 

and, p = , i = /1( for fuel 1 ) and /2(for fuel 2 ), The analytical solution to Eq. 3.26b 

for the liquid side temperature can also be obtained in a similar manner and is expressed 
as 

1 OO 

T {r,t) = Ct + - Y2 sin (pxr) exp {—p^'K^a tj (3.40a) 

^ n=0 

where, 

AnT = (“!)" + ^ ^in ^ sin (pirr) dr (3.406) 

Assuming that the surface gradients remain constant over the period t to {t + At), the 
quasi-linearized solutions at t = mAt are obtained using the superposition principle as 
follows : 


At t=0, ¥;•> = y;-;,. = C? and Aj, = 0 

Now, at t = At 


(r, t = At) = C- >intSin (pirr) exp At) 

n=0 

= y;^ + - y;^) 

1 OO 

= Y;° + (cl - cf) -f - 51 (^L- - ^m) sin (pTrr) exp (-xx) 

^ n=0 

where, 


^ni 


where, xx = p^ir^DAt. 
Similarly, at t = 2At 


2(-l)" 

2(-l)" 


(Vi- c’) 

y; (C? - C/) exp {-k xx) 


i;.*2(r,t = 2At) = 


1 °° 

c? -f- - 5D ^m-sin (pTTr) exp (-xa;) 

n=0 




68 


= + {cf - C}) + - £ [Ali - Al) sin (pTrr) exp (-xx) 

= Y*° + {cl - Cf) + - £ (a;„ - Ali) sin (pxr) exp (-xx) 

^ n=0 

1 oo 

+ (cf - cf ) + -Y1 {^i - ^ni) sin (p^^) 

^ ti=n 


where, 


4^. 


^^2 + (cf — Cf) exp (— xx)| 

y - cf-'') exp i-k xx) 

pV * 


And, at < = 3A< 


1 

y;*3 (r, t = 3 At) = Cf + - Af.sin (pxr) exp (-xx) 

= Y;^ + {Yr^-yr') 

= Y-^ + (cf - Cf) + i E - ^i) ““ 

^ ^ n=0 

= f;° + (cf - cf) + - E (^n.- - ^n.) sin (pTrr) exp (-xx) 
+ (cf - Cf) + - £ (aL- - ^ni) sin (pxr) exp (-xx) 

+ (cf - Cf) + - £ (Af, - A^.) sin (pxr) exp (-xx) 

^ f n=0 


where, 


4^. 

/lyj* 


: [ {cf- Cf) exp (-O.xx) + (cf - Cf) exp (-xx) + (Cf - Cf ) exp ( 2xx)] 

p2^2 L\ 

_ £ {Cf~^~^ - cf'’') exp (-fc xx) 

Thus, at t = mAt, the solution can be written in series form as 


y; =Y* ° + E^^' 

j=l 


* 3 


t*o+if:({oi - cr^)r+Xi<^ - 

» * r "=® 
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where, xx = p^z^DAt 

E (<?/-*-' - C/-‘) exp (-k m) (3.416) 

* fc=0 

Similarly, solution for the liquid side temperature can be expressed as 

A ^ ^ 

= T^ + 

i=i 

r"* = r° + - £ f r+ (Kt - Kt^) exp {-yy) sin (pxr) ) 

i=l X n=:0 / 

where, yy = p'^z^ctAi 

and, Ai^ = J2 - Cf’'') exp {-k yy) 

3.4 Numerical Evaluation 

The liquid phase temperature and mass fraction profiles were evaluated using Eqs. 3.41 and 
3.42 by truncating the infinite series to 100 terms. The semi-analytical solution was also 
compared with the finite difference solution of the same governing differential equations 
and excellent agreement (up to 4 places of decimals) was observed between them. It was 
found that the semi-analytical solution was about ten times faster than the numerical 
solution. For the finite difference procedure it was necessary to use at least a mesh size 
of 1/1000 within the droplet, for achieving a good match with the analytical results . 
Another advantageous feature of the semi- analytical procedure is that it could be applied 
to calculate only the surface quantities for each time step, if so desired. The process time 

required for evaluating the solution for 100 time steps by these three methods are given as: 

i) 250 seconds for numerical approach. 

ii) 30 seconds for full semi-analytical approach. 

iii) 4 seconds for evaluating only the surface parameters through semi-^atfalytical pro- 


(3.42a) 


(3.426) 


cedure. 
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3.5 RESULTS AND DISCUSSION 

Results have been predicted for the liquid phase temperature and concentration profiles and 
the fractional as well as the total gasification rates. The fuel combinations of n-heptane, 
n-octane, n-decane, undecane and dodecane have been considered, over the initial mass 
fraction range ( 0.1-0.99 ) of the more volatile component. In order to highlight the salient 
features of multicomponent fuel burning, results are presented for a binary mixture of 
heptane and decane with an initial proportion of ( 0.99 hep.- 0.01 dod ), in the following 
sections. 

3.5.1 Transient Features of Multicomponent Combustion 

The time history of mass-fraction and temperature profiles of a heptane- dodecane droplet 
is depicted in Figs. 3.1(a,i). Shortly after the introduction of the droplet, a concentration 
boundary layer is established near the surface due to rapid evaporation of the more volatile 
component. The concentration gradient near the surface is maintained almost throughout 
the droplet lifetime by the high liquid side diffusional resistance. This trend also implies a 
high surface regression rate compared to the liquid side mass diffusion rate. The boundary 
layer penetrates more into the core as time progresses. 

In contrast to the mass-fraction profiles, the temperature profiles do not show the 
existence of a steep thermal boundary layer. This can be attributed to the relatively 
higher liquid side thermal diffusivity, because of which heat is conducted very fast to the 
core of the droplet. Fig.3.1fe also indicates the presence of two distinct d^-law periods. In 
the first <P-law period, the droplet surface temperature remains nearly constant around 
the boiling point of the more volatile component. At the end of the first cP-law period 
(when the more volatile component is nearly depleted), rapid droplet heating occurs once 
again to raise it’s temperature closer to the boiling point of the less volatile fuel. This is 
followed by the second <P-law period for the less volatile fuel. If the degree of superheat 
for the more volatile component is high, microexplosion is likely to occur, which in the 
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Fig. 3.1 Time history of liquid side mass fraction and temperature profiles. 
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present work has been ignored. 

Figs. 3.2(a,6) shows four dinstinct periods namely initial transient , first cP — law 
period, transition and second — law period. In the initial transient period, the rate of 
evaporation is low and the droplet surface temperature rises very sharply due to high heat 
input for droplet heating. As a consequence of increasing surface temperature, evaporation 
of the more volatile fuel ( heptane ) and the rate of surface regression increase rapidly. Since 
initial mass-fraction of heptane is very high and the initial surface temperature is much 
below the boiling point of dodecane, the fractional mass gasification rate of heptane (e/i) 
is very high during this period. 

After the initial transients, the first cP-law period starts. During this period, the droplet 
surface temperature remains close to the boiling point of heptane. The total evaporation 
rate and hence the surface regression rate exhibit quasi-steady behaviour.For heptane, 
surface mass-fraction and fractional mass gasification rate are high, indicating complete 
domination of the combustion of more volatile fuel. Very gradual decline in the surface 
mass-fraction of heptane is observed due to the high liquid side mass diffusional resistance. 

The transition period starts at a time when the surface temperature crosses the boiling 
point of heptane. During this time, most of the heat is consumed in heating up the droplet. 
Enough heat is not available to maintain high rates of evaporation and surface regression, 
and hence these quantities decrease with time. This is also indicated by flame contraction 
(table 3.1 ). The surface mass-fraction and fractional mass gasification rate of heptane fall 
steeply because of high diffusional resistance as well as a. decrease in the surface regression 

rate. 

During the 2”'^<P-law period, as the droplet temperature approaches the boiling point of 
dodecane, droplet heating becomes negligibly small and the surface temperature increases 
asymptotically. As the entire heat is used up for evaporation only, surface regression 
rate increases again and the droplet radius decreases exponentially. However, the total 
evaporation rate increases initially, and decreases later because of smaller droplet volume. 










time (sec) Flame Stand Off Ratio (Non-Dim) 


0.11 

0.23 

0.31 

0.35 

0.39 

0.42 

0.49 

0.55 

0.61 

0.66 

0.76 

0.89 

1.13 

1.254 


33.62 
35.91 
36.27 
36.02 
35.34 
34.04 

29.62 
28.10 
30.08 
32.20 
35.33 
38.07 
40.85 
41.26 


Table 3.1 Variation of Flame Stand Off Ratio With Time. 
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It can be seen that during the initial transient and the intermediate transient periods, 
liquid droplet heating is very important and can not be neglected. 

3.5.2 Effects of Volatility Differential 

Figures 3.3 (a — e) show the two d^-law periods very prominently when the volatility differ- 
ential is maximum (for the heptane- dodecane combination). The intermediate transition 
period also become narrower.Distinction of the two d^-law regions gradually decreases for 
smaller volatility differential and the transition period broadens. This is an indication of 
approaching the single component fuel burning behaviour in the limit. Fig. 3.3b shows 
that vigorous droplet heating is completed in the initial transient stage itself for smaller 
volatility differential. Fig 3.3d shows that the evaporation rate has a longer quasi-steady 
period with a decrease in volatility differential. So, on the whole, the transient as well as 
the multicomponent nature of combustion are reduced if the volatility differential is less. 

3.5.3 Effects of Initial Mass Fraction 

The effects of initial mass fraction on the burning characteristics are shown in Figs. 3.4(o — 
e) for the fuel combination of heptane and dodecane. It is evident from the figures that 
the first d^-law period is almost absent when the initial mass-fraction of the more volatile 
fuel is very low and except for the initial transient period, the second d^-law period is 
present all through. This is because the droplet does not possess enough more volatile 
component (especially close to the surface) to control the evaporation rate at any time. As 
the initial mass-fraction of heptane is increased, more heptane is available at the surface 
for evaporation in the initial stage. So, the first d^-law period increases and the second d"- 
law period decreases gradually. It is evident that the transition from heptane to dodecane 
burning is merged with the initial transient for small initial liquid side mass fraction of 
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Fig. 3.3 Effect of volatility differential on the transient variations of square of the dianae- 
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heptane. For higher values of initial liquid side mass fraction of heptane, the first cP-law 
period occurs distinctly and the transition is also clearly observed. 

3.5.4 Comparison with Earlier Work 

Fig. 3.5a presents a comparison for the liquid phase mass fraction profiles between the 
predictions of the present study and those of Shaw ( 1990 ). Mixtures of heptane-dodecane 
with different initial mass fractions have been considered. The comparison is made at 
a time when the surface mass fraction of heptane is same for both the cases. A good 
match is observed between the two studies. In Fig. 3.5b', the normalised burning rate 

prediction, [k = {k - A:octanol)/(^undecane “ ^octanol)] octanol-undecane mixture has 
been compared with the experimental mesurements of Law & Law ( 1982 ). Here also good 
agreement is seen, with minor deviations which may be attributed to the approximations 
invoked in the theoretical model. 

3.6 Conclusions 

A semi-analytical model has been developed here which is computationally much faster 
compared to the numerical model. This study shows very clearly the two d^-law periods 
in multicomponent fuel burning, as observed in experiments. It is also found that the 
multicomponent nature is predominant for large volatility differential and a high initial 
mass fraction of the more volatile component. 
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Fig. 3.4 Effect of initial mass fraction on the transient variations of square of the diame- 
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Chapter 4 

Numerical Solution Using Adaptive 
Grid Generation 



4.1 Introduction 


In the previous chapter, stress was given to obtain semi-analytical solution to see the fun- 
damental aspect of multicomponent single droplet burning in quiscent atmosphere. In the 
present chapter, some of the important assumptions of the previous chapter are relaxed, 
and a more thorough and detaild study of spherically symmetric nqulticomponent single 
droplet burning is undertaken through elaborate numerical analysis. Thin flame approxi- 
mation is relaxed, and the non-linear generation term (tj) is evaluated directly. Unity Lewis 
number condition for the gas phase has been relaxed. Adaptive grid generation technique 
has been used to track the flame movement throughout the history of droplet burning. The 
flame chemistry for droplet combustion is very complicated because, in reality, a very large 
number of chemical reactions take place before final products are reached. For all practical 
purposes, it is reasonable to assume a single step glob.al reaction for each of the fuel. The 
contribution of each of the fuel has been considered separately. In reality, above a certain 
flame temperature (Kuo, 1986), the CO 2 and H 2 O dissociate to form CO, O 2 and H 2 - 
Since they involve endothermic reactions, the flame temperature is not allowed to increase 
beyond a certain limit. This dissociation aspect of combustion has also been included 
in the present analysis. This study includes two extremely transient processes, namely, 
ignition and extinction. Because of the adaptive grid generation technique, it is possible 
to make a detailed study of flame movement, ignition, and extinction processes. Constant 
properties for both liquid and gas phases are retained in general. However, for one case, 
variable gas phase specific heat ( as a function of temperature ) was also considered. The 
present model can be extended to include the effect of flow, as well. 

For the liquid and gas phases, governing differential equations along with boundary 

/ 

conditions and initial conditions were discussed in section 2.4.5 and 3.2. Before discussing 
the full numerical solutions, these equations along with the boundary and initial conditions 

are first sunamarized for the sake of completeness. 



4.2 Liquid Phase Transport Equations 

For the problem attempted here, since there is no generation inside the droplet, the net 
mass transfer inside the droplet is zero. But, since the moving boundary (the droplet sur- 
face) ia always fixed at r=l, this gives rise to a convective term in the liquid phase transport 
equations. In absence of any flow field, the liquid pheise equations (2.60a-c) for the species 
and energy transport equations for spherically symmetric burning, can be rewritten as 
( a ) Species transport equation 


dYn 

dt 


+ 


( b ) Energy transport equation 


dt 


+ 


* \ 

-t) 

ah. 

dr 

-5- 

(4.1a) 

y) 

ahi 

dr 


(4.16) 

2aLe^ 

\ dt 

d^t ^ 

(4.1c) 

r > 

) dr ■ 

- = n 


4.3 Gas Phase Transport Equations 

Similarly, and using Ud = r^Vr (mass conservation), the gas phase species and energy 
equations (2.61a-d) are 
( a ) Species conservation equation 

dYt 


dt 

dYi 


dt 

dK 


'n fl pr\dYn D d f jdYfA _ 


+ Sc Ud 


(4.2a) 

(4.26) 

(4.2c) 


(b) Energy conservation equation 
BYt 


dt 


Sc Ud 


\r'‘^ f ) dr H3r\, dr) r“ar V, 


p D C, 


dr j 


(4.2i) 
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where, 


Qfi 


PoqDqoP 


LOji 


Qr--^ 


Q.= 


Poo^OQp 


■OJ. 


X 1 


Qt = - 


/ iv \ 


Poo^^ooP \i—i 


Y,h'^d:i]/Ah 


Yt= f C^dT 
Jt° 


(4.2e) 


and ujfi, djf 2 and uj^ are the rates of consumption for fuel 1, fuel 2 and oxidiser respectively 
in mole/m^ / s. The process of calculating these rates and also the term, are 

described in the next section. Here, /I, /2 and x stand for fuell, fuel2 and oxidiser 
repectively. 


4.4 Boundary and Initial Conditions 


Boundary Conditions 
(a) Droplet centre (r = 0) 


(b) Droplet surface (r = 1) 
(i)mass balance 

{ii)gas side species balance 




/I 


dr 

ay/2 


dr 


= 0 , 
= 0 , 


dr 


0 


ay/x 


dr 

r= 

ay/2 


dr 

r= 

dY^ 


dr 


.dri 

pui = -pui = -p-^ 


+ ^ {^n - Yit.) = 0 

+ ^{-Y.,) = 0 


r=l 


(4.3a) 

(4.36) 

(4.3c) 


(4.4) 


(4.5a) 

(4.56) 
(4.5 c) 
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{iii)liquid side species balance 

dVn 


dr 

dYf, 


dr 


aa Sc 
:: 

r=l ^ 

aa Sc 


(e/i - = 0 


^ ua oc / - \ 

+ (^/2 - Yj 2 s ) = 0 

r=:l 


where, aa = uj pj p. 
{iv)heat balance 


dT 

%= 


c ^ r 

= Pr pud 2_^ eiU + A 

r=l t=l 


dr 


r=l 


[v)phase equilibrium at interface 

Xfis = XfisQXp 


where, 




Xf2s = Xf2seyiV 

WnLnC,„ 


/^/2 


1 1 


Ru 


and, (3j2 — 


Thfx Ts_ 

Tbf 2 Ts/j 

Wj2L 


Ru 


(c) At outer boundary (r oo) 


(4.6a) 

(4.66) 

(4.7) 

(4.8a) 

(4.86) 

(4.8c) 


Initial Conditions 
(a) Liquid phase 



(4.9a) 

1/1=0 

(4.96) 

O 

It 

(4.9c) 

y — y 

* a: — ■* XyOO 

(4.9d) 


f(r,0)=f.„. 

(4.10a) 

f/i(r,0) = K„./i, 

(4.106) 

Yf2(r,0) = Y^j2 

(4.10c) 

T{r,0) = Tin 

(4.11a) 

V>i(r,0) = Fi,,/i 

(4.116) 

I/2(r,0) = y;„,/2 

(4.11c) 

K(»',0) = Yin,x 

(4.1 Id) 


(b) Gas phase 


4.5 Calculation of Heat of Reaction 


The combustion reaction of the two fuels fl and /2 can be represented separately as 
VfiMjx + I'rfiO-z = vcof\(-'0 + vachs\(^'Oi + i'h2J\H2 + + UH2onH20 (4.12) 

and, 


}/f2M f2 4" l^rf202 ^'COf2(^0 + ^^002/2^02 + i^H2/2-^2 + *^02/2^2 + ^H20f2H20 (4.13) 

where, i/ are the stoichiometric coefficients corresponding to different reactants ( Mfi , Mp 
, O 2 ) and products ( (X) , CO 2 , H 2 , O 2 ■, H 2 O ) and are indicated by the subscripts. 
Now, the r(*action rat<‘ for tiie reactants as well as for the products (cu,) can be evaluated 
from the overall r^^action rates (wi and U 2 for the present case) as 




-l/flUJl 




and 


£a>C 02/1 = -t^C02/l‘^l 

W/2 = ~*^/2^2 ‘^xf2 — ~^®/2^2 

<^C0f2 ~ ■"*^CO/2<^2 ‘^C02f2 — — ^00212^^2 

U2H2f3 *= ->^H2f2<^2 <^H20J2 — -^'^ 20 / 2^2 


(4.14a) 


(4.146) 


w/i and W /2 can be calculated from Eqs. 4.14a and 4.14b, and can be calculated as 




fVli + '^ 1 , (4.Uc) 

where, Wj is the overall reaction rate for the reaction given by Eq. 4.12 and d >2 is the overall 
reaction rate for the reaction given by Eq. 4.13. These reaction rates for the two Eqs. 4.12 
and 4.13 can be calculated from the relations 

wi(mofe/m®/sec) ~ (4.15a) 



( 4 . 156 ) 


and, 


CjAmok/ni^/sec) = A2T^%^f,Y^[C,2fexp[—^] 

RuT 


Hc'at of road ion for fiiois 1 and 2 ran be calculated as 


AKji - + ‘kfiK ^'(■oj\h,Qo ~ ^cOijihQOi ~ - ^'Oj/i^Oo ~ ^H2on^h 


02 ~ ^HiOnrijj^o 
(4.16a) 


and, 


^hrJ2 — /^/2^/2 + Vxf'ih.j. — UaOfihcO "■ ^'C’Oz/l^COz “ ^■W2/2^H2 “■ ^C>2/2^02 — ^i?20/2^H20 

(4.166) 

Now, 


N 


■Y.h% 


isl 


— — /iyjW/i + 6.®2W/2 + hluxjl + h°U}xf2 

+ ^WjO^WaO/l + ^B20‘^H20J2 + ^C02^C'02/1 + ^C02‘^002f2\ 

= — + h^f2‘^f2 + + ^ 1 / 2 ) 

+ ^H20 i'^HuOfi + WHjO/a) + ^co2 (^^^02/1 + WCO 2 / 2 )] 

=s — -- h^pUji 2 f ^2 — K i^'xfl ‘^1 + ^xf 2 ‘^ 2 ) 

+ ^%20 (^ftO/l^I + ^H2012^2) + ^C02 (*^002/1^^1 + *^ 002 / 2 ^^ 2 )] 

= “ ^H20J1^H20 ~ ^C02Sl^%02 

+ W3 [I'/aAja + ’^#/2^2 "" ’^H20j2h%20 - ^coifzhcoi 
— UiiAhrJl + <jj2^Kj2 

So, - Eili in J/m^fs, can be calculated by the expression 

’*•* hfujj — u}%Ahfi + 1^2 Ah f 2 (4-17) 

»=i 

where, wj and «2 are calculated in molefm^fs, and AKji and Ahr,f2 are calculated m 
Jfmole. 



4.6 Finite Difference Method of Solution for Liquid 
phase 

Id the present, i>robknn of droplet combustion, once the droplets are sparyed into the com- 
bustion chamber, the droplet surfac^e immediately starts regressing, giving rise to moving 
boundary at the liquid-gas interface. For the sake of computational analysis, when this 
moving surface is converted into a stationary surface by a suitable transformation, a con- 
vective term appears in the transport equations of the liquid phase, although there is no 
motion inside the droplet. For a finite difference solution this convective term creates 
problem if th<‘ whole differential equation is discretized using central difference scheme. 
To eliminat<‘ this difficulty, upwind scheme is applied for discretizing the convective part 
and central diffi'rencing scheme is used for discretizing the the diffusive part. Thus, by 
discretization of the liquid phase equations, a set of algebraic equations is obtained. This 
set of algebric equation.s are put into a tri-diagonal-matrix form to obtain the solution. 
The program, developed for solving the liquid phase, is quite general in nature. It can 
handle implicit, explicit and, (Irank- Nicholson schemes separately. For the present case, 
implicit scheme has been used considering the freedom for choosing the time step, At. 

4.6.1 Domain Discretization 

It is seen that the eqiiations are very steep near the droplet surface. Therefore, a non- 
uniform exponential grid was applied to discretize the domain inside the liquid droplet. 
This exponential distribution of the grid points was obtained by applying the relation ( 
Sanyal, 1991 ), 

r“Bexp{AS) + C (4.18a) 

where, 

n,ln(e,) jg ^ C = 1 - J5 (4.186) 

nr - 1 exp(/4 - 1) 

nr = number of radial elements, 

^ _ maximum grid length 
minimum' '^id'lenpK ’ 
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n = starting of the domain for liquid phase, 
r 2 = end of the domain for liquid phase, 

and 6 varies from 0 to 1 in uniform step size (l/n^). For the present problem, = 30 and 
Cr = 500 were chosen. 

4.6.2 Algebraic Equations 

The set of governing differential Eqs. 4.1 for the liquid phase and the corresponding 
boundary conditions given by Eqs. 4.3 through 4.9 are discretized as follows. 


= 0 m- 

“1 Ti 

n m-\ 

1 _ _ ._j 

-1 r = 


^ Ar 1 ^ 

/\ 

+ 

< 

N/ 



Ar = I [Ar_i + Ar+i] 

Applying the central difference scheme, the derivatives can be expressed as 


_ ^+1 ~ ^-1 
dr 2Ar 

m 


(4.19a) 


cPY 


Yr 


m+l 


Ar, 


+1 


Ar. 


+ 


+1 


Ar_i 




Yrr^-X 


Ar_i 


J_ 

Ar 


(4.196) 


and, applying the upwind scheme, we have 


dr 


Ym-l 


Ar_i 


(4.19c) 


Equation (4.1a) can be discretized by applying the central difference for the diffusive term, 

and upwind scheme for the convective term as follows. 


yp+i VP 

At 


dYfx 


m 


p+i 


+ ii-e) 


dYfx 


dt 


6 


'2dYn &^Yfx 
r dr ^ dr^ 




Sc Ud p r \ dYfi 
0 ] dr 


tP +1 


+ (l-») 


r dr dr^ 


D 


Sc Ud p r\ dYfi 

. P 7- J 
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— of- fyfl,m+l f I I \ ^ 

Llr 1, 2Ar ) + + a;-) 

^/l,m-l \ ^ 1 A fScUdpr\ (^}l,m — iyi,m-l) 

Ar,^)Arf^~[ } ) a;:” 

I n / ^1,771+1 ““ ^/l,m-l\ fyfl,m+l f I 1 

^l,m-l \ 1 1 f) — ( ^ (Xfl.TO — y/l,m-l j 

Ar_i j Ar J \ p j Ar_i 

Finally, rearranging the above expression, we arrive at the following form of discretized 
equation for the liquid phase species, /I, as 

+ C/.KC„ = + C'l/i>7i,„« (4.20a) 

where, the coefficients can be expressed as 


„ . [ -D D Scud p r 

-^/l = —OAt — T h -r — 1 T~7 

^ rAr ArAr_i p Ar_i 


. 6 At b ( 1 , ^At Scud pr 

_ q. yAr+i Ar.ij p Ar_i 

-gAtP [1 1 ' 

~ Ar r Ar+i_ 


^l,.=(l-«)A( -^ + 


-£) I) Scud p r 


Bln = 1 


(1 ->e)AtD ( 1 


[rAr ArAr_i p Ar_i J 
1 1 \ (1 - 0)At Sc Ud p 

Ar+i Ar_i y ” p Ar_i 

(l-^)Af P [1 1 

Ar r Ar+i_ 


(4.206) 


(4.20c) 


(4.20d) 


(4.20e) 


(4.20/) 


(4.2O5) 


Sinularly, the energy transport equation for the liquid phase can discretized and expressed 


in the form 


+ Crf^'i = AlrfU + B1 A + CItI^+i (4.21a) 


Where, the coeffients are expressed as 


\-aLe ale Scuipr 
At == —&At ArAr_i ^ Ar_i 


(4.216) 
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Bt — 1 + 


d Le 
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1 1 
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VAr+i ' Ar_i 
■6 At a Le 


AIt = (1-9) At 


Ar 
-a Le 


BAt Sc Ud p r 
p Ar_i 
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r Ar. 


+ij 

a Le Sc Ud p r 
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[ rAr ArAr_i pAr_i 


BIt 
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Ar 

CIt 


+ 


1 


(1 — 6) At Scud p r 


\Ar+i Ar_i 
(1 — 6) At 6c Le 
Ar 


p Ar_ 


1 1 
— ). _ — 
r Ar+ij 


(4.21c) ■ 

(4.21d) 

(4.21e) 

(4.21/) 

(4.215) 


The above discretized Eqs. 4.20a and 4.21a are very general, and can be used for either of 


the three schemes, explicit, Crank-Nicholson or implicit. Here, 
9= 0 for explicit scheme, 

9 = I for Crank- Nicholson scheme, 

9— 1 for implicit scheme. 


4.6.3 Discretization of Boundary Conditions 


(i) for droplet center (r = 0) 

The spherical coordinate system can not be applied at r=0. So, at r=0, the differential 
equations are discretized in rectangular coordinate system as follows. 

Equations 2.32a and 2.34 for species (f/i) and temperature (f ) can be expressed in 3- 

dimensional rectangular coordinates as 
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■ SYfi ^ dYn dY, 

dx dy dz 
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df 
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d^Yn d^Yn . ^ 
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dx"^ dy'^ 
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d‘2f d^f d^f 
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= 0 


= 0 (4.22a) 


(4.226) 


where, aa = p ujp. Equations 4.22a and 4.22b can be diacretized in a manner aimilax to 
that done on Eq. 4.1a in the previous section, and the simplified form can be written as 
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The ri§ht h^^nd side of the s-bovo 6 <^ii 3 >tioii is expressed 

R.H.S. 

ip+i 

— (1 — 0) Sc aa 


as 


— 9 Sc aa 




dx ' “ dy dz 


J m,n,/ 


■ dYf, , dYn , dYj^ 


ip 


+ eD 


d%i . d^Yn . d^Y 


ip+i 


dx dy dz 

^ J 771, n. 


+ 


n 


+ 


'/I 


[ dx^ ' dy^ ' dz^ 
= — 9Sc aa 


+ {l-d)D 


Tn,nJ 


s^Y,^ 

i ^T":; — t- 


dx'^ dy"^ dz'^ 


J m^Ujl 
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+ (A^;f (Kif “ 


Now, assuming equal grid spacing in all the three directions, we have Ax = Ay = A^. 
Due to symmetric condition at r=0, 


17 l,Tn-l,n,/ = l/l,m,n,/j l/l,m,n-l,i: — Yfl,m,n,l, l/l,7n,n,/-l — 1/1, m,n,/ ('^•23) 


Therefore, the above discretized equation for species fl at r— 0 becomes 
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+ (1-9)B 
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or. 


[(1 + 6q9) Yji,rn - 69 ^ 1 /i,m+i]'’'^' = [(1 + 6?^) l/i,m - eqeYji,m+iY (4-24a) 


where, 


DAt 


(Ax) 


(4.246) 


Therefore, the co-efficients for Eqs. 4.20a at r— 0, are 


69 DAt -6.9 D At 

Afl=0 B/>=1+ (A,.)— (Arf 


(4.25a) 
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^1,1=0 Bi;, = 1 - Cl,. = 


(Ar) 


(Ar) 


(4.256) 


and, those for Eq. 4.21a at r=0, are 


Ar = 0 = 1 + Cr = 


(Ar)' 


_ , 6 (1 - 9) a ie At 

_ . f\ 


AIt = 0 BIt = 1 ^ — - -- " 2 — “ 

(Ar)" 


(Ar)" 

6 {1 — 9) a Le At 


(Ar)^ 


(4.26a) 

(4.266) 


(iii) for droplet surface (r = 1) 

Discretizing the boundary conditions ( Eq 4.6a ) at r=l, we get the following discretized 
equation 

^l,m— 1 . Sc 


2Ar 


dCL lj C / \ _ 

H j— (e/i - i/i.m) = 0 


or, 


^/l,7n+l 


— ~ 2 


D 

Ar aa Sc 


D 


(^/l ^l.m) 


(4.27a) 


Similarly, discretizing the boundary condition ( Eq. 4.7 ) at r=l, we get 

- 2Ar Hi 


^m+l 


+ 2m- 1 


(4.276) 


Now, substituting the expressions for ^/i,m+i and fm+i from Eqs. 4.27a and 4.27b into 
Eqs. 4.20a and 4.21a, we get 


, „ V 2Ar aa 5c\ ‘ .1 

{Afl + C/l) yfi,m-l+[^n ■! g j ^rn 

/„ 2Ar aa Sc\ , 2Ar aa Sceji 2Ar aa ScCfi 
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< 


(Ar + ar)r^\+Br2I+^ 


(Air + C^lj) 2^-1 + BlyT^ 


2Ar Hi 


-Ct 


P+i 2ArBj„, 

+ ::: -Oly 


(4.286) 


The whole set of algebraic equations can be represented in the form of tri-diagonal matrices 

and solved through standard algorithm. 



4.7 Finite Element Method of Solution for Gas Phase 

F'inite Element Method is basically an integral approach to solve governing differential 
equations representing a given process in a given domain. For the present study, Galerkin’s 
Weighted Residual Method is used to discretize the gas phase equations. In this method, 
the entire domain is discretized into a number of neighbouring elements. By substituting 
the approximate numerical solution back into the differential equation, one gets a residue. 
In order to minimise the residue in an integral sense, the residue is integrated with respect 
to a suitable weighting function, known as the interpolation function, over each elemental 
domain and then it is put to zero. As an outcome of the integration performed by applying 
a suitable quadrature rule, a set of algebraic equations is obtained for each element. All 
these equations are then assembled into a global matrix. Then Ihe boundary conditions 
are added. Later this global matrix is solved by a suitable algorithm. 

For the present analysis, the global matrix takes the shape of a tri-diagonal matrix 
(TDM) due to one dimensional spherical geometry. This TDM was solved by a standard 

TDM solver. 

4.7.1 Domain Discretization 

The domain discretization for the gas phase solution is very important for the present 
problem. This is because finite rate chemistry has been taken into account while the ^ 
problem is modeled. Here, an adaptive grid generation technique is applied for tracking the^ 
flame movement. Fbr this, adaptive grid generation technique is developed as a separate: 
module. As the solution scheme marches in time direction, after each time step, the 
adaptive grid generation module is called separately to check the necessity of changing gn<^ 
point distribution. This is based on an error calculation method which is done statisticallyJ 
This process is is separately discussed later. , 



4.7.2 Finite Element Formulation for the Gas Phase 


ri 

r2 

1 

2 

4=-l 

^=+1 


In finite element formulation, the interpolation functions, also known as the shape func- 
tions, are generally selected from polynomial functions like Lagrangian, Hermitian etc. For 
iso- parametric elements, the shape funtion for the variable and the mapping frmction for 
coordinate transformation are taken to be identical. Here, the iso-parametric feature of 
the elements are retained and the shape functions (Ni) are chosen to be linear, for the sake 
of simplicity, as 

Wi = 5(l + «i), ■ = 1.2 (^-29) 


where, ( is the local coordinate, 

The global coordinate, r, can be transformed into the local coordinate , by the relation 

r- fn±r2.'\ 

L±J- (4.30) 




r2-ri 

2 


where, r 2 = global coordinates for points 1 and 2 respectively as shown in the sketd 
above. Following the above iso-parametric formulation, the variables and the coordinate 
linearly within the elements and can be expressed as 


vary 


T = X] nNi , Yfi = Yfi,i , F/2 = XI Yf 2 ,iNi , Yx — Yx,iNi and Yt — ^ YT,iI^ 

i=l .=1 •=! ‘=1 


where Yr = /J. CpdT and i=l,2. By applying Galerkin weighted residual technique, the gas 
phase equations are weight®! by the linear shape functions, Ni . Applying this technique, 
the Eqs. 4.2a through 4.2d can be written as 
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ir r^ dr = 0 (4.32c) 
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(4.32d) 


where, " indicates the previously iterated value. In Eq. 4.32a, the integrals containing the 
double derivative term can be further simplified by applying integration by p^lrts to those 
highest order derivatives and consequently the double dervatives can be reduced to single 
derivatives as follows. 

R.H.S. of equation(4.32a) 
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Therefore, 


- f N,Q,ydr - NivlD ^ = 0 

fssra ^ Ifssri 

Boundary conditions appear as a part of the finite element expression given by Eq. 4.33a. 
Similarly, Eqs. 4.32b, c,d can also be simplified to get the following relations. 
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Reduction of highest order derivative is achieved through weak formulation f 
of multiple integrals. The advantage of such a reduction is that the order of t e s p 
functions can be reduced upto which this need be complete. By assembling ell the elemen a 
matrices and-vectors into one global assembly , we arrive at the following relations for e 


expressions 4.33a, b,c,<i. 
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It is to be noted here that when all the elements are assembled globally, the derivatives at • 
r = 1 and at r = r 2 cancel out for all the internal elements. Only the the derivative at r = 1 
and r = Too remain for all the expressions, 4.34a through 4.34d. At r = 1, l/i, Yf 2 and Yr 
have specified boundary condition, and at r = r^o, all Yfi, i/ 2 , Y^ and Yt have specified 
boundary conditions. Therefore, for Yji ,yf2 and ix, the derivative boundary conditions 
at r=l and at r = Too and that for i^ at r = Too that arise in the finite element expressions 
4.34a,b,d are ignored, because they do not contribute to the assembled matrix. For i^, the 
derivative boundary condition at r = 1, can be written by applying Y^ = as 

= NtUjSc = [iViiVjujSc] (4.35) 

r=:l 

Using the above conditions and Eq. 4.35 the finite element expressions, 4.34a through 
4.34d can be written in discretized form as 


NiD 


dY^ 


dr 


[^mass] {Yfi,j} + [gcon] {T/ij} = {ggen} 

(4.36a) j 

[gmass] + [gcon] {Yf 2 ,j} = {ggen} 

(4.366) 1 

[^imass] {iij} + [gcon] {Y^j} = {ggen} 

(4.36c) 1 

[^mass] {Yt]} + [gcon] {ir/} = {ggen} 

(4.36d): 


where, gmass = Eii emass; gcon = Ef=i econ; ggen = E^i egen . 


For Yn 


emass - J NiNjr^dr ; egen 


econ = Sc Ud 


JN, 


1 + 


= J^Qfir^Nidr 

dNi dN. 


dr 


dr dr 


dr 


(4.37a)i 
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For Yj 2 : 


emass 


[N.N.r^dr ; egen — iQnr^N^dr 


econ = Scu,l « (l + t^^?^dr + D [ (4.376) 


For K 


emass 


J NiNjV^ dr ; 
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jQ.r^Ni 
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^ iV, (l + ^dr + D I + u,Sc (4.37c) 


emass = J NiNjT^dr ; 
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(4.37<f) 


The above integrations were performed numerically after transforming the global coordi- 
nate to the local coordinate. The local coordinate, ^ , varies from -1 to -j-1. For such 
numerical integration, Gauss- Legendre method was applied. Now, 

y*+l yJk 
t fi i ^ H i 


.. . == " . l lhL 


Therefore, Eq. 4.36a can be discretized in terms of explicit, implicit or Crank-Nicholson 
scheme as 

[ymass] -b ^[ircon] {K/ij }*'*'* + (1 — ^) [?con] == (4-38) 

where, $— 0 for explicit scheme 

=1 for Crank Nicholson scheme 
= 1 for implicit scheme. 

Rearranging Eq. 4.38, • 

[[ffmass] -b e Af[ 5 cH] ~ bcon]] {K/ij}^ + {i'i? At 

(4.39a) 



Similarly, the diKrrrtizcd equations for K/j, K and Yt can be expressed as 

[[gma.s,s] 0 At [geon]] = [[.gma.ss] - (1 - 6) M\gcon]] + {ggen] At 

(4.396) 

[[gma.s.s] -f 0 At [.gcon]] = [[grnass] - (1 - 6) At \gcon]] + {ggen] At 

(4.39c) 

[[grna.s.s] + 0 At (f/con]] = [[gmas^s] - (1 - 0) At [gcon]] [Yrj]^ + [ggen] At 

(4.39d) 

4.8 Iterative Scheme for the entire Solution Proce- 
dure 

The transport equations for the ga^i phase and the liquid phase are highly non-linear. 
The liquid phase and the gas phase equations are coupled together by clausius-Clapeyron 
equations. The coupling parameters at the interface are the evaporation velocity at the 
surface , heat transfer rate from the gas phase to the liquid droplet, Hi and the fractional 
mass gasification rate, ({. These parameters appear in both the liquid phase and the gas 
phase equations, and they are functions of time. For the iterative scheme, these three 
coupling parameters are adopted to converge upon. Due to the highly non-linear reaction 
term, underrelaxation has to be used in the later part of the combustion period. This 
iterative scheme is based on a two-step sequential process as depicted in Fig. 4.1. Both of 
these sequential processes come under the outer time-step loop ( indicated by the loop 100 
in Fig. 4.1 ) i.e. for each time step, these two processes are carried out one after another. 

In the beginning, the required input data are fed to the program and the initial grid 
distribution is chc«en from experience. With this arbitrary initial grid, it goes to the first 
module or, step. In the first step, first starting with an initial guess of u^. Hi and €fi , the j 
liquid phase solutions are obtained by solving Eqs. 4.28a and 4.28b in subroutine funda . 
Once the liquid phase solutions are obtained, the droplet parameters, T/i*, T/ 2 s T, 
at the liquid side become known. Now, using Clausius-Clapeyron relation ( Eqs. 4.18a,b,c 



), the gas side surface parameters, T, are calculated. Now, using the 

values at the surface as the given bound ary conditions, the gas phase solutions are ob- 
tained by solving the Eqs. 4.39a-d. Them from Eqs. 4.5a,b and 4.7, the coupling surface 
paramet('rs (f/,. c/ 2 , ) are calrtrlated from Eqs. 4.42a-d and they are compared to 

previous guess values for convergence. IS it fails to converge then the coupling parame- 
ters are updated by choosing a suitable xiiiderrelaxation factor. This process of solution 
is repeated in the next iteration loop 13 ( see Fig. 4.1 ) till convergence. When these 
values converge within the tolerance limit;, it moves on to the second module ( subroutine 
’adgrid’ ). In the second module, first error in each element is calculated statistically by 
substituting the converged solutions for «acli node into the energy conservation equation 
of the gas phase. T’hen overall global error is calculated by taking root mean square of all 
the elemental errors. If this global error exceeds the global tolerance limit, then refinement 
and derefinement of grid takes place according to the relative errors of the different ele- 
ments. After adaption it goes back to thie starting of iteration loop 10 and starts the first 
module once again. In this way, the secjnential, operation of stepl and step2 is repeated 
within a certain time-loop till the globat,! error satisfies the global tolerance limit. When 
this criterion is fulfilled, it moves on to the next time step loop. In this way the whole 
solution marches in time-direction. Tbee second module which covers the entire adaption 
process is given separtately in details in Fig. 4.2. 

4.9 Processing of Parameters after each Time-step 

(i) Calculation of droplet radius ( n) 

The droplet radius (rrf) , after cacl ’tiiue step, is calculated explicitly. 

mass of droplet at a given time (t) is 

xmass = mass of droplet at the previo'us time step (t-l) 

- amount of mass evaporated at tke previous step (t-l). 

This mass of the droplet at a given time is known explicitly. Therefore, from this relation 




Fig. 4.1 Flow Chart for the Overall Numerical Scheme 










llu' droplet radius, can be calculated as 


r-d = 


^ 3 xmass \ ® 
Att 0 } 


(4.40) 


Thus, droplet radius, can be calculated at each time step explicitly. 

(ii) Calculation of Ud, Hi, c/i, e/2 

The gas phase boundary conditions for fuell, fuel2 and temperature are given by Eqs 

(4.41 h) 
(4.41c) 
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+ ^(«/i-ni.)=o 
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solving these three boundary conditions ( Eqs. 4.41a, b,c ), 


Ud = 


D ( 


dYn 

dr 


lr=l 


.+ 


dY, 


dr 
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c/i = - 


Sc {Yfu + yf 2 s-l) 


UdSc dr 


Is 


r==l 


e/2 = 1 — c/i 


and, 


Hi= X 


dr 
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Pr pudY2 


(4.42a) 

(4.42fe) 

(4.42c) 

(4.42d) 
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Assuming a linear distribution of F/i, V /2 o.nd T at the gas phase boundary, 

dYn 


dr 

dYn 
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r=l 


r=l 


and, 
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irssl 


1/1,2 ~ XflP 

Ar 

1/2,2 ~ 1/2,1 

Ar 


r 2 -ri 

Ar 


So, after solving the gas phase equations the three surface coupling parameters (ud, tn and Hi 
can be retrieved from the three Eqs. 4.42a, b,c,d given above. 



4.10 Adaptive Grid Generation 

4.10.1 Introduction 


1 he solutions for the governing differential equations of certain physical processes like com- 
bustion etc., contains regions in which the variables exhibit high spatial activity like steep 
fronts and sharp peaks. In such problems, these high activity regions need to be resolved 
adaptively in order to ensure an efficient solution to the problem. This feature of sharp 
peaks is predominant in the present problem of droplet combustion. This adaption of grid 
points in the region of high spatial activity zone is based on an attempt to equidistribute 
positive weight functions. In implementing the one dimensional adaptive grid strategy, 
first the problem is solved for an initial coarse grid. Substituting this coarse grid solu- 
tions back into the equations, local discretization errors (Eg) are calculated statistically ( 
Smooke et al., 1989; Rao, 1992 ). These discretization errors form the basis for the positive 
weight functions. Then it is checked whether the total error satisfies the global tolerance 
limit (gtol), if not, grid points are added to the high error zones and removed from low 
error zones. Then, once again the equations are solved by using the newly adapted grids, 
and again the global tolerance limit is checked. This process goes on till this tolerance 
limit is satisfied. After each adaption, the coarse grid solution is then interpolated linearly 
onto the new grid. The interpolated results serve as an initial solution estimate for the 
iteration procedure on the new mesh. For the present problem, this whole process of mesh 
adaption is repeated at each time step as the solution proceeds in time. Here the low 
activity zone and the high activity zone keep on changing from time to time. So, refine- 
ment and derefinement both take place according to the need depending on situations and 
circumstances. i 

4.10.2 Formulation for Adaptive Grid Generation 

Fbr the present problem, temperature profile has been used as an indicator for locating the J 
high activity zone and the low activity zone. So, here, the gas phase energy Eq. 4. 2d has 
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been used foi calculating the elemental errors, mean, standard deviation etc. as a part of 
mesh adaption. 

After substituting the approximate numerical solution, f , into the energy Eq. 4. 2d, 
the elemental error, Eq can be calculated by integrating the square of the residue [rf\ over 
the entire element as 

-Ee = = ^ / r^dr (4.43) 

TT Jri 

where, the residue, Vg is given as ( see Eq. 4.32d ) 


dYj f I p r\ dYx a Le d 

r-e = ^ + 5'C Ud - + V ' 


dt 


p 


dr 


J.2 Qj, 




dr 


(4.44) 


, hg— length of the element, e and ' indicates approximate numerical solution and In the 
present problem, linear shape functions as indicated in the previous sections, are used. For 
such a case, = 0. So the residue, rg is calculated as 


dYx 

dt 




(4.45) 


where, both the spatial and time derivatives are calculated as ( using Eq. 4.31 ) 


dYx 


dr 


' a a 




dr 

rri f 




and. 


dYx 

dt 


dt 

= Ni 


^ {N,Yn+N,Yn) 


'V-A'+l V^'' 

^Tl ~ ^Tl 


At 


+ N2 


'yK+l ■yK'' 
^T2 ~ -‘T2 

At 


where, 6g~ 0 for explicit scheme 

=1/2 for Crank-Nicholson scheme 
=1 for implicit scheme 
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The total global error {t}) is calculated by adding all the elemental errors and then 
taking square root of the summation as 


V = 


■ N 

El. 

Le=l 


, 1/2 


where, N — number of elements. 

The mean(/i), standard deviation (cr) and Jmax are calculated as 


(4.47a) 


N 2 

e=l ■''' 


(4.476) 


[E" . (E, - m)' 


1/2 


E, 


N-1 

in ax -£'7n,in 


J ' 7ft, vn , -j 

max — r i 


(4.47c) 

(4.47d) 


where, B,nax— maximum elemental error and, 

Emin— minimum elemental error. 

Here, Jmax is an indicator such that {Emax— Emin) < cr Jmax- H indicates that the elemental 
error. Eg cannot exceed {fJ.-\-Jmax cr). The algorithm for the adaption process is represented 
in the form of the flow chart given in Fig. 4.2. 


4.10.3 Testing of Adaptive Grid Algorithm 

Adaption of grid points are based on the energy equation for the gas phase ( see Eq. 4.2d 
). In order to test the adaptive grid generation algorithm, an analytical solution is chosen 
of the form 

Yt + [C 2 {r - 1) - a f] (4.48) 

where, Yf is a function of r and t and ci, C 2 and a are constants. Depending on the values 
of C 2 , the steepness of the analytical solution varies. So, this analytical solution containing 
tangent hyperbolic function, can be used for testing the adaption over a wide range of 
spatial steepness. In this method of assuming the analytical solution for the gas phase 
energy equation, the source term is calculated by replacing the analytical solution back 
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Fig. 4.2 Flow Chart for Adaptive Grid Generation 
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into the original governing equation as 
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dYi 

dt 


'2 oc Le 


~T + 


1 or\\ dYi 


d'^Yr 


dr dr^ 


(4.49) 


By using the expression for the assumed analytical solution, we get 
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dYx 

dt 


Sec^h [c 2 (r — 1) — a f] 


(4.50a) 

(4.506) 

(4.50c) 


Therefore, the corresponding source term for the assumed analytical solution is calculated, 
by using the above relations, as 
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2 a Le cj tank (xx) 
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1 N 

1 p r 
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se(?h (xx) (4.51) 


where, xx = C 2 (r — 1) — a t 

So, if one substitutes the source term given by Eq. 4.51 into the governing Eq. 4.49, then 
the assumed analytical solution for the gas phase energy equation remains valid and the 
effective governing equation becomes 


dYr 
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Sc Ud 
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d'^Yr 
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— Sc Ud 
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The authenticity of the grid generation can be checked by comparing the solution of the 
gas phase energy equation obtained through mesh adaption, using Eq. 4.39d, with the 
solution obtained from the analytical expression given by Eq. 4.48. Using this process, 
only the discretization error can be checked, so, this process gives a fairly good verification 
for the adaption of grids. Neither the boundary conditions, nor the original source term 
can be checked by this method. 



4.10.4 Results of the Test Problem 


Based on the analytical solution ( see Eq. 4.48 ), choosing suitable values of Cj, C 2 and a, 
various checks are carried out. Basically the goals of verification are to check 

(i) the efficiency of refinement and derefinement, and 

(ii) whether or not the solution is marching properly in time direction or not. 

Case I : Efficiency of refinement and derefinement are cheked in two ways. Figs. 4.3a-d 
shows the first process. Firstly, the adaption process is started with only 2 elements. Then 
through refinement and derefinement, the adaption process proceeds forward in a number 
of steps to the final grid configuration depending on the tolerance limit. Figs. 4.3a-d show 
that this adaption process is very smooth in reaching the desired condition. For this case, 
it is observed that after the 7'^ adaption, the numerical solution and analytical solutions 
are very close to each other. 

Secondly, the adaption process is started with a large number of uniformly distributed 
elements (30 elements) ( see Figs. 4.4a-d ). Here also, it is observed that during adaption, 
grids are removed from the low activity (spatially less steep) zone and, are added to the high 
activity(spatially very steep) zone. This removal and addition takes place very smoothly 
as seen in the figure . 

In both of the above situations, it is observed that the numerical solution perfectly 
matches with the analytical solution. This indicates that spatial discretization of the 
governing differential equations is fairly accurate. So, this case demonstrates the efficiency 
and smoothness with which the present algorithm works in spatial coordinates. 

CaseTI : This case ( Fig. 4.5 ) demonstrates the efficiency of the algorithm regarding 
how it moves in time direction. These results indicate that the transient part of the 
discretization is also fairly accurate. We can see from Fig. 4.5 that the numerical results 
and the analytical results match nicely with each other at all time steps. 

In verifying the above conditions, the values of the constants were taken to be Cj = 
1.0, C 2 = 0 5 and a == 1.0. The above conditions were also checked for different steep- 
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Fig. 4.3 Evolution of grids at a given time starting with only two elements. 
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Fig. 4.4 Evolution of grids at a given time starting with only thirty elements. 
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Fig. 4.5 Evolution of grids at different times. 
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ness of the analytical tangent hyperbolic solution. This was done by choosing = 

0.9, 0.7, 0.3, 0.1, 0.01. Similar matching was observed for all the cases. 

4.11 Closure 

In this chapter, both liquid gas phase transient equations are discretized for numerical 
solution for the case of spherico-symmetric burning of multicomponent fuel droplet. Liquid 
pha.se equations are discretized by Finite Difference Method ( FDM ) and the gas phase 
equations are discretized by using Galerkin Finite Element Method ( FEM ). Adaptive grid 
generation based on local elemental error calculation was used to determine the position 
of grid points. Gas phase energy equation was used as the basis for error calculation. 
The adaptive grid generation technique is tested with a known analytical solution which 
is tan hyperbolic in nature. Discussions of the results for these discretized equations are 
presented in the next chapter in details. 



Chapter 5 

Results and Discussions 
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5.1 Introduction 


This chapter, throws light on various important features of multicomponent single droplet 
fuel combustion. The elaborations include both liquid phase and gas phase profiles for 
the temperature and mass-fractions at different times throughout the history of droplet 
burning. They show how these profiles evolve initially during the highly transient ignition 
period, how they develop into the intermediate period ( first <P — law period, transition 
and second cP — law period ) comprising the major portion of the burning history, and 
finally how they move towards extinction. The structure of the flame is closely looked 
into taking the help of grid generation technique to track down the flame movement. In 
multicomponent fuel burning, the difference in volatility of the two fuels plays an important 
role. The effect of volatility differential on different burning characteristics has been looked 
into. How does any change in the initial mass-fraction affect the nature of burning, is 
also included in the present discussion. Several other important features, such as effect 
of dissociation, ambient temperature, variable and constant gas phase specific heats, gas 
phase and liquid phase diffusion coefficients etc. are also included in the discussion. Finally 
the comparison with the semi-analytical model as well as with existing experimental results 
are also made. The property values are used from different sources such as Reid et al. ( 

1988 ), penner (1957) etc. 

5.2 Surface Parameters 


This section discusses the transient variation of the surface parameters namely surface tem- 
perature, T., total volumetric gasification rate, V, flame stand-off ratio, r/, flame position 


from droplet center, r^, flame position from droplet surface, heat coming from the gas 
phase, heat input to the liquid droplet, Tf,, heat of evaporation, ifev, factional mass 

gasification rate, e,., liquid aide and gas side mass fractions. Y,u„ and droplet 

radius, r, and square of droplet diameter, with time. Here, heptane and hexadecane 
fuel mixture with an initial liquid side mass-traction combination of 0.9-0.1 is chosen for 
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discussing these parameters. It is seen from the figures that there are six distinct periods ' 
that cover the entire burning history namely (i) ignition delay period, (ii) ignition period 
(iii) - law period, (iv) transition ( from first d? - law period to the - law period 
) Period, (v) 2d— law period and finally, (vi) the extinction period. The extinction 
period is basically a volatility limited period during which the more volatile component is 
rapidly depleted from the entire droplet composition. After the droplet is introduced into 
the combustion chamber, there is a time lag before ignition takes place. This is the ignition 
delay period of droplet burning. Starting from the beginning, this period exists approx- 
imately up to 0.0015s, which is about 2.28 % of the total burning time for the present 
case and is 0.0013 s for the burning of heptane- do decane fuel with an initial mass fraction 
combination of 0.9-0. 1. During this period, the peak temperature starts rising very slowly 
at around ry=:4.2721, due to very slow reaction rate. Here, curves are not steep. Since the 
flame is almost absent, this period falls under the category of pure evaporation. 

At the end of this period, at r/=4.2721 the reaction starts gathering momentum and the 
initial transient marked by ignition, expansion of the flame towards the droplet surface and 
vigorous droplet heating sets in. As this initial transient sets in, the reaction at the flame 
region becomes vigorous and the flame temperature rises almost instantly to 3223 K at t = 
0.00166967s. At t=0.002s, heat coming from the gas phase ( Hg ), heat going into the liquid 
droplet ( Hi ) and the heat available for evaporation ( ) all reach their peak values. At 

this point, 65 % of the total heat available from the gas phase goes for droplet heating and 
the rest ( ) is available for fuel evaporation. Then the percentage of the heat available 

for the evaporation gradually increases and that available for droplet heating gradually 
decreases. As the flame temperature rises, the flame gradually moves and expands towards 
the droplet surface from r/ = 4.2721 to r/ = 3.31 because during this initial period, fuel is 
available only at the droplet surface. Therefore, the temperature gradient at the gas side 
of the droplet surface becomes very high. So, the droplet heating becomes very vigorous 
due to this high temperature gradient at the droplet surface. This results in a sharp rise 
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in surface temperature ( Z )■ As heat of evaporation, becomes high and the surface 
temperature rises sharply and approaches the the boiling point of heptane(371K), the 
evaporation of heptane also increases vigorously. Thus, the total volumetric gasification 
rate increases sharply. The liquid side surface mass-fraction, T/i/j, of the more volatile 
fuel is very high due to it s high initial value, but it falls steeply due to high diffusional 
resistance and low surface regression rate. Due to the high value of liquid side surface 
mass-fraction of the more volatile fuel(heptane), it’s mass-fraction on the gas-side, Yjug 
also increases very sharply. But, since the surface mass-fraction of the less volatile fuel on 
the liquid side is very small, it’s gas-side value(Y/ 2 s 5 ) remains very small. In this period, 
mainly the more volatile fuel(heptane) contributes to the total gasification rate, V as the 
surface is enriched with heptane and the temperature is much below the boiling point of 
the less volatile fuel ( hexadecane ). So, the fractional mass gasification rate, eji, for the 
more volatile fuel is also very high. Though, initially the flame is situated very near the 
droplet surface ( r/=3-31 ), but immediately after the flame is established, sharp increase 
in the gasification rate, V, causes fuel vapour accumulation between the flame and the 
droplet surface, and as a result, the flame starts moving away from the surface as indicated 
in Fig. 5.1c by r}. The initial transient period spans for about 0.0015 s, starting at 0.0015s 
and ending at 0.003s approximately. It covers 2.285 % of the total burning time. 

After the initial transient is over, - law period starts around 0.003s. This period 
covers nearly 35.63 % of the total burning time and is maintained till 0.2625s approxi- 
mately. During this period, as the flame moves away from the droplet surface, the gas 
side temperature gradient at the droplet surface gradually decreases and becomes almost 
steady. Therefore, the heat available from the gas phase also decreases, initially very fast, 
then gradually and finally becomes almost steady . The liquid side surface mass-fraction 
of the more volatile fuel, Y/Us, decreases at a constant rate. This decrease is due to high 
liquid-side diffusional resistance for which enough heptane cannot reach the surface for 
maintaining it’a surface mass-fraction. And, also, since the droplet radius regression rate 
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during the I'V - law period is neither too high nor too low, the slope of the Y,„, curve is 
maintained at an intermediate value. During this period, the gas side surface mass-fraction 
of heptane remains high due to it’s high liquid side initial value and consequently, the sur- 
face temperature remains in the vicinity of the boiling point of the more volatile fuel. In 
the initial part, Hi falls very fast and then decreases gradually and finally becomes almost 
steady. At about t= 0.015s, heat going into the droplet reaches a minimum and about 
69 % of the total heat available from the gas phase goes for droplet evaporation. The 
amount of heat available for evaporation, Hev, remains nearly constant. And since, the 
droplet surface area gradually decreases, the total gasification rate decreases quite steadily 
as shown in Fig. 5.1b. Due to the relatively large amount of the more volatile fuel present 
at the droplet surface, the fractional mass gasification rate of the more volatile component, 
Cfi, remains quite high and is maintained almost at a constant value during this period. 
The gas-side mass- fraction of the less volatile fuel remains low due to it’s low initial value. 
As the overall gasification rate decreases steadily, the fuel vapour accumulation becomes 
progressively lesser towards the end of this period. So, the flame movement slows down 
steadily. Because of the fast flame movement in the initial part of this period, the flame 
stand-oflf ratio remains quite steep. But, towards the end of this period, the flame stand-off 
ratio increases at a much slower rate more due to the slow rate of drplet radius regression. 

Due to the liquid side mass-diffusional resistance and low regression rate, the liquid 
side surface mass-fraction of heptane falls continuously and when it falls below 0.10325, 
it’s gas side surface mass fraction could not be maintained any more around 0.84 and 
eventually, it starts falling and the gas side surface mass fraction of hexadecane starts 
rising. Therefore, as a result of the above surface change, the transition from one cP — law 
period to the next starts which causes the surface temperature to cross the boiling point 
(371K) of the more volatile fuel. The time of such occurrence depends upon the liquid 
side mass-fractions of the fuels at the surface as mentioned above. For the present case, 
the transition starts approximately around 0.02575 s when Y/iag, F}i/s and T, are around 
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0.79, 0.10325 and 374 K respectively, and it ends around 0.039 s approximately when the' 
values of the above parameters become 0.57, 0.02 and 465 K. This transition period covers 
about 19.5 % of the total burning time. In this period, the gas side temperature gradient 
at the surface again increases resulting in an increase of total heat flow ( Hi ) from gas 
side. And also, the droplet heating is again increased as the droplet surface temperature 
makes a move towards the boiling point of hexadecane as indicated by the curve for Hi 
and Ts- At about t— 0.0315 s ( about half way through the transition ), the total heat 
coming frrom the gas phase ( ) as well as the heat going into the liquid droplet ( Hi ) 

reaches their maximum values and consequently, heat going for evaporation ( H^v ) reaches 
it’s minimum value; 56.2 % of the available heat goes for droplet heating. It is seen from 
Fig. 5.1b that in the first half ( from 0.02625 s to 0.0315 s ) of the transition period, the 
heat available for evaporation gradually falls. Therefore, during this period, the gradient of 
the gasification rate curve, V, increases relatively, indicating a decrease in the gasification 
rate ( Fig. 5.1b ) and consequently, surface regression rate also falls as indicated by the 
decrease in gradient for the curve of Fig. 5.1a. But, in the second half ( 0.0315 s to 
0.039 s ) of the transition period, the heat of evaporation again increases and also, the 
surface temperature approaches more towards the boiling point of hexadecane, thereby 
increasing the gasification rate for hexadecane as indicated by the curve for Yj 2 sg in Fig. 
5.1a, Therefore, the total gasification rate is increased during the later part of transition. 
So, the surface regression again picks up momentum as shown by the - law curve in 
Fig. 5,1a. During transition, fuel vapour accumulation becomes very near to zero. At this 
stage, the rate of fuel consumption at the flame is nearly equal to the rate at which it is 
being evaporated at the droplet surface. This transition from — /au) period to the 
- law period is gradual and is caused by the high droplet radius regression rate. The 
latter ( i.e. high droplet radius regression rate ) maintains the surface mass-fraction of the 
fuels by exposing the inner core to the surface, depending on their initial values. This also 

nullifies the effect of high diffusional resistance. 
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After the transition from - law period to the 2-^d? - law period is over, the 
second d-square law period starts at around t=0.039 s and it continues till the end of 
droplet burning. For the present case, the second cP - law period covers approximately 40 
% of the total droplet lifetime. As the droplet surface temperature rises above the boiling 
point of heptane ( during transition ) and eventually reaches an asypmtotic value, the 
droplet heating gradually decreases and becomes negligibly small. This asymptotic value 
of surface temperature depends upon the liquid side surface mass- fraction combination of 
the fuels. Most of the heat available from the gas phase is used up for evaporation only. 
At t=0.039 s ( i.e. at the start of the second (P — law period ), 58 % of the available heat 
goes for evaporation and it increases upto 100 % near the end of burning. So, the surface 
regression rate is maintained till the end. The gasification rate initially increases due to 
the inertial effect of the later half of the transition. After 0.04275 s, it falls initially less 
steeply, but later on, the steepness increases owing to the small droplet volume, though 
the heat for evaporation goes on till the end. The droplet radius also decreases very fast 
due to the insufficient volume of the droplet. Results show that only 18.5^ % of the initial 
droplet volume and 43 % of the initial surface area remains at the start of the second 
(P - law period. They become 5.7 % and 14.8 % respectively at t=0.057 s when both the 
gasification rate and the radius start to fall very steeply. During this period, the flame 
remains steady, maintaining almost a fixed distance from the droplet center. The flame 
stand-off ratio increases with lesser steepness during the initial part of the period indicating 
that the droplet radius fails less steeply. But later on, the flame stand-off ratio becomes 
very steep owing to the very high rate of fall of the droplet radius. During second dP — law 
period, the species are transported to and exposed at the regressing surface approximately 
at a constant rate. In this period, the droplet radius decreases so fast that both liquid side 
and gas side mass-fractions of heptane are not only maintained but, actually they increase 
very gradually ( Fig. 5.1a ) as a result of the droplet core being exposed to the surface. 
However, the liquid side and the gas side mass fractions of hexadecane fall very gradually 
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due to it’s lower initial liquid side value ( 0.1 ). The distance between the curves r', and 
r'j indicates the instantaneous droplet radius. 

At the end of the 2""d^ - leu, period, the highly transient phenomenon of extinction 
begins. This extinction starts at about 0.06648881 s, when the volume, surface area and 
diameter of the droplet becomes approximately 0.0257 %, 0.4044 % and 6.36 % respec- 
tively of their initial values. This is marked by a sharp decrease in flame temperature and 
is finished in approximately 260 fis s. Therefore, the time gradient for the flame tempera- 
ture like all other parameters e.g. droplet radius, overall gasification rate, both liquid-side 
and gas-side mass-fractions etc. become, practically, infinite during the extinction period. 
The flame also starts moving towards the droplet in search of fuel at an almost infinite 
rate. The droplet radius becomes very small at this stage. During this final period ( which 
may also be called volatility limited period ), the more volatile component ( heptane ) is 
rapidly depleted from the entire droplet composition. This happens because at this stage, 
the droplet size (less than 2-3 microns) becomes comparable to the characteristic mass dif- 
fusion lenth. Therefore, the liquid side mass diffusional transport becomes effectively very 
high causing perpetual uniformization of concentration. Thus, the volatility differential 
becomes the rate controlling factor for evaporation at the droplet surface. So, the surface 
temperature rises sharply, approaching the equilibrium gasification temperature of the less 
volatile component ( hexadecane ). This phenomenon was referred by Law ( 1982 ). 

5.3 Liquid-Side Mass-Fraction and Temperature Pro- 
files 

This section describes the time history of mass-fraction and temperature distribution within 
the liquid droplet for heptane-do decane combination with 0.9 initial mass-fraction of hep- 
tane ( see Figs. 5.2a and 5.3a ). Due to the non-dimensionalization of r-coordinate with 
respect to the instantaneous droplet radius, the droplet surface always remains at r— 1. 
Thus, the r-coordinate shrinks as the droplet surface regresses. 
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Figure 5.2a shows that as soon as the the droplet is introduced into the hot environment, 
a concentration boundary layer is established near the droplet surface on the liquid side. 
This happens because the more volatile component evaporates very rapidly from the droplet 
surface, and because of liquid side diffusional resistance, it cannot be supplied fast enough 
from the inner core to the droplet surface. Due to this , the liquid side surface mass- 
fraction of the more volatile fuel falls continuously and the concentration gradient at the 
liquid-side surface becomes progressively steeper. This happens up to 0.04 s ( indicated by 
curve 10 ) which is just before the beginning of — law period. During the — law 
period, the concentration gradient at the liquid-side boundary increases very very slowly 
and remains almost constant as indicated by the curves 10, 11, 13, 14 and 15 of Fig. 5.2a. 
This happens because the liquid side mass-fractions of the fuels are almost maintained at 

constant values due to high radius regression rate. 

It is observed from figure 5.2a that after 0.065 s i.e. towards the end of combustion, 
the concentration boundary layer falls apart and the concentration gradient at the droplet 
surface becomes lesser and lesser and, the liquid side concentration becomes perpetually 
uniform. As discussed in the previous section, this uniformization of concentration happens 
because mass diffusional transport inside the droplet becomes highly enhanced owing to 
the small droplet size towards the end of combustion. This is corroborated by Fig. 5.2b 
which shows that the mass fraction of heptane at the droplet centre is maintained around 
it’s initial value ( 0.9 ) till 0.065 s which is very near to the end of burning. Then the value 
( mass fraction of heptane at the droplet centre ) starts falling very fast, almost with an 


infinite gradient. 

It is seen that the temperature profiles ( see Fig. 5.3a ) do not approve the existence 
of any steep thermal boundary layer. This happens because of relatively high thermal 

diffusivity which causes the heat to be conducted very fast to the core. This is in contrast 

r i’ .....filrao wlif>rp the steep concentration boundary layer 

to the nature of the mass-fraction profile 5 • 

is present almost from the 
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Two distinct zones, corresponding to the two (P — law periods, are observed from Fig 
5.3a. In the initial transient, which exists approximately upto 0.005 s ( curve 3 ), the 
surface temperature makes a move towards the boiling point(371 K) of the more volatile 
fuel(heptane) and it remains near this temperature till the end of the (P — law period 
i.e. almost upto 0.02 s ( curve 6 ) and more heat is used up for fuel evaporation. During the 
Id— law period, the heat from the surface region gets enough time to be conducted to 
the droplet core. So, the temperature profiles become flatter compared to those prevailing 
during the initial transient. The second transient starts at the end of of the - law 
period ( between curves 6 and 7 ) and exists almost upto 0.04 s( curve 10 ) . During this 
period, the surface temperature rises above the boiling point (371K) of the more volatile 
fuel (heptane) and the surface is heated up more, leaving behind the droplet core. The 
period of the second transient is longer and less steep compared to the initial transient. 
This happens due to high radius regression rate for which surface mass-fraction of heptane 
at the surface is maintained resulting. in a gradual shifting of surface temperature from the 
boiling point ( 371 K ) of the more volatile component (i.e. heptane ) to a higher value 
beyond 371 K depending upon the liquid side surface mass-fraction combination of the 
fuels. This does not happen during the initial transient, as then the radius regresses very 
slowly. After the second transient, the 2"*^ <P — law period starts at around 0.04 s ( curve 
10 ) and it continues till the end of burning. During this period, again heat goes from 
the surface to the core for heating it up, whereas the surface temperature hardly changes 
because of almost constant values of gas and liquid side surface mass-fractions of the fuels 
(see also Fig. 5.1a). 

It is seen from Fig. 5.3a that after 0.065 s i.e. almost at the end of combustion, 
the entire droplet surface temperature starts rising. During this rise of droplet surface 
temperature, no gradient is observed near the droplet surface. This means that the droplet 
core is also heated up instantaneously along with the droplet surface. This instantaneous 
uniformization of droplet temperature happens due the very small droplet size ( less than 
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Fig. 5.2a Time history of mass-fraction profiles for heptane inside the fuel droplet for 
hcptane-hcxadecane burniiigi with Kn,/i= 0.9. 



Fig. 5.2b Transient variation of heptane mass fraction at the liquid droplet centre for 
heptane-hexadecane fuel burning with Jin, /I = 0-9 
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2-3 microns ) towards the end of burning. Owing to the small droplet size, heat transfer ‘ 
rate becomes highly effective and causes instantaneous uniformization of temperature. 
This effect is also seen from Fig. 5.3b where towards the end of droplet burning, the 
temperature of the droplet centre is seen to increase almost instantaneously. Figures 5.2b 
and 5.3b indicate that the effect of multicomponent nature of droplet burning is strong 
enough to the influence the temperature but not the concentration at the droplet centre 
for most of the droplet lifetime. The present work has ignored the possibility of micro- 
explosion. It is expected that if the degree of superheat for the more volatile fuel is high 
then micro-explosion may take place. 

5.4 Gas-Side Mass-Fraction and Temperature Pro- 
files 

Figures 5.4a-d, 5.5a-c and 5.6a-c depict the time history of mass-fraction profiles for the 
gas-phase species and also the temp'erature profiles for the same phase. For all theses 
figures, non-dimensional radial distance has been used as the abscissa. So, the movement 
and the position of the flame region indicated in all these figures are always in comparison 
to the movement of the droplet surface. However, Fig. 5.4d indicates both the actual 
position and movement of the flame region in relation to the droplet center and the flame 
stand-off ratio. 

Figure 5.4a explains how the temperature profile evolves during the initial transient 
period. In this initial period of burning, since the droplet surface hardly regresses, non- 
dimensional and normalised radial distances ( i.e. dimensional radial distance divided by 
the initial droplet radius ) are almost identical. Therefore, the flame movement as shown 
in Fig. 5.4a is almost identical with the actual flame movement. It shows that the ignition 
takes place very near the droplet surface at about r=4. Initially, the entire gas phase is 
maintained at 1000 K. The flame region is indicated by the steep temperature gradient 
region. After the ignition starts, the chemical reaction at the flame becomes vigorous 
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and as a result of the enormous amount of heat released during the chemical reaction, 
the temperature of the reaction zone becomes progressively higher. But, after some time, 
flame temperature reaches more or less a steady value. This is because, the heat evolved 
in the reaction is balanced by the amount of heat being conducted away on either side of 
the flame region as well as by the heat consumed by the dissociation reactions. When the 
flame temperature becomes high, the products like carbon dioxide and water dissociate 
into carbon monoxide, hydrogen and oxygen and form an equilibrium mixture depending 
on the temperature. These dissociation reactions , being endothermic in nature, do not 
allow the flame temperature to exceed a certain limit. The present model does not consider 
the effect of radiation and convection. It ( the present model ) uses both constant and 
variable specific heat in the gas phase. Figures 5.4a, b, c & d are for constant specific heat 
which also results in a slight overprediction of the flame temperature. The present results 
show the flame temperature to be as high as 3000 K. Whereas, for variable specific heat 
case, the flame temperature is predicted to be around 2950 K. 

During the initial transient, the amount of fuel vapour present, relative to the oxidiser, 
is quite small due to low evaporation rate. So, while the flame is being fully established, 
the flame region moves and expands towards the droplet surface in search of more fuel to 
sustain combustion. Thus, the gas-side temperature gradient at the droplet surface be- 
comes progressively higher. This enables vigorous droplet heating. So, during this period, 
the droplet surface temperature increases very rapidly, and as a result mass-gasification 
rate increases very fast. This results in fuel vapour accumulation between the flame and 
the droplet surface as a result of which the flame gradually moves away from the surface 
( Fig. 5.4a and 5.4d ). After the initial transient is over, the flame movement slows down 
and almost stops with respect to the droplet center after the transition period and remains 

in that region almost till extinction ( see Fig. 5.4d ). 

Figure 5.4c sheds light into the extinction phenomenon. It is a highly transient process. 

The figure shows that the flame gradually moves towards the outer boundary of the solution 



129 


domain where the free stream conditions are forced though the flame actually remainiUe ' 
same region w.r.t. the droplet center. This happens because the r-coordinate shrinks along 
with the droplet surface, although the flame remains almost at the same position. As the 
extinction starts, the sharpness of the flame progressively reduces and gradually the peak 
temperature comes down. 

Figures 5.4a,b,c also reveal the efficient adaptive nature of the grid generation technique. 
They show that the grids are removed from the low temperature gradient zone and are 
added to the high temperature gradient zone as the flame front moves one position to 
another. In this way, the tracking of the flame movement is done quite efficiently for the 
present study. 

Figures 5.5a-c show the mass-fraction profiles for heptane and oxygen. It is seen from 
Fig. 5.5a that the mass-fraction profiles for heptane are initially very steep near the droplet 
surface. The steepness at the flame becomes initially flatter ( curves 1, 2, 3 ) but then 
as the flame expands towards the droplet surface, this steepness at the flame increases ( 
curves 3, 4, 5 ) and then, finally, as the flame moves away from the droplet surface due to 
fuel vapour acummulation, the steepness at the flame falls gradually ( curves 5, 6 ). The 
mass fraction profiles of heptane at the droplet surface are initially quite steep ( curves 
1, 2 ), then becomes a little bit flatter for a small time ( curves 3, 4, around 0.001647 
s ), then again steeper ( curve 5, around 0.00174 s ) and then, finally, becomes flatter ( 
curve 6, around 0.005 s ). This is also evident from the curve of liquid side surface mass 
fraction for heptane in Fig. 5.1a. In case of oxidiser mass-fraction profiles, it is seen that 
initially ( when the flame is not fully established ), the profiles are flatter ( curves 1, 2 ), 
but, as the flame is fully established ( curve 3 and 4, around t 0.0016329 s ), steepness of 
oxidiser mass-fraction profiles at the flame region also increases till the flame starts moving 
outward. From Fig. 5.5b, one can see that as the flame moves away from the droplet 
surface, the mass-fraction profiles for heptane become progressively flatter in the flame 
region. The same thing happens for the oxidiser mass-fraction profiles also. It is clearly 
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understood from Fig. 5.5a-c that the fuels and the oxidiser are almost fully consumed at 
the flame zone. The figures certify that the flame is not thick. Figure 5.5b depicts mainly 
the intermediate period including first <P - law period, transition and second - law 
period, which is extended upto extinction ( refer Fig. 5.4c ). In the extinction period, it 
is observed that as the flame moves towards the end of the solution domain ( as a result 
of reduced droplet diameter ), the oxidiser mass fraction profiles at the outer boundary 
becomes steeper due to the enforcement of the boundary conditions. This happens up to 
0.0655396 s. Then again the oxidiser mass fraction profile becomes more and more flat at 
the flame zone and the flame moves inward towards the droplet surface at a rate faster than 
that of the surface regression rate. Therefore, the flame again moves inward relative to the 
outer boundary of the solution domain. The heptane concentration gradient at the droplet 
surface falls steadily because of the volatility limited mode of evaporation. From figure 
5.5c, it is quite evident that the consumption at the flame surface becomes progressively 
lesser and the oxidiser mass fraction profile gradually tends towards it’s initial uniform 
distribution which was prevailing before ignition. This feature is more clearly shown in 
Figure 5.6c. 

Figures 5.6a-c show the mass-fraction profile for the less volatile fuel ( hexadecane ). 
The trend of the profile histories for both fuels is of similax nature. The only difference is 
that the surface mass-fraction for the more volatile fuel grows much faster and is maintained 
at a higher value ( 0.84 ) due to it’s higher initial liquid side value ( 0.9 ) and this growth 
is over in the initial period only, whereas for the less volatile fuel, this change takes place 
more slowly over a long period and is maintained at a much lower value ( 0.16 ) owing to 
it’s lower initial liquid side maas-fraction ( 0.1 ). Figure 5.6c shows that towards the end of 
burning, the mass fraction of hexadecane rises very sharply both at the droplet surface and 
also away from the droplet surface. Increase at the droplet surface is due to the volatility 
limited mode of evaporation and this increase is also penetrated away from the droplet 
surface because the consumption at the flame becomes negligibly small. 
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Fig. 5.4 Evolutiou of gas side temperature profiles and transient variation of flame position 
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5.5 Effect of Initial Mass-fraction Variation on Sur- 
face Parameters 

Surface parameters namely T„ V, Hi, Yj.i,, Yj,,„ r,, <P, r^, 

Tf and were discussed in section 5.2. This section heiglights the effect of varying the 
initial mass- fraction combination on these parameters plus a few other related parameters. 
Figures 5.7 through 5.9 show the variation with time of Z, V, Hi, Yfn,, Yyf^^g, e/i, Yj 2 sg, 
^ 1 '^'fi various values of initial mass-fraction combination. 

For 0.9 initial liquid side mass fraction of heptane, the entire burning period is divided 
into six zones ( see Fig. 5.7a ) as discussed in Section 5.2 ( refer Fig 5.1a-c ). The two 
(f — law periods are clearly visible from Fig. 5.9a. The relative length and positions of 
the two dP — law periods and the transition zone in between depend on the initial mass- 
fraction of the fuels. It is seen from Fig. 5.9a that when the initial mass-fraction of the 
more volatile fuel ( heptane ) is very high ( 0.99 ), the 2"*^ <P — law period is almost 
absent and only the — law period is present all through before the extinction occurs. 
This happens because, for these cases, the droplet is always lacking in less volatile fuel ( 
i.e. dodecane ), specially at the surface ( Fig. 5.8a ). The transition zone in this case is 
shifted to the final stage of burning- But as the initial liquid side mass fraction of heptane is 
gradually decreased, influence of dodecane over the evaporation rate becomes progressively 
more significant. This is seen from Fig. 5.8c ( vs. t ). It indicates that as the initial 
mass-fraction of heptane is gradually decreased, domination of heptane evaporation is also 
gradually decreased and the transition becomes faster and is over in the early stage of 
burning. Thus, length of first <f-law period gradually decreases and length of the second 
dP ~ law period gradually increases. When the initial mass-fraction of heptane becomes 
very small (e.g. 0.1), the first (P - law period almost vanishes and the 2"“^ dP - law period 
covers almost the entire time length except for the initial transient period because, we can 
see from Fig. 5.8d that for this case, the surface is almost totally dominated by dodecane. 

The transition from the more volatile fuel burning to the less volatile fuel burning 
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becomes most distinct when the initial mass-fraction of heptane is 0.9 ( see Figs. 5.7a, 5.8 ' 
b d and o.9a ). As the mass-fraction of heptane is lowered, the transition zone becomes 
steeper and it shifts more towards the early part of droplet burning and vice versa. Figures 
5.7a ,5.8a-d and 5.9a clearly show this change of burning characteristics as the initial mass- 
fraction combination of the two fuels gradually changes from one extreme to the other. 

Figs. 5.7a and 5.7c indicate that droplet heating is more rigorous and it takes place 
in the early stage of droplet burning when initial mass-fraction of heptane is small (0.1). 
Since, the amount of heptane present inside the droplet is small, the surface is dominated 
by dodecane almost from the very beginning. Therefore, the surface temperature is also 
controlled by the boiling point of dodecane (472 K) from the very beginning. So, for this 
case, only one long initial transient exists. But, as the initial mass-fraction of heptane 
is increased gradually, the droplet surface becomes dominated with heptane for longer 
periods. So, the second transient gradually shifts to the later part of the droplet burning 
and essentially, the initial droplet heating as well as the rise of surface temperature take 
place in two stages depending on the initial value of liquid side heptane mass-fraction. 
This effect of initial liquid side mass-fraction on droplet heating and surface temperature 
is revealed in the Figs. 5.7a and 5.7c. 

Another interesting feature of the Figs. 5.7a through 5.9a is^that as the initial mass- 
fraction of heptane is gradually increased, the total burning time is gradually decreased and 
is least for pure heptane fuel burning. This happens because the droplet heating associated 
with the second transition becomes lesser due to the larger domination of heptane fuel ( 
see the area under the curves of Fig. 5.7c ) and also, more heat is available for fuel 
evaporation during the same duration. This is corraborated by Fig. 5.7d. The area under 
each volumetric mass gasification rate, V curve represents the relative amount of total 
heat available for evaporation and one can see from Figs. 5.7b,d that more heat becomes 
available for evaporation when the initial mass-fraction of heptane is higher within the 


same time limit. 



Heat going inside the droplet, H, ^ Surface Temperature (IC), 


: o 

1 .7 


Time in second 


.075 0 

0 

hep-dod 


Time in second 


.03 .045 

Tiiiie ia itcoad 

C 


.075 0 


.03 .045 

Time in second 


Fig. 5.7 Effect of initial mass-fraction on the transient variation of surface temperature, 
vol. mass gasification rate, droplet heating and heat available for evaporation. 
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Figures 5.8a-c show that as the initial mass-fraction of heptane is gradually increased, 
the curves become flatter and they are maintained at progressively higher values because, 
more amount of heptane becomes available at the droplet surface. Figure 5.8d shows that 
gas side surface mass-fraction profiles for dodecane are complementary to those of heptane 

mass fraction profiles. Fig. 5.9b indicates that the flame moves almost upto the the same 
distance for all the cases. 

5.6 Effect of Volatility-Differential on Burning Char- 
acteristics 

The nature of multicomponent combustion is highly dependent on the volatility-differential 
of the two individual fuels as can be seen from Figs. 5.10 through 5.12. Figure 5.12c shows 
that when the volatility-differential is very low ( e.g. for heptane-octane combination ), 
the two cP — law periods merge with each other and become one single (P — law period. 
As the volatility-differential between the two fuels gradually increases from heptane-octane I 
combination to heptane-dodecane , the two — law periods become more distinct and 
the transition from the (P — law period to the 2”'^ (P — law period becomes more : 
developed. The intermediate transition and thereby the multicomponent feature of fuel 
burning becomes more prominent when the volatility-differential is increased ( see Figs. ; 
5.10a, c; 5.11a-d; 5.12c ). 

One can see from Figs. 5.10a, 5.10c and 5.12c that When the volatility-differential is j 
low, most of the droplet heating is over during the initial transient and then it is followed | 
by a long quasi-steady period covering most of the droplet lifetime. As the volatility 
differential between the two fuels decreases, it approaches single component fuel burning i 
in the limiting case. When the volatility-differential becomes large, it is seen that the 
droplet heating takes place twice. Firstly, just after ignition takes place and secondly, when 
the surface temperature msices a move beyond the boiling point of the more volatile fuel. 
This second transient, basically, corresponds to the transition zone where the combustion : 
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process moves from - law period to the - law period. 

Figure 5. lOb indicates that the overall evaporation rate has a longer quasi-steady period 
(similar to what is seen for a single component fuel) when the volatility-differential is small. 
As the volatility-differential is gradually increased, two quasi-steady periods separate out 
corresponding to the two d? - law periods. 

B igures 5.11a-d show how the liquid side and gas side surface mass-fractions of the two 
fuels as well as their relative strength of evaporation depend upon their volatility differen- 
tial. Figure 5.11a shows that the liquid side surface mass-fraction of the more volatile fuel 
(heptane ) falls more steeply as the volatility differential between them increases. Figure 
5- lie indicates that the relative strength of evaporation for the two fuels remains almost 
constant for low volatility differential, whereas, for high volatility differential, two distinct 
regions separate out corresponding to the two d? — law periods. Similar trend is also ob- 
served for Yfiag a-nd Yf 2 sg curves in Figs. 5.11b and 5. lid. The gas-side mass-fractions for 
heptane ( Fig. 5.11b ) initially increase and then decrease continuously, whereas, for the 
less volatile fuel ( Fig 5.1 Id ), they increase continuously. This happens because, when the 
volatility differential is very low ( e.g. heptane-octane), both the fuels evaporate with equal 
preference. Therefore, almost constant values are maintained for all the curves for the low 
volatility differential case. Whereas, when the volatility differential is progressively higher, 
more volatile fuel is evaporated from the surface with higher preference. Therefore, the 
liquid side surface mass-fraction of the more volatile fuel falls more steeply. As a result of 
which relative evaporation strength, c/i, as well as it’s gas side surface mass-fraction also 
fall Consequently, as the surface temperature makes a move towards the boiling point of 
dodecane, the gas side surface mass-fraction for dodecane increases. During the 2'^^d?law 
period, Rractional mass-gasification rate and, liquid and gas-side mass-fractions of both 
fuels are maintained almost at constant values. As the volatility-differential is decreased, 
the distinction between the two periods vanishes into almost one long period. 
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Fig. 5.10 Effect of volatility-differential on the transient variation of surface temperature, 
vol. mass gasification rate and droplet heating ’ 
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Fig. 5.11 Effect of volatility-differential on the transient variation of liquid and gas phase 
surface mass-fractions of heptane, gas phase surface mass-fraction of do decane and frac- 
tional mass gasification rate of heptane. 
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Fig. 5.12 Effect of volatility-differential on the transient variation of flame position, flame 
stand-off-ratio and square of the diameter. 







5.7 Comparison of Analytical Models 

This section describes the characteristics of the following three models: 

(i) Quasi-steady model 

(ii) Model neglecting the effect of surface regression 

(iii) Present model. 

Figures 5.13 through 5.15 show the variation with time of T/i/s, Yfisgi F/ 2 S 35 Hi, Tg 
rj, r'j and V for two values of initial mass-fraction combination for each of the above 
mentioned models. In the full quasi-steady model, the surface mass-fractions of both fuels 
are assumed to be maintained at the initial mass-fraction values throughout the entire 
burning history. For the present study, the full quasi-steady model (i) is valid as a special 
case when surface regression is assumed to be very high from the very beginning. Effect of 
maas diffusion is very small for this model. The continuous surface regression gives rise to a 
moving boundary problem. When the surface is fixed at r=l, the transformed coordinate 
becomes moving in nature. This transformation gives rise to a convection term in the 
transport equations for the liquid and phases as discussed earlier. In the second model, 
this convection term is neglected, considering that the regression rate of the droplet radius 
is very small This case may be applicable in case of pure vaporization where the rate of 
evaporation is very low. But, in practical combustion, situations one can neither neglect 
the liquid side mass-diffusion ( specially towards the end of burning ) nor surface regression 
rate, because both are significant. Therefore, the present model ( i.e. the third model), 
includes both mass-diffusion and surface regression. From the results also, one can see that 
the practical situations are combinations of these two extreme limits. 

For elaborating upon the comparison, two cases of 0.9-0.1 and 0.5-0.5 initial mass- 
fraction combinations of heptane-dodecane mixture are chosen. From Fig. 5.13a, one can 
see that for the second model in which the effect of surface regression is neglected, the 
mass-fraction of the more volatile fuel at the surface falls very steeply. This happens due to 
the high liquid side diffusional resistance. One can see that the present model resembles 
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the second model only in the initial stage of droplet burning when the surface regression 
rate is very low due to the low evaporation rate. As evident from the figure, the degree of 
this tendency depends on the initial mass-fraction of the fuels. They closely approach each 
other if the initial mass-fraction of heptane is low. As the initial mass-fraction of heptane is 
increased, the present model increasingly deviates from the second model. As the burning 
proceeds beyond this initial period, surface regression rate gathers momentum due to more 
evaporation at the droplet surface. Therefore, the inner core is exposed to the surface very 
rapidly, supplying the necessary amount of fuels for maintaining the surface mass-fraction. 
So, for this later part of the burning, the present model behaves like the quasi-steady model 
in which the surface mass-fraction of the fuels remains almost constant. This behaviour 
continues till the end of the burning. The similarity of the first and third models in the 
later part of droplet burning is only qualitative. As evident from the Figs. 5.13 through 
5.14, the qua.si-steady model (1) basically behaves like a single component fuel burning 
model in the sense that it has only one oP — law period ( see Fig. 5.14d ). But actually, the 
burning is controlled by the surface mass-fractions of the two fuels which, in turn, depend 
on both liquid phase mass- diffusion and surface regression rates. Therefore, when surface 
regression is not neglected, the role of initial mass-fraction, as evident from Fig. 5.13a, 
on burning characteristics becomes more important. For the quasi-steady model, the total 
droplet heating is much less which is represented by the area under the curves of Fig. 5.14b. 
Consetiuently, for this case, more heat becomes available ( as a percentage of total heat 
available ) for evaporation. It is seen from Fig. 5.14d that when the surface regression 
is neglected ( model 2 ), the transition from one (P - law period to the next becomes 
very sharp and the duration of the 1** <P - law period as well as the total burning time 
are underpredicted. Compared to quasi-steady model, for models 2 and 3, a progressively 
larger amount of heat goes into the droplet and consequently, lesser amount of heat is 
available for evaporation- Hence, for all practical purposes, the regression rate cannot be 
neglcct^id. Figures 5.15a,b show fiame movement for theses three models. It is seen that 
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when the convective terms for the liquid and gas phase equations ( arising out of surface 
regression ) are neglected ( in model 2 ), the flame comes back towards the droplet surface 
very gradually as against the movement of the flame for other two cases where the flame 
remains steady during the later half of the burning. And it comes back towards the droplet 
surface very steeply at the end period of the burning. This difference of flame movement 
occurs due to the neglecting of the convective term for the gas phase equations for model 
2 because, the gas phase does not receive any information regarding the surface regression. 

5.8 Gas Phase Dissociation 

Various combustion parameters are strongly affected by gas phase dissociation. Figures 
5.16 through 5.18 show the variation with time of T/i/s, Yfisg-, ^fi, Yf 2 ag, Hi, rj, r//, H^v, 
T,, d^, Tf and V with and without gas phase dissociation for heptane-dodecane burning 
with 0.9-0. 1 initial mass-fraction. 

Figures 5.16a and 5.16b show the effect of dissociation reaction on flame temperature, 
flame position and fuel vapour acummulation. When gas phase dissociation is included, 
flame temperature is considerably reduced due to the endothermic nature of dissociation 
reactions. The flame, however, moves in similar manner for both the cases ( i.e. with and 
without dissociation as seen from Fig. 5.16b ), the only difference being that the flame 
without dissociation, because of higher flame temperature which leads to slightly higher 
rate of fuel vapour acummulation, moves faster and also farther away from the droplet 
surface, and finally, comes back faster and vanishes into the droplet. Higher value of heat 
generation at the flame causes larger rate of heat flow from gas to liquid phase ( Fig. 5.17a 
}. This results in a higher rate of surface temperature rise ( Fig. 5.17b). Larger rate of 
heat is also available for evaporation at the droplet surface ( Fig. 5.17a ). Consequently, 
gasification rate is also increased ( Fig. 5.17c ) which results in faster rate of fuel vapour 
.ccummulation md Same movement ( Fig. 5.16a ) and also faster regression of the droplet 
surface (Fig. 5.17d ). Liquid side surface mass-fraction of heptane falls more steeply ( 
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Fig. 5.16 Effect of gas phase dissociation on the transient variation of flame temperature 

flame position and flame stand-off-ratio . 
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Fig. 5.17 Effect of gas phase dissociation on the transient variation of droplet heating, 
heat available for evaporation, surface temperature, vol. mass gasification rate and square 


of the diameter. 
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Fig. 5.18 Effect of gas phase dissociation on 
phase surface mass-fractions of heptane, fracti 
gas phase surface mass-fraction of dodecane. 







Fig. 5.18a ) due to higher initial gasification rate and reaches quasi-steady state faster 
due to relatively high regression rate in the earlier part. The temporal variation of gas- 
side mass-fraction of heptane and dodecane, and also that of fractional mass-gasification 
rate of heptane ( Figs. 5.18b,c,d ) behave in a similar manner. Liquid side temperature 
profiles do not indicate any significant change. Any change in the gas-phase is scaled down 
considerably in the liquid phase due to higher inertia of the liquid phase. 

5.9 Effect of Ambient Temperature 

In this section, Figs. 5.19 through 5.20 show the variation of different surface parameters 
such as T/, Hi, T,, V, YfUs, Yfug, ^fi and Yf 2 sg have been plotted against time for three 
different ambient temperatures namely 300 K, 500 K and 1000 K for heptane-dodecane 
burning with 0.9-0. 1 initial mass-fraction. When the ambient temperature is 1000 K, 
the flame temperature is maintained around 3000 K and for 500 K and 300 K, the flame 
temperatures are maintained around 2760 K and 2665 K respectively ( see Fig. 5. 19a ). This 
happens because the reaction is more vigorous for higher ambient temperature. Therefore, 
more heat becomes available for both droplet heating ( see Fig. 5.19b ) and for fuel 
evaporation. Consequently, the surface temperature, T,, increases at a faster rate and 
it reaches a higher value ( Fig. 5.19c ). The volumetric mass gasification rate, V, also 
becomes faster ( Fig. 5.l9d ) due to the easy availability of heat for fuel evaporation, H^ 
and therefore, the burning is completed earlier ( Fig. 5.19e ). 

Figures 5.19a-e and 5.20a-d reveal that the transition from l*‘d^ - law period to the 
— law period takes place earlier and the curves are shifted slightly to the left for 
the higher ambient temperature. Figures 5.20a and 5.20b indicate that for the higher 
ambient temperature, the liquid side mass-fraction for heptane falls more steeply before it 
reaches a constant value and it’s gas-side mass-fraction reaches higher value faster and this 
higher value is maintained till the end of burning. As explained already, the volumetric 
mass gasification rate is higher for higher ambient temperature. Therefore, for this case, 
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Fig. 5.19 Effect of ambient temperature on the transient variation of flame temperature, 
droplet heating, surface temperature, vol. mass gasification rate and square of the diame- 
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heptane evaporates at a faster rate and as a result, it’s gas side surface mass-fraction 
increases at a faster rate and is maintained at a higher value. Figures 5.20c and 5.20d 
show that the variation in ambient temperature does not influence the fractional mass 
gasification rate for heptane and the gas side surface mass-fraction for dodecane during 
the initial and final part of burning. The influence is mainly in the transition region. 

5.10 Variable Specific Heat for the Gas Phase 

In the present section, the surface parameters, viz. Tj Hg, Hi, Hgv, Ts, V, r^, ry/, YfUs 
and (P are chosen to highlight the effect of variable specific gas phase heat as against 
three constant values of specific heat for heptane-do decane burning with 0.9-0. 1 initial 
liquid side mass-fraction. The values chosen for the three constant specific heat cases are 
1.0xl0^ 1.25x10^ and 1.5x10^ J/kg-K. 

It is observed that when variable specific heat is used, the flame temperature comes 
down to around 2950 K as compared to the constant specific heat ( with value 1.0x10^ 
J/kg-K )case when the flame temperature remains around 3000 K ( see Fig. 5.21a ). But, 
for all other surface parameters, the variation for the two cases ( curves 1 and 2 ) are very 
small. The high specific heat at the flame region ( for variable Cp case ) does not allow the 
flame temperature to exceed a certain limit and that’s why the flame temperature remains 
lower for this case. But, the temperature near the droplet surface and other places except 
for the flame region are not very high. Therefore, near the droplet surface and other places 
excepting the flame region, the value of specific heat for the above two cases are very close 
to each other. Therefore, the temperature gradient at the surface and consequently, the 
available heat from the gas phase for the two cases remains close to each other ( 5.21b | 
). Therefore, almost the same amount of heat are available for the two cases for droplet 
heating ( Fig. 5.21c ) as well as for fuel evaporation ( Fig. 5.21d ). So, all other parameters | 
such as surface temperature (Fig. 5.22a ), volumetric mass gasification rate ( Fig. 5.22b ), : 
flame movement and flame stand-off ratio ( Fig. 5.22c ), liquid side surface mass fraction / 
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Fig. 5-21 Effect of variable specific licat a,s well as tliffcrciil values of constant specific heat 
on the transient variation of flame temperature, lieat coming from the gas phase, droplet 
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Fig. 5.22 EfTect of variable specific heat as well as different values of constant specific 
heat on the transient variation of surface tenii)eraturc, vol. mass gasification rate, flame 

position, flame stand-off- ratio, liquid |)ha.sc surface mas.s-fractions of heptane, and square 
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Fig. 5.22d ) vary negligibly for these two cases. For the constant specific heat cases, we find 
that as the gas side specific heatis gradually increased, the flame temperature gradually 
decreases ( Fig. 5.21a ) and since specific heat is increased all over the domain uniformly, 
the temperature decreases over the entire domain. Therefore, the heat available from the 
gas phase becomes gradually less ( Fig. 5.21b ) due to the lower temperature gradient at 
the gas side of the droplet surface. Consequently, heat available for increasing the surface 
temperature ( Fig. 5.21c ) and fuel evaporation ( Fig. 5.21d ) also becomes less. Therefore, 
overall volumetric mass gasification rate falls ( Fig. 5.22b ) and consequently, evaporation 
rate constants become gradually lower ( Fig. 5.22e ) and the liquid side mass fraction 
of heptane falls with lower steepness ( Fig. 5.22d ). The nature of the flame movement 
remains similar ( Fig. 5.22c ). As a result of lower rate of fall of liquid side heptane mass- 
fraction at the droplet surface, surface temperature ( Fig. 5.26a ) increases more gradually 
over a larger time span. Therefore, the total burning time increases ( Fig. 5.22e ). 

5.11 Gas Side Diffusion Coefficient 

Figures 5.23 shows the the effect of varying the gas phase diffusion coefficient on the over- 
all burning characteristics for heptane-dodecane burning with 0.9-0. 1 initial mass-fraction 
combination. Diffusion coefficient of value 2.248 x 10~® m^/s corresponds to unity gets 
phase Lewis number. Values of gas phase diffusion coefficients used for the four curves in 
these figures are 1.0 X 10"®, 2.248 X 10"®, 3.5 X 10"® and 5.0 x 10"® mVs- It is observed 
that as the diffusion coefficient is gradually decreased, the total burning time gradually 
increases and the intermediate transition becomes more gradual ( Fig. 5.23a-e ). As the 
diffusion coefficient of the species in the gas phase becomes progressively lower, the dif- 
fusion of the evaporated fuel vapour away from the droplet surface becomes progressively 
lesser. Therefore, effectively, the vaporization of the more volatile component ( heptane ) 
gets relatively suppressed. So, the liquid side surface mass-fraction of heptane ( Fig. 5.23b 
) falls less steeply zmd since the surface temperature depends mainly on the surface mass- 
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fraction of the two fuels, the surface temperature ( Fig. 5.23c ) rises with lesser steepness. 
For the lower value of the gas phase diffusion coefficient, more heat is available for droplet 
heating ( Fig. 5.23d ), and therefore less heat is available for the fuel evaporation ( Fig. 
5.23e ). Consequently, it takes longer time for combustion due to lower evaporation rate 
constant ( 5.23a ). 

5.12 Liquid Side Diffusion Coefficient 

Figures 5.24 through 5.26 show the effect of liquid side diffusional resistance on the overall 
gasification behaviour for heptane-hexadecane fuel with 0.7-0. 3 initial mass-fraction com- 
bination. Here the curves are plotted for different liquid side diffusion coefficients having 
values 5.0 x 10“^°, 5.0 x 10~°^, 5.0 x 10“°® and 5.0 x 10"°® m^/s. When the diffusional 
resistance is low, the first cP — law period ( controlled by the burning of the more volatile 
fuel ) exsists for a longer period due to it’s supply from the interior of the droplet ( Fig. 
5.24a ). But as the diffusional resistance becomes higher, the length of the first <P — law 
period decreases because the mass-fraction of the more volatile fuel decreases at the surface 
due to it’s insufficient supply from the droplet core. Therefore, for this case, the second 
(P — law period exists for major part of the burning. Here, one interesting point to be 
noted is that for the case when the liquid side diffusion coefficient is very high ( greater 
than 5.0 x 10"® xrP/s ), as shown for curve 4 of the Figs. 5.24 through 5.25, there exists : 
two distinct batch-distillation mode of burning regimes corresponding to the two different I 
fuels. These two regimes are separated by one extremely narrow zone when very high rate | 
of droplet heating takes place ( Fig. 5.24d and 5.25a ). This is also very distinctly observed? 
in all the Figs. 5.24 through 5.25 for curve 4. Due to the very high value of liquid side? 
diffusion coefficient the more volatile fuel ( heptane ) is constantly supplied to the surface: 
as evident from Fig. 5.24b. So long as it is available inside the droplet interior ( Fig.; 
5.26 ). Therefore, batch-distillation mode of burning for heptane occurs during this first 
^ ^law period. This is also very clearly shown in Fig. 5.24c where the fractional mass 
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Fig. 5.24 Effect of liquid side mass diffusion coefficient on the transient variation of square 
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Fig- 5.25 Effect of liquid side mass diffusion coefficient on the transient variation of droplet 
heating, heat available for evaporation, vol. mass gasification rate, flame position, and 
flame stand-off- ratio. 










gassification rate of heptane is almost 1. During this period, surface temperature remains 
near the boiling point of heptane. Then when the droplet interior is almost completely 
devoid of heptane fuel ( This takes place at around t=:0.06 s for I) = 5.0xl0“^mVs ( curve 
3 for Figs. 5.24 and 5.25 ) and at around t=0.035 s for D = 5.0xl0““^mVs ( curve 4 for 
Figs. 5.24 and 5.25 ) ), the liquid side surface mass-fraction of heptane falls very sharply 
( see also Fig. 5.26 ) and becomes almost zero. Therefore, the surface temperature, being 
controlled by the surface fuel mass fraction combination, rises very steeply and reaches 
very near to the boiling point of hexadecane ( 560 K ) in a very short time ( Fig. 5.24d ). 
This results in vigorous droplet heating ( Fig. 5.25a ) and consequently, heat available for 
evaporation falls steeply ( Fig. 5.25b ) during this short time span. Therefore, volumetric 
mass gasification rate also falls very steeply ( Fig. 5.25c ) for which sudden shrinkage of 
flame takes ( prominent for D = 5.0xl0~®m^/s, curve 4 ) place as indicated in Fig. 5.25d. 
This is indicated by the sudden lowering in slope of curve 4 in Fig. 5.24a. It is to noted 
that if thermal expansion of the liquid fuel were taken into account, then this change in 
slope for those two curves of Fig. 5.24a would be even sharper. After this transition, I 
second (P — law period corresponding to the second batch-distillation mode of hexadecane I 
burning takes place ( Fig. 5,24a ) and it continues till the end of burning. The domination | 
of hexadecane is clearly visible in Fig. 5.24c where the fractional mass-gasification rate of i 
heptane is shown to be practically zero during this period. From this paxametric study,: 
one can see that for a very high value of liquid side mass-diffusion co-efficient, two batch 
distifiation modes of burning can exist for a multicomponent fuel. 

5.13 Comparison with the Semi Analytical Model 

Figure 5.27 shows a comparison of the present numerical model with the semi-analytical 
model of chapter 3. The results are shown for heptane-undecane with 0.5-0. 5 initial mass- 
fraction combination. Curve 1 represents the full numerical model where transient feature 
of the gas phase as well as the effect of surface regression are retained. Curve 2 represents 
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Fig. 5.27 Comparison of the semi-analytical model with two czises of the numerical moc 
where for one case, the effect of surface regression is neglected and for the other case, t 
effect of surface regression is considered. 
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the numerical model where the effect of surface regression is neglected. Curve 3 represents 
the semi-analytical model where the effect of surface regression is neglected and also, the 
gas phase is assumed to be quasi-steady. Curves 1 and 2 matche perfectly in the initial 
stage of burning, but deviate during the later part. This happens because, initially, the 
effect of regression of droplet diameter is small and can be neglected. But, later on, the 
droplet diameter regression rate becomes quite high and therefore, it can not be neglected. 
The curve 3 deviates from curves 1 and 2 in the initial stage because, the semi-analytical 
model doesnot include the transient feature of the gas phase. Later on, it deviates more 
from curve 1 but, matches with curve 2 because the latter ( i.e. curve 2 ) neglects the effect 
of surface regression rate. Therefore, the transient feature of the gas phase is prominent in 
the initial stage of droplet burning and consequently, the assumption of gas phase quasi- 
steadiness can well be applied for the later part of the burning. 

5.14 Comparison with Earlier Work 

Different values of evaporation rate constants for heptane burning at room temperature 
are available in literature ( Aldred and Williams, 1966; Graves, 1952; Smith and Graves, 
1957; Monaghan et ah, 1968; Goldsmith, 1956 ). They are reported to vary from 0.8428 
mm^/s to 1.0 mm^/s according to Fig. 10 of Williams ( 1973 ). Kumagai and Isoda ( 1957 
), in their experimental results, reported the value of the evaporation rate constant to be 
0.49 mm^/s for zero gravity case for suspended droplet burning. Theoretical value of this 
constant calculated by Goldsmith and Fenner ( 1954 ) was reported to be 0.86 mm^/s for 
steady state combustion. Kumagai and coworkers ( 1956, 1971 ) found the evaporation rate 
constants ( for a 1 mm initial diameter ) as 0.75 mmVs for a suspended droplet and 0.79 
mm^ /s for a free droplet of heptane. Figure 5.28 shows comparison of the present numerical 
model with the experimental results of Okazima and Kumagai ( 1975 ) and Wang et al. ( 
1984 ) under two different conditions for heptane burning. Okazima and Kumagai ( 1975 
) presented results ( curve 1 of Fig. 5.28 ) for spherical combustion under zero gravity 
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condition in nornaal ambient conditions. According to Okazima and Kumagai ( 1975 ), 
evaporation rate constant is 0.78 mm^/s. The present model predicts the evaporation rate 
constant to be 0.70175 mm^/s ( see curve 2 ) and 0.6837 mm^/s ( see curve 3 ) when liquid 
side specific heats are 1.2 x 10^ and 1.75 x 10^ J/kg-K respectively. It is seen that the 
liquid side specific heat generally varies from 1.2 to 2.0 J/gm-K ( Reid et al., 1988 ) The 
variations are 9.1 % and 12.3 % respectively. The initial droplet diameter for Okazima and 
Kumagai was 1.327 mm whereas, for the present problem, the initial diameter is 0.2 mm. 

Experimental results of Wang et al. ( 1984 ) are shown as curve 4 for heptane droplet 
moving at a low Reynolds number ( < 1 ) where the ambient temperature is as high as 1200 
K. Curves 5 and 6 show the result obtained from the present model for 1200 K ambient 
temperature. The experimental curve 4 predicts the evaporation rate constant to be 0.9836 
mm^/s, whereas the present model estimates it to be 0.7565 mm^/s ( curve 5 ) and 0.73059 
mm^/s ( curve 6 ) for 1.2 x 10^ and 1.75 x 10^ J/kg/K liquid side specific heat respectively. 
The variations are 23 % and 25.7 % respectively. It is seen that the curves 4, 5 and 6 
show good agreement in the initial stage of the burning, but deviate later. Experiments 
of Wang et al. ( 1984 ) show that the droplet velocity was small in the initial stage 
and gradually increased towards the the end of the experiment under natural convection 
condition. Therefore, curves 3 and 4 show better agreement in the initial stage of burning. 
Okajima and Kumagai ( 1975 ) reported that the evaporation rate constant increases as 
the external flow ( due to forced convection ) was increased and this increase was 10 % at 
Uoo==2 cm/s when the ambient temperature was around 300 K and then remained almost 
constant. Their initial droplet diameter was 1.4 4^0.2 mm. Kumagai and Isoda ( 1957 ) in 
their suspended drop experiment reported that the evaporation rate constant varieed from 
0.49 mm^/s to 1.02 mm^/s for zero gravity and full gravity conditions respectively. Thus, 
they indicated the effect of natured convection to be quite significant for the evaporation 
rate constant. The droplet velocity reported in the experiment of Wang et al. ( 1984 
} under natural convection condition, was 90 cm/s on an average and they used droplet 



diameter in the range 25 to 35 ixm. Therefore, the deviation of curves 5 and 6 from curve 
4 in the later part of the burning could be attributed to the external flow as well as greater 
initial droplet diameter. 

Apart from the above mentioned reasons, other property variations specially in the 
gas phase ( where the temperature varies over a wide range ) and neglecting radiative 
heat-transfer can also be attributed to the cause for the overall deviation of curves 2 and 
3 from curve 1, and also of curves 5 and 6 from curve 4. In addition, the reaction kinetics 
involved at the flame region is extremely complex. While modelling, these different complex 
situations are simplified to a great extent. Also, it is extremely difiicult to estimate the 
accurate property values. As shown earlier, the dependence of burning characteristics on 
property values is sensitive. Keeping all these assumptions in account, curves 2 and 3 show 
a reasonably good agreement with the experimental curve 1. 

Fig. 5.29a shows the comparison of the flame stand-off ratio ( ratio of flame radius to 
the droplet radius ) vs. t/ro curve for two cases. Curve 1 shows the experimental results : 
of Okazima and Kumagai ( 1975 ) and curve 2 represents the present model. They agree i 
very well upto 0.75 s/mm^ and then deviate from each other. Similar agreement is also ; 
observed for Figs. 5.29b, c,d. For the above two cases. Fig. 5.29b shows the comparison : 
for the variation of normalised flame radius ( flame radius divided by the initial droplet : 
radius ) vs. t/r^ and Fig. 5.29c shows the comparison for the normalised flame radius ! 
versus normalised droplet radius. Fig. 5.29a show the flame movement in relation to the 
movement of the droplet radius and Fig. 5.29b show the flanne movement with respect to | 
the droplet center. In Fig. 5.29c, curve 1 shows the experimental results of Okazima and ; 
Kumagai ( 1975 ) and curve 2 is plotted for the present model. This figure shows that j 
the flame radius is proportional to droplet radius, only towards the end of droplet burning 
when extinction occurs. Therefore, the steady state burning reported by Goldsmith and j 

Fenner ( 1954 ) is not supported here. j 

Figure 5.30a shows a comparison of the present model with the experimental results j 
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of by Wang et al. ( 1984 ) for the burning of multicomponent ( mixture of heptane and 
hexadecane with an initial mass-fraction combination of 0. 7-0.3 ) fuel droplet. The result 
shows an excellent agreement in the initial part of the droplet burning. But the two 
results deviate as burning proceeds in time. The present numerical model predicts the 
evaporation constant to be 0.65 mm^/s for the first (P — law period and 0.7 mm^/s for the 
second <P — law period ( curve 2 ). Whereas the experimental results ( curve 1 ) predict 
the evaporation rate constant to be 1.022 rrmp/s for the first (P — law period, and 0.971 
mm^/s for the second cP — law period. For the curve 2, Cp = 1.75 x 10^ J/kg-K, A = 0.1 
W/m-K and .D = 5.0 x 10“^ rxP/s were used. Curve 2 shifts towards curve 1 when the liquid 
side specific heat is decreased. One of the causes attributed to the variation in the total 
burning time can be the natural convection ( average droplet velocity = 90 cm/s ) present 
in case of Wang et al. ( 1984 ) and this effect of natural convection is not considered in the 
present model. As already discussed, according to Kumagai and Isoda ( 1957 ), this effect 
of natural convection on evaporation rate constant is very strong. For the experimental 
results, the droplet burns out much faster owing to larger value of the evaporation rate 
constants. The experimental result shows that the droplet diameter does not change in 
between 0.053714s and 0.059428s and remains at 0.29179mm. Wang et al. ( 1984 ) have 
claimed that this sharp change when burning shifts from the first <P — law period to the 
second is caused mainly due to lowering of fuel evaporation rate as well as the thermal 
expansion of the liquid droplet which becomes prominent due to the intermediate stage of 
vigorous droplet heating during transition. Thus, their experiment showed three periods, 
the two periods correspond to the preferential gasification of the two fuels and the third 
period corresponds to the short intermediate transition in between as mentioned above. 
The present model also exhibits three dominant zones, two <P — law periods separated by 
one transition zone. But, this transition from one (P — law period to the next is gradual 
and takes place over a longer period of time for the present model. 

After the initial transient is over, the temperature at the droplet surface will tend 
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Fig. 5.30 Comparison with the exprimental results of Wang et al. ( 1984 ) for heptane 

hexadecane fuel burning with O. 7 .O .3 initial mass-fraction combination. 
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to approach closely the equilibrium boiling temperature of the mixture at the droplet 
surface. The boiling point of the mixture , in turn, depends upon the relative strength 
of the constituents at the surface. Now, since the liquid side mass diffusional resistance 
is high and initially, the surface regression is low, the liquid side mass fraction of the 
more volatile fuel falls sharply. Consequently, the gas side surface mass-fraction of the 
more volatile fuel also starts falling quite early in the droplet lifetime. But because of the 
continuity of the burning process and also due to the fact that as the burning proceeds, 
more amount of the more volatile fuel is supplied to the surface as a result of droplet 
radius regression, the change in the relative strength of the constituent fuels at the droplet 
surface takes place gradually. Therefore, the boiling point at the droplet surface, being 
dependent on the relative strength of the constituent fuels, also starts rising gradually. So, 
the droplet heating becomes gradual during this transition. It may be possible that the 
thermal expansion of the liquid droplet is quite high during the transition period and since 
the present model does not take into account the effect of thermal expansion, this may be 
responsible for the deviation in the predictions of the present model from the experimental 
observations of Wang et al. ( 1984 ). 

However, it is worth mentioning here that this sudden droplet heating is very prominent 
for the case where, two distinct batch-distillation mode of burning take place. This is 
possible when the internal transport within the droplet is very high. This may be possible 
for the case when the external velocity is quite high ( i.e. for high Reynold’s number flow 
) which may give rise to a very high value of internal circulation ( Law, 1976 ). Prakash 
and Sirignano ( 1978 ) have shown that spatial uniformity is not achievable even with an 
infinite circulation rate inside the droplet. 

It is seen from Fig. 5.7a that when the volatility differential is sufficiently large ( e.g. 
for heptane-dodecane mixture ), sudden droplet heating during transition is observed for 


two cases: 
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(i) when the initial naass fraction of the low volatility fuel is very low ( e.g. 0.01 as ' 
shown in Fig. 5.7a ), the sudden droplet heating during transition occurs near the end of 
the burning, and it’s intensity being less, the rise in surface temperature is much below 
the boiling point of the less volatile fuel and it’s liquid side mass fraction remains much 
below 1.0 ( Fig. 5.8a ) and 

(ii) when it’s initial mass fraction is very high ( e.g. 0.9 ), the sudden droplet heating 
during transition takes place at the very beginning of burning and subsequently, it’s inten- 
sity being quite high, the surface temperature rises up pretty close to the boiling point of 
the low volatile fuel. For such a case, it’s liquid side surface mass fraction reaches a value 
very close to 1.0 ( Fig. 5.8a ). 

But for an intermediate value e.g. for the case of 0.1 initial mass fraction of the low 
volatile fuel which reveals multicomponent burning feature more distinctly, the droplet 
heating takes place over a longer period of time and the temperature rise lies in between the 
above mentioned two cases. Shaw and Williams ( 1990 ) also reported the sudden droplet 
heating to be present in their analysis for an n-decane droplet undergoing combustion 
with an initial low volatile impurity mass fraction of 0.1 or 0.02, but, the variation in 
temperature rise for different initial mass fraction combinations is not exhibited by the 
model of Shaw and WilHams ( 1990 ). They showed that for both cases, the liquid side 
surface mass fraction of the low volatility fuel reaches almost unity after certain time. This 
means that their model does not take into account the effect of fuel supply from droplet 
interior to the surface. Fig. 5.8a shows that the highest liquid side mass fraction reached by 
the low volatility fuel varies depending on it’s initial mass fraction and it remains around 
0.82 and never reaches unity for the 0.1 initial mass fraction case. This happens because of 
the high droplet radius regression in the later part of droplet burning for which the more 
volatile fuel is constantly supplied to the droplet surface. Fig. 5.13a and 5,13b supports 
the above reasoning because Fig. 5,13a shows that when the effect of surface regression is 



neglected, the liquid side surface mass fraction of the low volatile fuel reaches almost unity ■ 
and correspondingly, the surface temperature (Fig. 5.14a) reaches the boiling point of the 
low volatile fuel earlier or later, depending on it’s initial mass fraction. 

5.15 Conclusions 

In the second part of the thesis, a full numerical model has been developed for the case of a 
transient multicomponent single droplet burning in an otherwise quiescent atmosphere. An 
adaptive grid generation technique has been used based on the equidistribution of a positive 
weight function such that the local discretization error is minimised. This adaptive grid 
generation technique helps to record accurately the flame structure as well as the position 
and movement of the flame. This study gives a critical analysis of how a diffusion flame 
around the liquid droplet evolves during the ignition period and then how it moves away 
from the droplet, and finally, how it is extinguished. The study reveals that the flame 
zone is initially established very near the droplet surface and then as the chemical reaction 
gathers momentum, the flame temperature rises very fast and almost instantly reaches it s 
peak temperature. Almost simultaneously, the flame expands towards the droplet surface 
insearch of fuel vapour for sustaining the enhanced chemical reaction at the flame. As 
the flame moves and expands towards the droplet surface, the temperature gradient at the 
surface becomes very high as a result of which the droplet surface is heated up very fast and 
the surface temperature approaches the boiling point of the more volatile fuel. So, the fuel 
starts evaporating vigorously. For this high evaporation rate, the flame starts moving away 
from the droplet surface due to fuel vapour accumulation. But as the burning proceeds, 
the droplet volume becomes progressively lesser for which the overall mass-gasification rate 
falls gradually. Therefore, fuel vapour accumulation between the droplet surface and the 
flame becomes progressively lesser. At about mid way of burning, the fuel evaporation 
rate at the surface almost equals the consumption rate at the flame and the flame becomes 
almost steady at one place. It observed that the flame remains steady, maintaining a 



fixed distance from the droplet center till the end of droplet burning. The flame exists 
till the droplet becomes very small. The total burning time for a heptane-hexadecane 
droplet for 0. 9-0.1 initial mass fraction combination is about 65.625 milli seconds. Ignition 
becomes complete within 3 milli seconds and extinction takes about 31.9 ^is. Present model 
identifies six distinct periods for the entire burning process. They are pre-ignition, ignition, 
first <P — law period, transition, second — law period and extinction. Results show that 
the pre-ignition period, ignition period, first (P - law period, 2"“^ transient ( transition ), 
second P - law period and extinction covers approximately 2.28 %, 2.28 %, 35.6 %, 19.5 
%, 40 % and 0.36 % respectively of the total burning time. These six zones or periods 
depend on three aspects: volatility differential of the species, their initial mass fraction 
combination and the prevailing ambient conditions. Liquid side surface mass fractions of 
the species are observed to be the direct causes that control the two P — law periods and 
their transition in between. From the results, it is seen that it is very difficult to point out 
one single parameter or the root cause which controls the distinct features of the 2 (P — law 
periods during multicomponent fuel droplet combustion. Though apparantly, the liquid 
side mass fraction combination of the fuels seems to the controlling factor. It is seen from 
the present study that during extinction, the droplet radius regresses very fast. As a result 
of this, the transformed non-dimentional coordinate system also shrinks very sharply. On 
the other hand, the flame position remains quite steady. Due to the shrinking coordinate 
system, the far stream boundary conditions comes very close to the flame, which does 
not happen in practical situations. Thus, error can creep in while studying the extinction 
phenomenon. 

The study reveals that the multicomponent feature of burning is predominant when 
the volatility differential of the two components is high, and also, when the initial liquid 
side mass fraction of the more volatile fuel is high ( around 0.9 ). Therefore, volatility 
differential and initial mass fraction combination are the two key parameters for studying 
multicomponent fuel droplet combustion. For the present study, liquid side Peclet number 



is found to be around 150 which is very high. Therefore, the assumption of diifusion limit 
model is justified. This study brings out the two (P — law periods quite explicitly. The 
second (P — law period is identified with the ’quasi-steady’ mode of burning during which 
the droplet composition and temperature remain almost constant. It is also seen that this 
period covers almost 95 % of the total burning time for heptane-dodecane burning with 
0. 1-0.9 initial mass fraction combination. This value decreases to about 27 % for the 0.9- 
0.1 combination case. This value increases when the boiling point difference of the two 
species is more and vice-versa. Results show that the flame is not thick. Comparison of 
the three analytical models ( see Figs. 5.13 through 5.15 ) show that the surface regression 
rate is neither too high nor too low. Therefore, one can neither neglect the effect of 
surface regression, nor can we assume the quasi steady mode of burning to be valid from 
the very beginning. Droplet radius regression can be neglected only in the initial part of 
the droplet burning. But, later on, as the radius regression rate gathers momentum, the 
quasi steady mode of burning becomes more and more established. Gas phase dissociation, 
variable gas phase specific heat consideration and fuel vapour accumulation are important 
considerations and they should not be neglected. Through parametric study of important 
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parameters such as D, D, Cp, Cp and ambient temperature, it is found that the total 
burning time as well as other burning characteristics are very much sensitive to the property 
values being used. Comparing with the quasi-steady model of chapter 3, it is observed 
that the quasi-steady assumption of the gas phase is quite reasonable barring the initial 
transient and the final extinction. But the effect of surface regression can not be neglected. 
The present result shows very good agreement with the experimental results of Okazima 
and Kumagai ( 1975 ) for hepane fuel burning under zero gravity condition with normal 
ambient conditions. When the present model is compared with the experimental results of 
Wang et al. ( 1984 ) for the burning of heptane-hexadecane fuel with 0. 7-0.3 initial liquid 
side mass fraction combination, the compaxison is quite satisfactory considering the fact 
that the present model deals with spherically symmetric combustion whereas, Wang et al. 



( 1984 ) performed their experiment on freely falling droplets under natural convection 
conditions. Their experiment showed three periods, the two periods correspond to the 
preferential gasification of the two fuels and the third period corresponds to the short 
intermediate transition in between during which vigorous droplet heating takes place. The 
present model also exhibits three dominant zones, two cP — law periods separated by one 
transition zone. But, this transition in case of Wang et al. ( 1984 ) is reported to be very 
sudden as against a gradual change or transition exhibited by the present model. The 
present model does not include the effect of thermal expansion for the liquid phase which 
may cause this deviation. 

5.16 Contributions to Knowledge 

i. The entire burning period is identified with six distinct zones. 

ii. Duration and nature of the — law periods ( hence, the burning rate constants 
) and the transition in between are observed to be directly dependent upon the gas side 
mass fraction of the fuels and indirectly dependent upon the initial mass fraction as well 
as the volatility differential of the fuels. 

iii. Transition from one <P — law period to the other is seen to be gradual. It becomes 
sharper when D high and also when the effect of surface regression is neglected. 

iv. Quasi steady mode of burning ( when the droplet composition hardly changes ) is 
found to correspond to the second — law period . 

v. Flame initially move towards the droplet surface, then it moves away from the 
surface and finally, at the end, it comes back and vanishes into the droplet surface. 

vi. Variable specific heat and gas phase dissociation lowers the flame temperature 
considerably. 

5.17 Recommendations for Further Study 

L Variable property inodel should be incorporated to get more realistic results. 



2. The present spherically symmetric model can be relaxed to include the effect of flow 
around the droplet. This can be done in two ways: 

(i) The full Navier-Stoke’s equations can be used for solving both the liquid and the gas 
phase flow situations. This will give us the complete picture regarding the the importance 
of both batch distillation limit and diffusion limit models. 

(ii) A known low Reynolds number flow solution for the gas phase can be first assumed 
and then the present model can be extended to include the effect of internal circulation 
through spherical vortex model and see how this affects the droplet burning. 

This will give us a clear picture regarding the relative importance of batch distillation 
and diffusion limit models. 

3. The present adaptive grid generation technique automatically generates the grid spacing, 
but can not choose the time interval for the numerical solution on it’s own. Error due to 
improper choice of the time interval selection can be included as an additional part of the 
adaptive grid generation process. 

4. radiation loss is important and it can be included to get more realistic results. 

5. micro-explosions can be looked into. 

6. The effect of the shrinking coordinate on the study towards the end of droplet burning 
is quite significant and it demands special attention. This can be looked into for more 
accurate analysis. 

7. Effect of convective terms arising in both the liquid and the gas phases on different 
parameters can be looked into in more details. 
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/alues of Constants and Properties; 


i/alue of universal gas constant (R ) = 8.134 J/Mol-K. 

u 

\/alue of constants used in Arrhenius equation (2.76) are: 

11 4 

h = 4X10 ; a = 0.0; E = 12.6 X 10 J/mol . 



heptane 

decane 

undecane 

dodecan 

Latent heat of 
evaporation 
in J/kg 

3 . 182,- 

X 10^ 

2.8 

X 10^ 

2.6 

X 10^ 

2.55 _ 

X 10 

Density in 
kg/m^ 

684 

730 

740 

748 

Boiling pt. 
in K 

371 

447 

469 

489 

Molecular Wt . 

100 

142 

156 

170 


For the gas phase : 

3 

Density ip) = 1.1614 kg/m . 

Binary diffusion coefficient (D) = 2.25 x 10 m /s. 

Specific heat (C^) ® 1.007 x lO^J/kg-K. 

Thermal conductivity (X) = 26.3 x 10 ^ W/m-K. 

For the liquid phase: 

^ ”“9 2 

Binary diffusion coefficient (D) = 5.0 x 10 m /s. 

Specific heat (C ) = 1.75 x lO^J/kg-K. 

Thermal conductivity (X) =0.1 W/m-K. 

For the variable specific heat case, the data were taken from 
Penner (1957). Other references consulted for the property values 
are Reid et al. (1988) and Shaw (1990) . 





